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INTERNATIONAL JOURNAL OF ABSTRACTS 
STATISTICAL THEORY AND METHOD 


COVERAGE OF JOURNAL 


THE aim of this journal of abstracts is to give complete coverage of papers in the field of statistical 
theory and new contributions to statistical method. Papers which only report applications or examples 


| of existing statistical theory and method will not be included. There are approximately two hundred 


and fifty journals published in various parte of the world which are wholly or partly devoted to 
the field of statistical theory and method and which will be brought within the scope of this journal of 
abstracts. In due course it will be possible to issue a complete list of these journals. In the case of the 


_ following journals, however, being those which are wholly devoted to statistical theory—all contributions, 


- 


whether a paper, note or miscellanea, will be abstracted : 
Annals of Mathematical Statistics, 
Biometrika 
Journal, Royal Statistical Society (Series B) 
Bulletin of Mathematical Statistics 
Annals, Institute of Statistical Mathematics 


_ Within the larger group of journals, which are not wholly devoted to statistical theory and method, there 


are some journals which have the vast majority of their contributions in this field. These journals, there- 
fore, will be abstracted on a virtually complete basis : 

Biometrics 

Metrika 

Metron 

Review, International Statistical Institute 

Technometrics 

Sankhya ~ 
After experience in the publication of this journal of abstracts it may be found desirable to add further 
journals to these lists. In any case readers may be assured that all papers properly to be included in this 
journal of abstracts will be included irrespective of such notification. If any reader of this journal discovers 
a paper which happens to have been overlooked, the General Editor will be pleased to be informed so 
that the appropriate abstract can be made: always provided that the date of publication is after rst 


- October 1958, when the abstracting for this journal commenced. 


In addition to the ordinary journals, there are two kinds of publication which fall within the scope 
of this journal of abstracts. They are the experiment and other research station reports—which occur 
particularly in the North American region—and the reports of conferences, symposia and seminars. 


"Whilst these latter may be included in the book review sections of journals it is unusual for any individual 


contribution to be noted at any length. These publications are, in effect, special collections of papers 
and for this reason the appropriate arrangements will be made for them to be included in this journal. 


By the same token, abstracts of papers given at conferences and reproduced in an appropriate journal will 


be disregarded until the definitive publication is available. 


FORMAT OF JOURNAL 


The abstracts will be about 400 words long—the recommendation of UNESCO for the “long” 


’ abstract service: they will be in the English language although the language of the original paper will 
be indicated on the abstract together with the name of the abstractor. In addition, the address of the 


author(s) will be given in sufficient detail to facilitate contact in order to obtain further detail or request 
nually in a supplement. 
lil 


The scheme of classification has been developed upon lines that will facilitate the transfer to punched 
cards of the code numbers allocated to each abstract: to allow for future development it is suggested that use 
is made of 4-column fields. Each abstract will have two classification numbers: the primary number in 
heavy type to indicate the basic topic of the paper and the secondary number in brackets to take account of 
the most important cross-reference. The sheets of the journal are colour-coded according to the twelve main 
sections of the classification and it should be noted that it is the main section number of the primary 
classification which determines the colour code for each abstract. It is believed that this method of 
colour-coding the pages will provide a distinctive visual aid to the identification of abstracts both when 
the journal is in bound form or dismantled and filed on cards or in binders. The format and simple 
binding allows of the following alternative treatments by users of the journal : 

(a) Leave intact as a shelf-periodical. 

(b) Split and filed in page form according to the main sections of the classification. 

(c) Split and guillotine (single cut) each page ready for : 

(i) pasting onto standard index cards. 
(ii) filing in loose-leaf or other binders—for which the appropriate holes are punched. 


It will also be possible for those users who need a completely referenced file to insert skeleton cards or 
sheets according to the secondary classification number provided on each abstract. 


ae 

Q Primary and 

secondary 

1 classification 

! 

Holes for filing 

into binders 

| 

| 

i 

L 

GC) Serial number 

ik st abstract 
[es] 
feet Guillotine guide marks — — — — — — | 
stale 
Primary and 


secondary 
classification 


= 36) 


Holes for filing 
into binders 
| 


| 
1 
v 
O Serial number 


ays of abstract 


peaeerest| 


INTERNATIONAL JOURNAL OF ABSTRACTS 
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LIST OF PUBLICATIONS FOR REGULAR SCANNING 


The following list of publications in the statistical and allied fields will be regularly examined 
for papers which should be abstracted for this Journal. Other journals will also be scrutinised 
from time to time and a check maintained by scanning other relevant abstracting journals. 
In addition, abstracts will be included of special collections of papers as published in reports 
of conferences, symposia and seminars together with the published reports of experiment and 


other research stations. 


Acta Acad. Cient., Lima 
Acta Math., Stockh. 
Acta Math. Acad. Sci. Hung. 
Acta Math. Sinica 
Acta Pontif. Acad. Sci. 
Advanced Management 
Advance. Enzymol. 
Advanc. Protochem. 
Advance. Vet. Sci. 
Agric. Engng., St Foseph, Mich. 
Agriculture, Lond. 
Agron. fF. 
Allg. Statist. Arch. 
Amer. Econ. Rev. 
Amer. Heart F. 
Amer. F. Clin. Nutr. 
Amer. }. Clin. Path. 
Amer. }. Digestive Diseases 
Amer. }. Hum. Genet. 
Amer. F. Hyg. 
Amer. F. Med. 
Amer. F. Physiol. 
Amer. F. Publ. Health 
Amer. F. Sci. 
Amer. F. Sociol. 
Amer. Ff. Vet. Res. 
Amer. Nat. 
_ Amer. Statistician 

Analyt. Chem. 
An. Soc. Cient. Argent. 
Ann. Acad. Sci. Fenn. 
Ann. Appl. Biol. 
_ Ann. Fac. Econ. Com. Catania 


Ann. Fac. Econ. Com. Palermo 


Ann. Fac. Sci. Porto 

Ann. Hum. Genet. 

Ann. Inst. Poincaré 

Ann. Inst. Statist. Bari 

Ann. Inst. Statist. Math., Tokyo 

Ann. Inst. Sup. Sci. Econ. 
Financ. 

’ Ann. Math. Statist. 

Ann. N.Y. Acad. Sct. 

. Ann. Sci. Ec. Norm. Sup., Paris 

Ann. Sci. Econ. Appl. 

Ann. Soc. Sci. Bruxelles, I 


publication 


Acta de la Academia Cientifica, Lima 

Acta Mathematica 

Acta Mathematica Academiae Scientiarum Hungaricae 

Acta Mathematica Sinica 

Acta Pontificiae Academiae Scientiarum 

Advanced Management 

Advances in Enzymology and Related Subjects 

Advances in Protochemistry 

Advances in Veterinary Science 

Agricultural Engineering 

Agriculture (Min. of Agric.) 

Agronomy Journal 

Allgemeines Statistisches Archiv 

American Economic Review 

American Heart Journal 

American Journal of Clinical Nutrition 

American Journal of Clinical Pathology 

American Journal of Digestive Diseases 

American Journal of Human Genetics 

American Journal of Hygiene 

American Journal of Medicine 

American Journal of Physiology 

American Journal of Public Health 

American Journal of Science 

American Journal of Sociology 

American Journal of Veterinary Research 

American Naturalist 

American Statistician 

Analytical Chemistry 

Anales de la Sociedad Cientifica Argentina 

Annales Academiae Scientiarum Fennicae 

Annals of Applied Biology 

Annali della Facolta di Economia e Commercio, 
Catania (University) 

Annali della Facolta di Economia e Commercio, 
Palermo (University) 

Annaes da Faculdade de Ciéncias do Porto 

Annals of Human Genetics (London) 

Annales de I’Institut Henri Poincaré 

Annali dell’Instituto di Statistica, Bari (University) 

Annals of the Institute of Statistical Mathematics 

Annaes do Instituto Superior de Ciéncias Economicas 
e Financieras 

Annals of Mathematical Statistics 

Annals of the New York Academy of Science 

Annales Sciéntifiques de l’Ecole Normale Supérieure 

Annales de Sciences Economique Appliquées 

Annales de la Société Scientifique de Bruxelles— 
Serie 1 

Vv 


The coverage of this Journal amounts to at least 400 sources of 


Peru 
Sweden 


Argentina 
Finland 
Great Britain 
Italy 


Italy 


Portugal 
Great Britain 
France 

Italy 

Japan 
Portugal 


U.S.A. 
WES-AG 
France 
France 
Belgium 


Ann. Statist., Rome 


Ann. Univ. Lyon, A 


Ann. Univ. Sci. Budapest. (Sect. 


Math.) 
Annu. Rev. Physiol. 
Aplik. Mat. 
Appl. Sci. Res., Hague 
Appl. Statist. 
Archimede 
Arch. Biochem. Biophys. 
Arch. Math., Karlsruhe 
Ark. Mat. 
Atti (mem.) Accad. Lincei 


Atti Riun. Sct. Soc. Ital. Statist. 


Aust. F. Agric. Res. 
Aust. J. Appl. Sct. 
Aust. F. Biol. Sci. 
Aust. F. Phys. 
Aust. F. Sct. 

Aust. ¥. Statist. 


Beil. Mber. Ost. Inst. 
Wirtschaftsf. 

Bell Lab. Record 

Bell Syst. Monographs 

Bell Syst. Tech. F. 

Biometrics 

Biometrika 

Biom. Zeit. 

Bla. Dtsch. Ges. Versich.-math. 


Bol. Fac. Filos. Ciént., S. Paulo 


Bol. Inst. Actuar. Port. 
Bol. Soc. Mat. Mexicana 
Boll. Cent. Ric. Operat. 
Boll. Un. Mat. Ital. 

Brit. FJ. Nutr. 

Brit. F. Prev. Soc. Med. 
Brit. ¥. Sociol. 

Brit. F. Statist. Psychol. 
Bull. Acad. Sci., 111 

Bull. Amer. Soc. Test. Mat. 
Bull. Ass. Actuatr. Suisses 
Bull. Calcutta Math. Soc. 
Bull. Calcutta Statist. Ass. 
Bull., C.S.I.R.O., Aust. 


Bull. Indian Soc. Qual. Contr. 
Bull. Indian Stand. Inst. 

Bull. Inst. Egypte 

Bull. Jap. Statist. Soc. 

Bull. Math. Statist. 


Bull. Nat. Inst. Sci. India 
Bull. Oxf. Univ. Inst. Statist. 
Bull. Res. Inst. Kerala 

Bull. Sci. Math. 

Bull. Soc. Math. Belg. 


Annali de Statistica (Instituto Centrale di Statistica, 
Rome) 

Annales de l’ Université de Lyon—Section A 

Annales Universitatis Scientiarum Budapestiensis de 
Rolando Eétvés nominatae—Sectio Mathematica 

Annual Review of Physiology 

Aplikace Matematiky 

Applied Scientific Research, The Hague 

Applied Statistics 

Archimede 

Archives of Biochemistry and Biophysics 

Archiv der Mathematik 

Arkiv for Matematik 

Atti della Accademia Nazionale dei Lincei 

Atti della Riunione Scientifica Societa Italiana di 
Statistica, Rome 

Australian Journal of Agricultural Research 

Australian Journal of Applied Science 

Australian Journal of Biological Sciences 

Australian Journal of Physics 

Australian Journal of Science 

Australian Journal of Statistics 


Beilagen zu den Monatsberichten des Osterreichischen 
Institutes fiir Wirtschaftsforschung 

Beil Laboratories Record 

Bell System Monographs 

Bell System Technical Journal 

Biometrics 

Biometrika 

Biometrische Zeitschrift 

Blatter der Deutschen Gesellschaft fiir Versicherungs- 
mathematik 

Boletim da Faculdade de Filosofia, Ciéncias e Letras, 
Sao Paulo 

Boletim do Instituto dos Actuarios Portugueses 

Boletim de la Sociedad Matematica Mexicana 

Bollettino del Centro per la Ricerca Operativa 

Bollettino dell’ Unione Matematica Italiana 

British Journal of Nutrition 

British Journal of Preventive and Social Medicine 

British Journal of Sociology 

British Journal of Statistical Psychology 

Bulletin de l’Académie des Sciences (Classe 3) 

Bulletin, American Society for Testing Materials 

Bulletin de l’Association des Actuaires Suisses 

Bulletin of the Calcutta Mathematical Society 

Bulletin, Calcutta Statistical Association 

Bulletins, Commonwealth Scientific and Industrial 
Research Organisation 

Bulletin, Indian Society of Quality Control 

Bulletin, Indian Standards Institution 

Bulletin de l'Institut d’Egypte, Cairo 

Bulleting of the Japan Statistical Society 

Bulletin of Mathematical Statistics, Research Associa- 
tion of Statistical Sciences 

Bulletin, National Institute of Sciences of India 

Bulletin of the Oxford University Institute of Statistics 


Bulletin, Central Research Institute, Kerala (University) 


Bulletin des Sciences Mathématique 
Bulletin de la Société Mathématique de Belgique, 
Gembloux 
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Italy 


France 
Hungary 


U.S.A. 


Czechoslovakia 


Netherlands 
Great Britain 
Italy 

Wis. A? 
Germany 
Sweden 

Italy 

Italy 


Australia 
Australia 
Australia 
Australia 
Australia 
Australia 


Austria 


U.S.A. 
U.S.A. 
U.S.A. 
U.S.A. 

Great Britain 
Germany 
Germany 


Brazil 


Portugal 
Mexico 

Italy 

Italy 

Great Britain 
Great Britain 
Great Britain 
Great Britain 
Poland 
U.S.A. 
Switzerland 
India 

India 
Australia 


India 
India 
Egypt 
Japan 
Japan 


India 

Great Britain 
India 

France 
Belgium 


Bull. Soc. Math. France 
Bull. Statist. Soc., N.S.W. 


Canad. }. Econ. Polit. Sci. 
Cancer, New York 

Cas. pést. math. 

Chem. Engng.. Progr. 

Clin. Chem. 

Colect. Estud. 

Coll. Mat. 

Collog. Math. 

Comment. Math. Helvetia 
Comment. Pontif. Acad. Sci. 
Compos. Math., Groningen 
C.R. Acad. Sci., Paris 
Cornell Vet. 

Current Science 

Curso Mira Fernandes 
Czech. Math. 7. 


Dacca Univ. Statist. Bull. 

Defense Sci. F. 

Dokl. Akad. Nauk. SSSR 

Dopov. Akad. Nauk. Ukrain, 
SSR 


Econ. Appl. 

Econ. Internatzionale 
Econ. f. 

Econ. Leban. Arabe 
Econ. Record 
Econometrica 
Economica 

Edinb. Math. Notes 
Egypte Contemp. 
Egyptian Statist. 7. 
Ekon. Tidskr. 

Emp. #. Exp. Agric. 
Endocrinology 
Engineer, Lond. 
Engineering, Lond. 
Ergebn. Naturw. Grenzgeb. 


Esta 

Estadistica 
Euclides, Groningen 
Eugen. Rev. 
Evolution 


Food Res. 
Food Tech. 


Ganita 

Gaz. Mat. 

G. Economisti 

G. Ist. Ital. Attuari 
Genetics 


Genus 
Growth 


Harvard Busin. Rev. 
Helvetica Phys. Acta 


Bulletin de la Société Mathématique de France 
Bulletin of the Statistical Society, N.S. Wales 


Canadian Journal of Economic and Political Science 
Cancer (New York) 

Casopis pro péstovani mathematiky 

Chemical Engineering Progress 

Clinical Chemistry 

Colectanea de Estudos 

Collectunea Matematica 

Colloquium Mathematicum 

Commentarii Mathematica Helvetia 
Commentationes Pontificiae Academiae Scientiarum 
Compositio Mathematica, Groningen 

Comptes Rendus de |’Académie des Sciences, Paris 
Cornell Veterinarian 

Current Science 

Curso Mira Fernandes 

Czechoslowak Mathematical Journal 


Dacca University Statistical Bulletin 
Defense Science Journal 

Doklady Akademii Nauk, SSSR. 

Dopovedy Akademii Nauk Ukrainskoj, SSR 


Economie Appliquée 

Economia Internatzionale 

Economic Journal 

L’Economie Lebanaise & Arabe, Beirut 
Economic Record 

Econometrica 

Economica 


' Edinburgh Mathematical Notes 


L’Egypte Contemporaine, Cairo 

Egyptian Statistical Journal 
konomisk Tidskrifts 

Empire Journal of Experimental Agriculture 

Endocrinology 

Engineer 

Engineering 

Ergebnisse der Naturwissenschaften und ihrer Grenz- 
gebiete : 

Esta 

Estadistica 

Euclides, Groningen 

Eugenics Review 

Evolution 


Food Research 
Food Technology 


Ganita, Lucknow-Allahabad (Universities) 
Gazeta de Matematica 

Giornale degli Economisti e Annali di Economia 
Giornale dell’Istituto Italiano degli Attuari 
Genetics 

Genus 

Growth 


Harvard Business Review 
Helvetica Physica Acta 
vii 


France 
Australia 


Canada 

USA. 
Czechoslovakia 
Use 

U.S.A. 

Spain 

Spain 

Poland 
Switzerland 
Italy 
Netherlands 
France 

U.S.A. 

India 

Spain 
Czechoslovakia 


Pakistan 
India 
W.SsS2Re 
SSR. 


France 

Italy 

Great Britain 
Lebanon 
Australia 
U.S.A. 

Great Britain 
Great Britain 
Egypt 

Egypt 
Sweden 
Great Britain 
WUsSz4e 

Great Britain 
Great Britain 
Germany 


Spain 

WEseAe 
Netherlands 
Great Britain 


U.S.A. 


U.S.A. 
U.S.A. 


India 
Spain 
Italy 
Italy 
U.S.A. 
Italy 
U.S.A. 


WS: As 
Switzerland 


Heredity 
Hum. Biol. 


IFO-Studien 
Inc. Statist. 
Indag. Math. 


Indian Coun. Agric. Res. Bull. 
Indian Coun. Agric. Res. Statist. 


Newsletter 
Indian Econ. f. 
Indian Econ. Rev. 
Indian F. Agric. Econ. 
Indian }. Phys. 


Indian #. Power & River Devel. 


Indian F. Soc. Work 
Indian Math. Soc. 7. 
Industr. Engng. Chem. 
Industr. Qual. Contr. 
Industria, Milano 

Towa Acad. Sct. F. 

Towa Exp. Sta. Publ. 
Izv. Akad. Nauk, SSSR 


Izv. Akad. Nauk, Uzbek., SSR. 


3b. Nat.-Okon. Statist. 
Jber. Dtsch. Mathver. 
F. Agric. Food Chem. 
F. Agric. Sct. 

F. Amer. Pharm. Ass. 


J. Amer. Statist. Ass. 
JF. Amer. Vet. Med. Ass. 
FJ. Animal Sci. 

JF. Appl. Physiol. 

JF. Ass. Comp. Mach. 

F. Aust. Inst. Agric. Sct. 


F. Aust. Math. Soc. 
F. Biol. Chem. 

F. Chronic Diseases 
F. Clin. Invest. 

F. Dairy Sci. 

F. Ecology 

J. Econ. Ent. 

F. Exp. Educ. 

F. Exp. Med. 

F. Forestry 

F. Gen. Physiol. 

F. Genet. 

F. Heredity 

F. Hyg., Camb. 

F. Industr. Econ. 
J. Industr. Engng. 


F. Indian Soc. Agric. Statist. 


F. Inst. Actuar. 

F. Inst. Actuar. Stud. Soc. 
JF. Madras Uuiv., B 

JF. Mammalogy 

F. Marine Res. 

F. Math. Pures Appl. 

F. Math. Soc. Japan 

F. Nutrition 

F. Operat. Res. 


Heredity 
Human Biology 


IFO-Studien 

Incorporated Statistician 

Indagationes Mathematicae, Amsterdam 

Indian Council of Agricultural Research: Bulletin 

Indian Council of Agricultural Research: Statistical 
Newsletter 

Indian Economic Journal 

Indian Economic Review 

Indian Journal of Agricultural Economics 

Indian Journal of Physics 

Indian Journal of Power and River Development 

Indian Journal of Social Work 

Indian Mathematical Society Journal 

Industrial and Engineering Chemistry 

Industrial Quality Control 

L’ Industria 

Iowa Academy of Science Journal 

Iowa State College Experiment Station Publications 

Izvestiya Akademii Nauk, SSSR 

Izvestiya Akademii Nauk, Uzbekskoj, SSR. 


Jahrbucher fiir Nationalékonomie und Statistik 

Jahresbericht der Deutschen Mathematikervereinigung 

Journal of Agricultural and Food Chemistry 

Journal of Agricultural Science 

Journal of the American Pharmaceutical Association 
(Scientific Edition) 

Journal, American Statistical Association 

Journal of American Veterinary Medicine Association 

Journal of Animal Sciences 

Journal of Applied Physiology 

Journal of the Association of Computing Machinery 

Journal of the Australian Institute of Agricultural 
Science 

Journal of the Australian Mathematical Society 

Journal of Biological Chemistry 

Journal of Chronic Diseases 

Journal of Clinical Investigation 

Journal of Dairy Science 

Journal of Ecology 

Journal of Economic Entomology 

Journal of Experimental Education 

Journal of Experimental Medicine 

Journal of Forestry 

Journal of General Physiology 

Journal of Genetics 

Journal of Heredity 

Journal of Hygiene, Cambridge 

Journal of Industrial Economics 

Journal of Industrial Engineering 

Journal of the Indian Society of Agricultural Statistics 

Journal of the Institute of Actuaries 

Journal of the Institute of Actuaries Students’ Society 

Journal of Madras University, Section B. 

Journal of Mammalogy 

Journal of Marine Research 

Journal de Mathématiques Pures et Appliquées 

Journal of the Mathematical Society of Japan 

Journal of Nutrition 

Journal of the Operations Research Society of America 

Vili 


Great Britain 
UIS. A. 


Austria 
Great Britain 
Netherlands 
India 

India 


India 
India 
India 
India 
India 
India 
India 
USIAS 
WESrA.. 
Italy 
U.S.A. 
U.S.A. 
WESsSaEe 
U.S:S.R. 


Germany 


. Germany 


Wss-A; 
Great Britain 
WESEAG 


U.S.A. 
WiS:As 
U.S.A. 
Wes Ae 
U.S.A. 
Australia 


Australia 
U.S.A. 
U.S.A. 
U.S.A. 
W.S-As 

Great Britain 
WeSzAG 
WSsAG 
U.S.A. 

WES AG 
WESEA 

India 

WS eA. 

Great Britain 
Great Britain 
WiSvAe 

India 

Great Britain 
Great Britain 
Great Britain 
WES vAs 

WES sAe 
France 
Japan 
U.S.A. 
U.S.A. 


F. Pharmacol. 


F. Plant Physiol. 

F. Polit. Econ. 

Jj. reine angew. Math. 

F. Res., Nat. Bur. Stand. 


JF. Roy. Soc., N.S.W. 
F. R. Statist. Soc. (A & BY 


F. Sct. Engng. 

F. Sci. Industr. Res. 

F. Sct. Industr. Res., India 

F. Soc. Statist. Paris 

F. Statist. Social Inquiry Soc. 


J. Univ. Baroda 


Kodai (Tokyo Inst. Tech.) 


Magy. Tud. Akad. III Oszl. 


Kozl. 
Management Sci. 
Manager (India) 
Manchester School 
Matematica 


Mat. Lapok. 

Math. Ann. 

Math. Gaz. 

Math. Nachr. 

Math. Stud. 

Math. Tab., Wash. 
Math., Tech., Wirtschaft 


Math. Zeit. 

Meded. Landbhoogesch., 
Wageningen 

Mem. Coll. Sci., Kyoto 

Mem. Fac. Sct., A, Kyusyu 


Mém. Soc. Sci. Liége 

Metabolism 

Met. Mag., Lond. 

Metrika 

Metroeconomica 

Metron 

Mitt. Ver. Schweiz. Versich.- 
Math. 

Monatsh. Math. 

Monthly Not. R. Astr. Soc. 


Nat. Bur. Stand., Statist. Rep. 


& Publ. 
Nature 
Ned. Verzek. Actuar. By. 


Nieuw Arch. Wisk. 
Nieuw Tijdschr. Wisk. 


J ournal of Pharmacology and Experimental Thera- 
peutics 

Journal of Plant Physiology 

Journal of Political Economy 

Journal fiir reine und angewandte Mathematik 

Journal of Research, National Bureau of Standards, 
Washington 

Journal of the Royal Society, N.S. Wales 

Journal of the Royal Statistical Society (Series A and 
B) 

Journal of Science and Engineering 

Journal of Scientific and Industrial Research 

Journal of Scientific and Industrial Research 

Journal de la Société Statistique de Paris 

Journal of the Statistical and Social Inquiry Society of 
Ireland 

Journal, Maharaja Sayajiroa, Baroda University 


Kodai (Mathematical Seminar reports) Tokyo Institute 
of Technology 


Magyar Tudomanyos Akadémia III Matematika és 
Fizikai Oszlalayanak Kozleményei 

Management Science 

Manager 

Manchester School 

Matematica, Notas del Instituto de Matematica 
Facultad de Ingenieria de Rosario 

Matematikai Lapok 

Mathematische Annalen 

Mathematical Gazette 

Mathematische Nachrichten 


‘ Mathematics Students 


Mathematical Tables and other aids to Computation 
Mathematik, Technik, Wirtschaft (Zeitschrift fur 
moderna Rechentechnik und Automation) 
Mathematische Zeitschrift 

Mededeelingen van de Landbouwhoogeschool, 
Wageningen 

Memoirs of the College of Science, Kyoto (University) 
Memoirs of the Faculty of Science (Series A), Kyusyu 
(University) 

Mémoires de la Société Royale des Sciences de Liége 
Metabolism 

Meteorological Magazine 

Metrika 

Metroeconomica 

Metron 

Mitteilungen der Vereinigung Schweizerisher Ver- 
sicherungs-Mathematiker 

Monatsheft fiir Mathematik 

Monthly Notices of the Royal Astronomical Society 


National Bureau of Standards, Washington : Statistical 
Reports and Publications 

Nature (London) 

Nederlands Verzekeringsarchief, Actuarieel Bijvoegsel, 
The Hague 

Nieuw Archief voor Wiskunde, Groningen 

Nieuw Tijdschrift voor Wiskunde, Groningen 
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U.S.A. 
Germany 
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Australia 
Great Britain 


India 
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France 
Treland 


India 


Japan 


Hungary 


WS AN 

India 

Great Britain 
Argentine 


Hungary 
Germany 
Great Britain 
Germany 
India 

U.S.A. 


Germany 


Germany 
Netherlands 


Japan 
Japan 


Belgium 
WSs 

Great Britain 
Germany 
Italy 

Italy 
Switzerland 


Germany 
Great Britain 
U.S.A. 


Great Britain 
Netherlands 


Netherlands 
Netherlands 


Operat. Res. Quart. 
Osaka Math. fF. 
Ost. Ingen.-arch. 
Oxf. Econ. Papers 


Pakistan Econ. f. 

Phil. Mag. 

Phil. Phenom. Res. 

Phil. Trans. 

Population Stud. 

Portug. Math. 

Poultry Sct. 

Proc. Amer. Soc. Hort. Sct. 


Proc. Amer. Soc. Qual. Contr. 
Proc. Amer. Soc. Test. Mat. 


Proc. Camb. Phil. Soc. 
Proc. Indian Acad. Sct. 
Proc. Indian Ass. Cult. Sci. 


Proc. Indian Sci. Congr. 

Proc. Inst. Statist. Math. 

Proc. fapan Acad. 

Proc. Kon. Ned. Akad., 
Wetensch. A 

Proc. Lahore Phil. Soc. 

Proc. Lond. Math. Soc. 

Proc. Math. Phys. Soc. Egypt 


Proc. Nat. Acad. Sct. India 
Proc. Nat. Inst. Sct. 
India (A & B) 
Proc. Pakistan Statist. Ass. 
Proc. Qual. Contr. Ass., 
Bangalore 
Proc. R. Irish Acad. 
Proc. Roy. Soc. A & B 
Proc. Roy. Soc. N.S.W. 
Proc. Roy. Soc. N.Z. 
Proc. Soc. Exp. Biol., N.Y. 


Proc. Soil Sci. Soc. Amer. 
Przegl. Statyst. 
Psychometrika 

Public Health Rep. 

Public Opinion Quart. 
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ABRHAM J. (Charles University, Prague) 


An approximate method for non-linear programming—ZIn Czech 


Cas. pest. math. (1958) 83, 425-439 (5 refs.) 


The paper gives a solution of the problem of finding the 
minimum of a strictly convex function on the set of all 
non-negative solutions of a system of linear equations. This 
is a generalisation of the basic problem of linear programming. 

The author constructs an iterative process for solving 
the problem and shows that it converges to the true solution. 
He shows further how this method can be applied to quadratic 
programming. He discusses a special case where the func- 
tion to be minimised is a sum of a linear form and a positive 
definite quadratic form and shows that in this case the 
iterative process can often be modified so that it stops after 
a finite number of steps. In the last section of the paper the 
author gives a numerical example with some practical advice 
for getting the solution more quickly. 

The iterative process formulated by the author is the 
following: let us have a system of m-+k linear equations 
in m unknown quantities, let M/ be the set of all non-negative 
solutions of this system. By a vector of the set M shall be 
meant a solution of the corresponding homogeneous system. 
For every point X in M there can be chosen a vector basis 
(uy, ...) Ux) Satisfying the following conditions: (a) the 
vectors u; are linearly independent, (6) any u; has zeros 
at all or at all but one components which are zero for X, 
(c) any component which is zero for X is non-zero for exactly 
one u;. 


AKAIKE, H. (Institute of Statistical Mathematics, Tokyo) 


0.8 (0.1) 


Let f be the strictly convex function we have to mini- 
mise on M. We start with an arbitrary point X, of WM and 
choose a corresponding basis (u;, ..., u,) and then a point 
X, = X,+tu, such that f(X,) be minimum under the 
condition that X,<¢M. If the given basis is a basis for Xo, 
too, we seek X,; = X,+tu,. in the same manner, and we 
continue so for X4, ..., taking the vectors u; cyclically one 
after the other. If the given basis cannot be used for the 
new point, a new basis is chosen and we proceed as when 
starting with Xj. 


(F. Zitek) 


0.6 (6.4) 


On a computation method for eigenvalue problems and its application to statistical 


analysis—In English 


Ann. Inst. Statist. Math., Tokyo (1958) 10, 1-20 (8 refs., 2 tables) 


The paper deals with the numerical method for solving the 
N-dimensional eigenvalue problem Ax = ABx: where A 
is a real symmetric matrix and B a positive definite real 
symmetric matrix, both N by N, and x is an N-dimensional 
column vector (N being large). The paper commences with 
a discussion of a general approximation procedure followed 
by individual cases: first for one vector (k = 1) and then 
for two specified forms of the functions é. A computation 
scheme is presented which is cyclic and converging with 
checks for the errors of rounding. In the actual computa- 
tion a finite digital representation of each quantity has to be 
made and because of this the problem sometimes reduces 
to a linear form after some steps in the approximation. The 
paper proceeds to discuss an accelerated scheme with par- 
ticular regard to the convergence property and rounding 
errors. The theoretical exposition of the papers concludes 
with a generalisation of the scheme—after noting the con- 
vergence to the eigenvector with the maximum eigenvalue, 
and that in practical applications it is frequently necessary 
only to determine this single eigenvector. 

The scheme of computation given in this paper for a 
high-speed computer and the numerical examples given in 
the paper were all worked out on an automatic relay com- 
puter and the technical details will be the subject of a 
further paper. The paper concludes with some statistical 


applications in the field of multivariate analysis : canonical 
correlations, the linear discriminant function, the quantifica- 
tion problem and trend analysis of multiple time series. 
For canonical correlation the ordinary formulation as well 
as the least squares method is briefly considered. 


(C. Hayashi) 
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COX, C. P. (Nat. Inst. for Res. in Dairying, Univ. of Reading) 0.6 (6.1) 
A concise derivation of general orthogonal polynomials—In English od 
JF. R. Statist. Soc. B (1958) 20, 406-407 


Orthogonal polynomials for describing the regression of y 
and x are easily derived from first principles. The author 
gives an alternative derivation by generalising the poly- 
nomials from the general regression equation in its deter- 
minantal form using the Schweinsian expansion of a 
determinant. The partitioning of a sum of squares into 
independent parts corresponding to successively fitted 
orthogonal polynomials is also derived in a similar way. 


(F. N. David) 


FURST, D. (University of Rome) | 0.7 (0.8) 
An application of the Theory of Games—In Italian 
Archimede (1958) 10, 211-223 (4 refs., 4 figs.) 


The author gives the principles and the basic methods of 
the theory of games in an elementary form. Starting from a 
particular example (a three by three game), he introduces 
the general problem of two-person zero-sum games and, 
after showing its links with linear programming, describes 
the method of solving it. The third section is devoted to 
the solution and discussion of an example. After solving 
the game, the author studies how the value of the game 
and the optimum strategy vary when the pay-off matrix 
changes according to given rules. 

The problem is also treated from a geometrical point of 
view, and this gives to the author the opportunity of pointing 
out the geometrical interpretation of the fundamental 
theorem of Von Neumann. 

The last part of the paper is devoted to studying the 
case in which a player can modify, within certain restraints, 
the pay-off matrix. Since the example examined in the 
paper is of this type, the author treats the problem from 
the point of view of that player, i.e. finding the pay-off 
matrix for which the value of the game is a maximum. 


(G. Dall’Aglio) 
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GOOD, I. J. (Royal Naval Scientific Service, Cheltenham) 0.5 (9.2) 
The interaction algorithm and practical fourier analysis—In English 
J. R. Statist. Soc. B (1958) 20, 361-372 (9 refs.) 


The interactions of a 2” factorial experiment can be calcu- 
lated using a simple algorithm suggested by Yates. Box 
suggested an extension of the algorithm of Yates for the 3” 
experiment. The author points out that these algorithms 
imply the existence of linear formulae for the interactions 
and proves a theorem whereby these formulae can be written 
down. 

He considers the situation where the interactions are 
known and it is desired to “‘ recover” the observations. 
An inverse algorithm is given for this reverse process. The 
uses of this inverse algorithm are discussed. Apart from 
the obvious use that it enables the forward calculation of 
the interactions to be checked the author suggests that 
non-significant high order interactions may be put equal to 
zero and the inverse algorithm used to smooth. the observa- 
tions with the effect of making them closer to those expected 
on hypothesis. 


(F. N. David) 


KABE, D. G. (Karnatak University, Dharwar, India) fea 0.4 (3.9) 
Applications of Meijer-G functions to distribution problems in statistics—In English 
Biometrika (1958) 45, 578-580 (5 refs.) 


The author points out that, by contrast with the extensive 
use of the Fourier transformation in statistical distribution 
theory, the Mellin transformation appears to have been 
neglected. He notes that the Mellin transformation appears 
particularly useful in those cases where the moments of a 
statistic can be expressed in terms of the product and the 
ratio of I-functions. As illustration he derives the distri- 
bution of the generalised correlation ratio and the joint 
distribution of two variables, each a function of sample 
canonical coefficients, given in 1939 by Girskick [Ann. 
Math. Statist., 10, 203]. 


(F. N. David) 
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LEVI, E, (Inst. of Financial Mathematics, University of Catania) 
Fitting of life tables by sums of exponential functions—ZJn Italian 


Statistica (1959) 19, 3-19 (15 refs.) 


The problem involving the theory of life table functions in 
connection with values of life annuities by means of the values 
of annuities-certain has recently been treated by several 
authors. The present method establishes a necessary and 
sufficient condition under which the l,, function is expressed 
as the general solution of an homogeneous linear differential 
equation of the kth order with constant coefficients. The 
above solution is given by 


k-1 
l= 2 buf 
r=0 
where the 6, are constants and the u, are real or complex : 
the cases for which the u, are all real are known. 

Some mathematicians—cited by the author—are 
interested in suggesting particular applications of the expres- 
sion given above which show some convenient methods 
for determining the parameters. To overcome difficulties 
concerning the method of least squares, the author follows 
two procedures : one based on ‘“‘ points”’ and another based 
on “‘sums’”’, that is to say “interpolation by points’ and 
“interpolation by sums ’’. 

By analogy with the Steffensen system of calculating 
life annuities, it is possible to unite the two systems in one. 
In the above unified system of equations the determinant of 
coefficients is the well-known Hankel determinant. In order 


VARADARAJAN, V. S. (Indian Statistical Institute, Calcutta) 


A remark on strong measurability—In English 
Sankhyd (1958), 20, 219-220 (3 refs.) 


A mapping of a measure-space into a metric space is called 
measurable if the inverse image of every Borel set of the 
metric space is a measurable set. A measurable mapping is 
called strongly measurable if it is the almost everywhere 
limit of a sequence of measurable mappings each one of 
which assumes only a finite number of values. It is well 
known that if the metric space is separable, then every 
measurable mapping into it is strongly measurable. In this 
paper the author discusses this proposition when no separ- 
ability assumptions are made on the metric space. Assuming 
the continuum hypothesis or, more generally, assuming 
only that the power of the class of all subsets of any un- 
countable set is greater than that of the continuum, the 
author demonstrates that every measurable mapping of a 
measure-space into a metric space is strongly measurable 
provided that the o-field in the measure-space 1s generated 


by denumerably many sets. 


(R. Ranga Rao) 


0.3 (0.2) 


to assist the solution, the author suggests a useful trans- 
formation of this determinant. Furthermore, he suggests 
that the parameters of the above expression should be fitted 
using two iterations of the procedure. 

After considerable discussion concerning the roots of 
the characteristic equation, the author gives an example 
concerning the interpolation of functions of the above type 
on the life table published in the November 1957 issue of 
Bolletino Mensile di Statistica, Rome. This set of demo- 
graphic data comes from the Italian Census of Population, 
1951, and from the mortality of the male population during 
the years 1950-1953. 


(V. Amato) 
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BARTON, D. E. & DAVID, F. N. (Uni ity C 
Bune ie Rite tn ae (University College, London) 1.3 (1.4) 


Biometrika (1958) 45, 572-578 (5 refs., 6 tables) 


It is supposed that there are r identical beads of & colours 
r;, being the number of beads of the ith colour. PHase 
beads are supposed to be threaded on a string to form a 
bracelet which may not be turned over although it may be 
rotated. The authors point out that two different problems 
are involved. First, the beads may be imagined threaded 
on a string and the string then tied to form a ring and 
second, the beads may be imagined placed in a circle in which 
symmetries have to be allowed for. Classical writers such 
as Whitworth discuss the first problem and assume that it 
is the second. 

The first problem of the line bent to the circle follows 
previous work on the alternations of colour in a line. Tables 
are given showing the exact number of alternatives for two, 
three, and four colours and up to twelve beads. The 
moments of these distributions are noted, as is also the fact 
that the positive binomial! with the correct first two moments 
is a suitable approximating function. : 

The alternations in the second problem are given the 
name of Jablonski runs after the mathematician who, among 
others, gave a partial solution for the total number of 
arrangements of this kind. Two different methods of 
enumerating these alternations are given and a further set 
of tables has been calculated for those cases in which the 
number of Jablonski runs is different from those of the first 


problem. 
(F. N. David) 
9 
BASU, D. (Indian Statistical Institute, Calcutta) I.0 (4.1) 


On statistics independent of sufficient statistics—In English 
Sankhyda (1958) 20, 223-226 (2 refs.) 


The author in an earlier paper [Sankhya 15, 377] stated 
that any statistics independent of a sufficient statistic must 
have the same distribution for all values of the unknown 
parameter. In this paper the author gives an example to 
show that this proposition is not true in general. Con- 
ditions under which this is true are discussed. For this, 
the author needs the concept of “‘ connectedness ” between a 
pair of parameter points which he introduces as follows : 
8 ¢ Sand 6’ « Q are said to be connected (by a statistic T) 
if there exists a finite number of parameter points 
0,, 95, ..., 0, which along with @ and 6’ form a chain : 

Bo Origen oy Ons Ones (putting 0 = 4) and 6’ = 4,1) 
such that the distributions of JT’ overlap for any two con- 
secutive members of this chain. ‘Two distributions are said 
to overlap if there is some interval to which they assign 
probabilities less than unity. 

The author shows that if all pairs of parameter points 
in 8 are connected by the sufficient statistic, then any 
statistic independent of the sufficient statistic has the same 
distribution for all values of the unknown parameter. 


(S. K. Mitra) 
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BAYES, T. 


Essay towards solving a problem in the doctrine of chances—In English 


Biometrika (1958) 45, 296-315 
[Reproduced from Phil. Trans. (1736) 53, 370-418] 


This is a reproduction of the famous paper of the Rev. 
Thomas Bayes, F.R.S. It is in the form of a letter from 
Richard Price to John Canton, and begins, ‘‘I now send 
you an essay which I have found among the papers of our 
deceased friend Mr. Bayes, and which, in my opinion, has 
great merit, and will deserve to be preserved”. Price 
explains what he thinks Bayes has done but is not entirely 
explicit about his own additions and subtractions. Bayes 
begins with “ Problem: Given. The number of times in 
which an unknown event has happened and failed: Re- 
quired. The chance that the probability of it happening in 
a single trial lies somewhere between the two degrees of 
probability that can be named.” He defines his terms and 
in Section I gives a series of seven propositions regarding 
the addition and multiplication of probabilities 

Section II is devoted to the solution of the Problem. 
Balls are thrown randomly on a square table and the 
probability of their resting place is discussed. The argu- 
ment leads to three rules. “‘ Rule 1. If nothing is known 
concerning an event but that it has happened p times and 
failed q, in p+q or n trials, and from hence I guess that 
the probability of its happening lies somewhere between 
any two degrees of probability as X and x, the chance that 
I am right in my guess is...” Rules 2 and 3 are variants 
of Rule 1. ‘“‘ The first rule gives a direct and perfect solu- 


II 


BRAUMANN, P. (Faculty of Sciences, Lisbon) 


Generalisation of two theorems for multidecision statistical games—In Spanish 
Rev. Fac. Cienc. Lisbon, 2a serie, A-Ciéndias Matemdaticas (1959) 7 


The author deals with the problem of extending some results 
contained in Blackwell & Girshick’s Theory of Games and 
Statistical Decisions. 

He defines a symmetric multidecision problem as one 
with a finite number of decisions in which the integration 
in the sample space and the decision space may be inter- 
changed. The paper then discusses the result that, for this 
kind of problem, whatever is the strategy of N ature for each 
second kind of strategy of the statistician their exists a first 
kind of strategy with the same (equivalent) risk. ay 

Finally the author characterises some decision functions 
with minimum a-posteriori risk as invariant Bayes strategies 
in some subspace of the strategies. 


(J. Tiago de Oliveria) 
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1.0 (11.9) | 


tion in all cases: and the two following rules are only particular | 
methods approximating to the solution given in the first rule.” | 

In the Appendix applications of the rule to particular | 
problems are given. The probability of the Sun rising, a 
problem supposed to originate with Laplace, is discussed | 
among others. 


(F. N. David) 
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ERDOS, P., RENYI, A. & SZUSZ, P. (Eétvis University, Budapest) 


On Engel’s and Sylvester’s series—In Hungarian 


1.5 (10.1) 


Ann. Univ. Sci. Budapest—Sect. Math. (1958) 7-32 (16 refs.) 


The authors are dealing with metrical properties of Engel’s 
series of the first kind and Sylvester’s series (also called 
Engel’s series of the second kind). Let the Engel’s and 
Sylvester’s series of the real number X (o< X<1) be 


X= 1/4 +1/qige+...+1/q192---dn+-.- (1) 
and 

X = 1/Q,+1/Q,+...... +1/On+...... (2) 
respectively. Here g, and Q, are defined as follows: 


Qn41 = {1/Tny1} where ty =x, Tayr = Ty(r,) and T(t) 
= t{1/t}—1 (o<t<1) and Q,,, = {1/R,}, where Ry = &, 
Ri4i = T(Rn) and T,(t) = t—1/{1/t}. Here {z} denotes 
the least integer which is =z. ‘The following theorems are 
proved: (1) the probability that the digit k occurs 
exactly r times in the sequence gq, is equal to k—1/k’—}. 
Here, and in what follows, the interval (o, 1) and the Lebesgue 
measure is considered as the space of elementary events 
and the probability measure respectively. 


lim Pog an Vi<y) = (2/29) |" edt (2) 


n> +0 


n 
| lim Vq, = z) = y (3) 
n> +0 
P( lim log q,—n/V 2n log log n = r) = 1 (4) 
n+>+0 
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-FRANCKX, E. (Ecole Royale Militaire, Bruxelles) 


jim P Uos{(Qn/Q102---On-1) —n}/Vn}<y] 
= GIVan)[ ettay (5) 


lim log Q,/2” (6) 


N> +o 
exists and is finite and positive for almost all X. 

Theorems 2, 3 and 4 have been found earlier by P. 
Lévy and E. Borel [E. Borel, Sur les développements uni- 
taires normaux, Comptes Rendus (1947) 225, 51, E. Borel, 
Sur les développements unitaires normaux, Ann. Soc. Polo- 
naise de Math. (1948) 21, 74-79, P. Lévy, Remarques sur 
un théoreme de M. Emile Borel, Comptes Rendus (1947) 
225, 918-919]. The authors give new and detailed proofs 
for these known results—as well as for their new results. 
The proofs of the theorems on Engel’s series are based on 
the fact that the sequence of random variables q,(x) forms a 
homogeneous Markoff chain and the random variables 
log Gn/n—1 are almost independent. 


(P. Révész) 


1.6 (-.-) 


Strong Law of Large Numbers for uniformly limited variables—In Spanish 


Trab. Estadist. (1958) 9, 111-115 


The strong law of large numbers holds for a sequence X; of 
' aleatory variables if the difference between the arithmetic 
mean 5S, of the first n variables and the mathematical expecta- 
tion of this mean converges to zero with probability = 1, 
when n goes to infinity. The characteristic sub-set is defined 
as a sub-sequence Sy, such that both the Lim w,4,/u; = 1 


and convergence to zero hold when the centres of the vari- 
ables are at the mathematical expectation. 

If any assumptions on the dependence or independence 
of the variables are made, considering only bounded variables, 
the following criterion is established: the strong law of 
large numbers holds for the uniformly limited variables if 
the sequence S, has a characteristic sub-set. The variables 
X; are uniformly limited if the absolute value of the differ- 
ence between X; and its mathematical expectation m, 1S 


smaller than a certain number L which is independent of 7. 


It can be demonstrated that this condition is necessary by 
the fact that the set in which S,, does not converge 18 con- 
tained in the set in which S, does not converge. é 

It is not so easy to show that the condition is sufficient. 
For this the author enunciates three previous lemmas : 
Lemma A. If »» is a sub-sequence of positive numbers 
uniformly limited by L, and m; is an increasing sequence of 
positive integers such that Lim ui4,/u; = 1 then Sn—S,p F 
becomes as small as possible choosing ¢ sufficiently large. 


S,, denotes the arithmetic mean of the n first numbers of 
the sequence A, (nm must satisfy the condition n;=n<n,444). 
Lemma B. A similar condition is established when the 
numbers X, are replaced by uniformly limited aleatory 
variables X,—mgn. Lemma C. It is proved by means of 
lemma B that it is equivalent to asserting that the event 
which states that the variable S,,—m,, is greater than the 
absolute value of « for some k>n is incompatible with the 
event which expresses the absolute value of ¢/z for alln;>n, 
always provided that n be sufficiently large. Hence, it is 
easy to assume that the existence of the characteristic sub_ 
set is sufficient. In fact, from the condition that the pro_ 
bability of convergence of Sng—MSy, be equal to 1, it 


follows that the probability of divergence of S,—ms, must 


be equal to zero. 
hold. 


So the strong law of large numbers must 


(P. Zoroa) 
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MOTOO, M. (Institute of Statistical Mathematics, Tokyo) 1.5 (10.1) 
Proof of the Law of Iterated Logarithm through diffusion equation—In English 
Ann. Inst. Statist. Math., Tokyo (1958) 10, 21-28 (5 refs.) 


The aim of this paper is to prove by a unified method a 
group of theorems of the iterated logarithm type for the 
absolute value of the k-dimensional Brownian motion 
(Rass %. 250.) 

The author first considers a group of four classes (C) of 
increasing (decreasing) functions in the interval (0, ©) of 
a diffusion process {X(z)}. If the process is recurrent and 
the expectation of its recurrence time is finite—as in the 
Uhlenbeck process—these four classes of function are com- 
pletely determined. One theorem and three lemmas are 
given in this section of the paper. The next section applies 
the theorem to the k-dimensional Brownian motion—a one- ) 
dimensional stationary diffusion process which does not 
satisfy the above condition of recurrence. By using a change | 
of scale and of time it is shown how the theorem of the | 
previous section can be rendered applicable to this case. | 
‘Two theorems on the law of the iterated logarithm are proved | 
and the paper concludes with a brief application to a general | 
process subject to certain regularity conditions. 


(C. Hayashi) 
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NASR, S. K. (Faculty of Science-Moharem Bey, Math. Dept., Alexandria) 1.5 (1.6) | 
A law of large numbers for abstract random variables—In English | 
Metrika (1958) 1, 89-98 (15 refs.) | 


- The author considers random variables with values in an Z, is a sufficient estimation of the expectation within a | 
abstract set X in which a separating family D of pseudo- family of normally distributed random variables with 
metrics is given. He recalls and discusses his definition, common variance. Finally, some particular cases are 
based on the ideas of S. Doss [Bull. Sci. Math. (1949) (2), considered. 
73, 48-72], of the expectation of a random variable and the 
arithmetic mean Z, of random variables X,, ..., Xn. The 
sequence Z,, is called strongly stable relative to D if there 
is a sequence c, in X such that d(Zp, Cn) converges stochastic- 
ally to zero for every din D. This stability is said to be 
normal if, in addition, the expectation of every Z, exists 
and can be chosen for c,. If X denotes a Banach space and 
D the family of semi-norms defined by the elements of the 
dual space of X, the strong normal stability of Z, for any 
sequence X,, of independent random variables with a common 
distribution follows immediately from the classical weak law 
of large numbers. In comparing this result with the one 
of Fortet and Mourier [C.R. Acad. Sci., Paris (1953), 231, 

18-20] the author remarks that the latter impose restrictions 
on X and the X, ; however, their conclusions are, of course, 
much stronger, viz. almost certain strong convergence instead 
of stochastic weak convergence. Next the author combines 
his result with the one of Doss [Publ. Inst. Statist., 
Paris (1954), 3, 135-142] asserting that in any Banach space 


(K. Krickeberg) 
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OKAMOTO, M. (Osaka University) 


Some inequalities relating to the partial sum of Binomial probabilities—In English 


I.0 (2.4) 


Ann. Inst. Statist. Math., Tokyo (1958) 10, 29-35 (5 refs.) 


This paper gives a simplified proof for a somewhat stronger 
inequality than that given by J. V. Uspensky, at page 102 
of his book, Introduction to Mathematical Probability (New 
York, 1937) relating to the partial sum of binomial pro- 
babilities. If x is a random variable following a binomial 
distribution B(n, p), Uspensky states that 


P{| x/n—p | 2c} <2 exp (—ne*/2) 
where c>o and o<p<t. 


After considering two preliminary lemmas, four theorems 
are stated of the form P(«/n—p=c) < exp (—anc?) (Theorem 
1(i)). Theorem 1(ii) reverses the inequality within the 
bracket and states ¢ with a negative sign. Theorems 2(i) 
and (ii) modify the exponent to exp (—nc?/2pq) where 
p2o-'5 and pSo's for the negative value of c. Theorems 
3 and 4 take the root of x/n and p in the basic inequality 
with the exponent unchanged in Theorem 3 but halved in 
Theorem 4, which deals with the negative value of c in a 
similar way to Theorems 1(ii) and 2(ii). 

The paper concludes with an application of Theorems 
3 and 4 to Matusita’s multinomial distance—a test of fit 
for an empirical multinomial distribution [Ann. Math. 
Statist., 26, 641-640 and Ann. Inst. Statist. Math. 7, 67-80]. 


(C. Hayashi) 
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RENYI, A. (Math. Inst., Hungarian Academy of Science, Budapest) 


1.4 (11.7) 


On a one-dimensional problem concerning random place filling—Jn Hungarian 
Mag. Tud. Akad. Kut. Mat. Intézet. Kézleményet (1958) 109-127 (7 refs.) 


The following problem is treated : how many unit intervals 
can be placed on the interval (0, x) in the following manner : 
the first unit interval is placed at random on the interval 
(0, x) (this means that its centre J is uniformly distributed 
in the interval (4, «—4); similarly, asecond unit interval is 
chosen at random (in the same sense) independently from the 
first. If the second interval intersects the first it is discarded 
and the choice repeated till the chosen interval does not 
intersect the first one. The process is repeated in the same 
way, i.e. if there have been chosen already n (disjoint) 
unit intervals, we choose another interval at random, but 
retain it only if it does not intersect any of the previous 
intervals; otherwise the choice is repeated, etc. The 
process is finished when there is no possibility of placing a 
further unit interval without intersecting any of the previous 
intervals. ‘The number v, of the unit intervals thus obtained 
is a random variable ; this is studied in this paper. 

The author derives recursive formulae for the moments 
and for the distribution function of vz. Introducing Laplace- 
transforms, and applying a Tauberian theorem, he deduces 
that, denoting by M(v,) the mean value of v,, 


es t 
[e (af 1 —e~“/n)du~o.748, 


ro) ° 


lim M(v,)/« = C 


x co 


- and that v,/x tends stochastically to C for x>oo. 


The solution of the problem has the following practical 
application, remarked by N. G. de Bruijn: there is a parking 
place consisting of a pavement of L metres length. Every 
car arriving places itself randomly in the intervals being still 
free (it does not strive to close up). ‘The problem is how 
many cars can park if each one requires a free interval of 
length 1 ? It follows that with a high probability the number 
of these cars will be approximately 0.748 L/I. 


(K. Sarkadi) 
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RENYI, A. (Math. Inst., Hungarian Academy of Science, Budapest) 


Investigation by means of probability theory of the conductivity 


of certain resistances—In Hungarian 


Mag. Tud. Akad. Kut. Mat. Intézet. Kézleményei (1958) 247-256 (3 refs., 1 fig.) 


Account is given of some probabilistic considerations con- 
cerning the conductivity of resistances which are manu- 
factured by pressing the mixture of carbon-powder with the 
powder of some non-conducting material (e.g. phenoplast). 
Such resistances are, as it is well known, used in radio- 
technics. The propagation of the current in the resistance 
mentioned above is considered as a chain reaction. The 
calculations have been carried out under the supposition 
that the carbon-powder is distributed with uniform density 
(according to Poisson’s law) in the whole body of the resist- 
ances. ‘The numerical results were in striking contradiction 
with the experimental facts. This has led Mr J. Katona 
(of the Industrial Research Institute for Telecommunica- 
tion) to investigate the microstructure of the mentioned 
resistances. ‘The microscopic pictures (one of which is 
reproduced in the paper) have thrown light on the reason 
of the contradiction between theory and facts, as they have 
shown that the carbon-powder is not uniformly distributed, 
but is contained in thin layers. This observation gave a 
possibility to improve the quality of the mentioned resist- 
ances. 


(Author’s summary) 
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THIONET, P. (Inst. Nat. Statist., et Etudes Econ., Paris) 


I. (10.0) 


On the average interval between two consecutive Friday’s 13th—In French 


J. Soc. Statist. Paris (1958) 99, 284-293 


This paper is an amusing study posing (without resolving ; 
the problem is unsolvable) the question of applying the 
calculus of probabilities to the study of an isolated event. 
The author refers to the central point of Borel’s views on this 
subject—that each case should be considered on its merits 
and without preconceived ideas. 

After references to groups of equally possible events and 
that the observations are not always a realisation of a stoch- 
astic process, the author states that the case to be considered 
is not a chain process or especially ergodic. He reassures 
the reader as to all this mathematical vocabulary : the paper 
is concerned with the appearance of Friday, 13th on the 
calendar. It is not the approach of a new year, and the 
arrival of calendars, that raises this problem afresh but a 
_paper given by Nemchinov at the ISI Conference in Brussels 
(September 1958). The author makes use of an application 
of arithmetically ordered probabilities and the following 
sentence summarises his point of view : 


“Tt seems to us that, in certain cases, events deter- 
mined rigorously but by a sufficiently complex calcu- 
lation so as to reveal the nature of the simple prob- 
ability consideration, permitted at least a partial explana- 
tion.” 


The paper then deals with the statistical regularities 
observed in the data and the author concludes that the 


sequence of events includes (qa) statistical regularities which 
can be explained by a simple probability rule, (6) statistical 
regularities which are unexplainable and (c) a tendency to 
irregular oscillations which also need explaining. After a 
section on the problem of a causal explanation, the author 
concludes with a critical commentary on the explanation 
based upon probability considerations. 


(D. Dugué) 
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VARADARAJAN, V. S. (Indian Statistical Institute, Calcutta) 


A useful convergence theorem—ZJn English 
Sankhya (1958) 20, 221-222 (2 refs.) 


This paper discusses the convergences of a sequence X,, 
of k-dimensional vectors. It is well known that in order that 
this Sequence may converge to a given vector X in distri- 
bution, it is necessary and sufficient that for every linear 
function I, the sequence of random variables (Xp) converge 
in distribution to the random variable W(X). In this paper 
the author considers the case when the existence of Y is 
not known beforehand, and proves that in order that there 
may exist a vector X towards which the X,, are convergent 
in distribution, it is necessary and sufficient that the distri- 
butions of the random variables 1(X,,) converge to some 
distribution for every linear function J. As a corollary, the 
author deduces that in order that a sequence of vectors be 
asymptotically normal, it is necessary and sufficient that 
every linear function of the vectors be asymptotically normal. 


(R. Ranga Rao) 
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ZOROA TEROL, P. (Instituto de Investigaciones Estadisticas, Madrid) 


1.5 (-.-) 


1.8 (10.1) 


A problem on sequential decision in economics—In Spanish 


Trab. Estadist. (1958) 9, 103-110 


A trader has to sell some goods in one of n successive 
opportunities that can be called years 1, 2,...”. ‘The profit 
he can get in the year 7 is a random variable represented by 
x;: the joint distribution of (x, ..., x,) is known by assump- 
tion. ‘The trader has to choose the year for selling, con- 
sidering the profit he may get for this year and the precedent 
ones. Let a be the profit obtained when m non-selling years 
had passed. Then the sale must.be performed if this year 
is reached without selling and x,5=a. So the variable x» 
is changed into another one X,, which represents the best 
profit the trader can get if the year m is reached without 
selling. Its value is x» when x,2a and a when x»<a. 

The best profit X,,_, obtained if the year m—1 is reached 
without selling can be deduced from the knowledge of Xm. 
Its value is either x_, when xm_; is bigger than the profit 
expected for the year m or the value expected for the year 
m when xm_1 is smaller than the profit expected if the sale is 
realised during the year m. In the first case the sale must be 
performed in the year m—1 and in the second one it would 
be better to sell in the next year. Following on this argu- 
ment a recurrent law may be found in order to establish the 
best strategy for selling. Another similar situation is stated 
‘assuming that the last year in which the sale is possible 
is not known for sure, but this case will be ruled by a random 


principle. 
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The first problem is generalised to the case of a con- 
tinuum of possibilities, i.e. that a customer may arrive at 
any instant in an interval of time o<t<1. The arrival of 
the customer obeys a Poisson’s process. If there is a possi- 
bility of performing the sale then the distribution of the 
variable x; which represents the profit at this moment is 
known. The condition satisfied by the solution of this 
problem is found and is applied to a simple case. 


4 


(S. Rios) 
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CROW, E. L. 


(Nat. Bureau of Standards, Boulder, California) 


2.5 (2.1) 


The mean deviation of the Poisson-distribution—In English 


Biometrika (1958) 45, 556-559 (4 refs., 1 fig.) 


After deriving the mean deviation of the Poisson distribution, 
the author considers the behaviour of the ratio R(m) of the 
mean deviation to the standard deviation as the mean m 
increases from zero. R is shown to follow a continuous 
damped oscillation with infinitely many maxima and minima, 
the minima occurring at integral values of m and the maxima 
being equally spaced between the minima. The slope of 
R is always finite but is discontinuous at integral values of m. 

The value of R for a normal distribution is (2/7)? and 
this is attained by the Poisson R once and only once between 
each adjacent pair of maximum and minimum values. It is 
shown that R is a poor criterion for distinguishing between 
the Poisson and normal distributions, or for determining 
whether a Poisson distribution is approximated by the 
normal. 

An asymptotic expression for R(m) is obtained and this 
gives the pair of points between successive minima at which 
R takes the normal value. The accuracy of the asymptotic 
expression is discussed with reference to the maximum and 
minimum values it gives for R(m) as m varies, and also with 
reference to the asymptotic ratio of the maximum positive 
difference from the normal R, to the average of the absolute 
values of the adjacent maximum negative differences. 

The ratio R for the binomial distribution is also con- 
sidered, using an asymptotic formula given by Johnson 
(Biometrika, 1957) in his note on the mean deviation of 
the binomial distribution. 
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The behaviour of R(np) is similar to that of the Poisson 
R(m), and when np is varied by varying p with n fixed, 
analagous results are obtained approximately. Since these 
are confirmed for small n by exact calculations, it is con- 
cluded in particular that for any fixed there are exactly 2n 
values of p for which the binomial R takes the normal value. 

For fixed p and varying n, R is not continuous but its 
behaviour can be investigated by regarding n as continuous. 
Some points of this R are plotted in a diagram for two 
values of p. This diagram also shows the behaviour of the 
binomial ratio for two fixed values of n, and the Poisson 
and normal ratios for m between o and 10. ‘The upper 
bound of R for any distribution is attained with p, n having 
values 4 and 1 respectively. 


(D. Allibone) 


2.6 (3.8) 


Approximate formulae for the statistical distributions of extreme values—In English 
Biometrika (1958) 45, 447-470 (8 refs., 1 table, 3 figs.) 


A sample of size is supposed with the observations ranked in 
order of magnitude, say x,, <%s, <..., <Xn. Approximate 
formulae are derived for the distribution of the mth value, 
i.e. the mth from the right. The limiting distributions of 
extreme values, given by Fisher-Tippett and Gumbel, 
are well known. The author generalises the extreme value 
techniques to get his approximate formulae. Three types of 
functional form are supposed for the initial distribution of 
any x,;, viz.: exponential, Cauchy, and finite range. ‘The 
basic conditions required to be imposed on these initial 
distributions are discussed. A formula is obtained for what 
is designated “‘ the excess ”’ i.e. the difference between unity 
and the cumulative distribution function of any x. 

The general form of the distribution of the mth largest x 
is then considered. Formulae are derived whereby the mode 
and the maximum value of this distribution can be deter- 
mined. The behaviour of the maximum value for varying 
n and m is studied and approximate formulae for the distri- 
bution of the mth largest x are deduced. The range of appli- 
cation of each formula is delineated. The limiting formulae 
for mth largest x for the three types of initial distribution 
are found and the well-known limiting extreme-value 
distributions are also achieved. 


(F. N. David) 
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FERRERI, C. (School of Statistics, Palermo) 
An important new model of dependence—In Italian 


2.5 (1.1) 


Ann. Fac. Econ. Com. Palermo (1958) 12, 221-238 (6 refs.) 


The author has constructed an important theoretical model 
founded on the simple assumption of a constant dependency 
of successes: the unfavourable outcomes are supposed 
equally probably but without effect. 'The mathematical 
properties of the model are explained by a series of drawings 
of black and white balls. The author forms the one-sided 
dependency by replacing 1+A white balls for every white 
ball that is drawn or only the drawn ball when the ball is 
black. He has found, for the case of m (the number of 
trials) very large and rare events with slight dependency, 
the probability distribution function P, of x successes in 
n trials 
P, = [e*/I(A/d)] [(1 -e-9)* P(xe-+A/d)/x 1] 

where d = log (c?/M) is the dependency parameter, and 
A = np, (pz being the probability of drawing a white ball 
at the nth trial when x successes have been obtained in the 
preceding m—1 trials). 

The new model is a generalisation of the Poisson, 
negative binomial and Polya distributions. The parameter 
A is related to the variance o? and to the mean M by the 
relation d = M*d/2—M. 


(F. Galantino) 
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JOHNSON, N. L. (University College, London) 
The mean deviation, with special-reference to samples 


2.7 (4.3) 


from a Pearson Type III population—IJn English 


Biometrika (1958) 45, 478-483 (1 ref., 6 tables) 


A random sample of nm independent observations from a 
Type III distribution with a single parameter, «, is, at first, 
considered. The expected value of the sample mean devia- 
tion, m, (written m in the text), is found after making certain 
linear transformations on the observed variables : for larger n, 
two approximations are suggested. The population value 
of the mean deviation is the limiting value of the expected 
value of the sample mean deviation as n>. A table is 
provided showing the exact values of the ratio &(m,/o), 
o (= v/a, for the Type III population) being the standard 
deviation of the distribution, for different n and « = 1(1) 
8, co. A comparison of the values for different « shows that, 
for n=15, a small error in « will not have any serious effect 
on the unbiased estimator, m,[o/@(m,)] of o. A second 
table, calculated from the second approximation—which is 
better than the first, shows the values of &(m,/c) thus calcu- 
lated are in good agreement with the exact values, especially 
for nZ 10. 

The relative efficiency of two estimators of o, m,[o/&(mn)] 
and 5S, the sample standard deviation is discussed in general 
for any parent population. This is achieved by making 
certain large sample approximations to the expectations 
and variances of m, and S; the adequacy of the approxi- 
mations to the variance is investigated for the particular 
case when the parent population is Normal. A comparison 
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of approximate coefficients of variation of m, and S leads to 
the conclusion that m,[o/&(m,)] will have a smaller standard 
error than S if 8, > —3-+407/[@(m,)]?, where B, is the second 
moment ratio of the distribution of the observed variable. 
Specific applications of the above result are made to Type 
III, VIL and II populations, each having a single parameter, 
denoted by «. Three tables, one for each Type, are pro- 
duced showing values of 8, and 407/[&(m,)]? —3, for different 
a. It is concluded that sample mean deviation may, to the 
degree of approximation used, provide an unbiased estimator 
of o, with a smaller standard error than sample standard 
deviation, even for moderately leptokurtic populations. A 
modified form of the above criterion is suggested for finite n. 


(B. K. Ghosh) 


KAMAT, A. R. (Fergusson College, Poona) 
Hypergeometric expansion for incomplete moments of the 
Bivariate Normal Distribution—ZJn English 
Sankhyd (1958) 20, 317-320 (2 refs., 2 tables) 


2.3 (2.1) 


This paper can be regarded as a continuation of the author’s 
previous paper in Biometrika (1953) 40, 20-34. In this new 
paper, the author expresses the incomplete moment [m, n| 
of the standardised trivariate normal distribution in terms 
of the well-known hypergeometric series. The incomplete 
moment [m, n] is defined to be the contribution to the 
expectation of xy" from the positive quadrant. This is a 
simple multiple of the m, nth moment of the distribution 
when the variables are confined to the positive quadrant. 
The expressions for the absolute moments are at once 
deduced from those of the incomplete moments. The 
author has tabulated the incomplete moment [m, n] for all 
integers m, n such that m+n<4 and p = —1 (o-r)+1. 


(J. Sethuraman) 
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KAMAT, A. R. (Fergusson College, Poona) 
Incomplete moments of the trivariate normal distribution—In English 
Sankhya (1958) 20, 321-322 (1 ref.) 


2:3) (221) 


In this paper, which may be regarded as a sequel to the 
author’s other paper in this issue of Sankhyad—and both 
being in the nature of a continuation of the author’s previous 
paper in Biometrika (1953) 40, 20-34, the author defines 
_ the incomplete moment [/, m, n] for the standardised tri- 
variate normal distribution as the contribution from the 
positive octant to the expectation of x’y"3". Exact expres- 
sions of [/, m, n], when J, m and n are integers such that 
l+m+nS4, are given without proof. For the method of 
derivation the reader is referred to the author’s paper in 
Biometrika (1953). 


(J. Sethuraman) 
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KUDO, A. (Kyushu Univ., Japan, & Indian Statist. Inst., Calcutta) 
On the distribution of the maximum value of an equally correlated 


sample from a normal population—ZIn English 
Sankhyd (1958) 20, 309-316 (5 refs., 2 tables) 


In this paper a series expansion in terms of successive 
derivatives of the standardised normal distribution function 
is given for the probability integral of a multivariate normal 
distribution in which each variable has the same mean zero, 
the same variance o” and the same correlation p with every 
other variable. For the distribution of the maximum of an 
equally correlated sample from a normal population, what is 
required is the above integral with all limits of integration 
equal. In this case the series expansion can be more simply 
expressed in terms of certain functions Ay,(x) which are 
defined by means of recurrence relations involving de- 
rivatives of the standardised normal distribution function. 
Ay,»\(x) is tabulated for N = 1(1)10, m = 0(2)8, x = 1°8(-1)3°5. 

In the final section the author mentions the following 


application: 4%, %,..., Xn,.N, are N,+N, observations 
from a normal population with unknown mean m and known 
variance o*. Let xjg be the maximum of x, ..., xy, and x 


the mean of x, ..., Xy,,N,- The author has shown elsewhere 
[Sankhya 17, 67] that the statistic (x,)g—%)/o is optimum for 
testing the significance of an outlying observation provided 
it is known a priori that there may be just one x; among 
%4, ..., *N, Which is from another normal population with the 
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PARSONS, P. A. (Dept. of Genetics & ‘Trinity College, Cambridge) 


2.3 (5.5) 


same variance but with a greater mean value. As 
3 = Xj—K (7 = 1, ..., Ny) are equally correlated, the previ- 
ously developed theory gives the necessary distribution. 


(S. John) 


2.4 (1.2) 


Equilibria in auto-tetraploids under natural selection for a 


simplified model of viabilities—In English 
Biometrics (1959) 15, 20-29 (2 refs., 4 tables) 


It is pointed out that the study of gene frequencies for 
disomic populations has an immediate analogy with that for 
tetrasomic populations when the latter is considered in terms 
of frequencies of the gametic genotypes. The resulting 
recurrence relations are established before and after selection 
when there are two alleles present. 

The author states that the ‘‘ determination of equilibrium 
points for five different viabilities is difficult’’ and the rest 
of the paper is concerned with the simplified model men- 
tioned in the title. Verbally, this model is characterised 
by the fact that selection is symmetric for the two alleles, 
i.e. is determined by the ‘‘ degree”’ of heterozygosity of 
the organism at that locus. By standardisation of the 
viability of the most heterozygous (G2:g2) organism to unity 
the number of parameters is thereby reduced from five to 
two. 

This situation is then examined under selfing for : 

(i) Equilibria 

(ii) Stability of these equilibria. 
The two trivial equilibria of homozygosity are noted : these 
are both stable if the viability of the homozygous (G;,) 
organism is greater than } and also than that of the (G,g) 
organism. 

There is always a symmetric equilibrium whose fre- 
quencies the author expresses in terms of the heterozygous 


gamete frequency. This frequency is tabulated for viabilities 
[o, + (X2) 256] of the heterozygous organisms. 

There may be a pair of asymmetric equilibria, the con- 
ditions are established and the equilibrium frequencies are 
tabulated (where existent) for the same viabilities as above. 

The method of determining stability is displayed, the 
results are tabulated and discussed briefly. 


(C. W. Clunies-Ross) 
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RAMASUBBAN, T. A. (London School of Economics) 


The mean difference and the mean deviation of some discontinuous distributions—In English 


Biometrika (1958) 45, 549-556 (1 ref., 5 tables) 


The author considers the following discontinuous distribu- 
tions: the negative and positive binomial, Poisson, the log- 
arithmic and the geometric. For these distributions he finds 
algebraic expressions for the mean difference, A. For the 
hypergeometric, the negative binomial, Poisson, the logarithmic 
and the geometric distributions he finds algebraic expres- 
sions for 5, the mean deviation about the mean. The author 
shows that the usual definitions of the mean difference and 
the mean deviation (about the mean) can be considered as 
particular cases of more generalised forms. Furthermore, 
associated with these more generalised measures are two 
ratios involving the variance of the distribution as a standard- 
ising factor: one of these ratios is an extension of Geary’s 
ratio [F. R. Statist. Soc. (1930), 93, 442]. No question of 
sampling is involved and a further paper is envisaged ex- 
tending this study. 


(F. N. David) 
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STUART, A. (London School of Economics) 


Equally correlated variables and the multinormal integral—In English 


G. R. Statist. Soc. B (1958) 20, 373-378 (8 refs.) 


The joint distribution of random variables all equally 
positively correlated and with correlation equal to p, can be 
constructed from a set of +1 random variables. ‘The 
following result is demonstrated. “‘ Any statement con- 
cerning the means of m equally correlated variables with 
arbitrary means, equal variances, o*, and correlations 
-p,(—1/(n—1)S p<1), is precisely equivalent to the same 
statement concerning uncorrelated variables with the same 
means and equal variances o%(1—p). Any statement con- 
cerning o? for the equally correlated variables is directly 
translatable into one concerning o?(1 — p) for the uncorrelated 
variables.” The result is applied (i) to obtain a ¢ test for 
equally correlated variables (ii) to write down an integral 
for the evaluation of the positive secant of the multivariate 
normal surface where the variables are equally correlated 
(iii) to remark on the fact that problems regarding random- 
isation sets can be reduced to problems regarding independent 
normal variables. Equally negatively correlated variables 
are also discussed. 


(F. N. David) 


TALACKO, J. (Marquette University, Milwaukee) 
A note about a family of Perks’ distributions—In English 
Sankhya (1958) 20, 323-328 (4 refs., x table) 


It is well known in the classical probability theory that there 

exist bell-shaped distributions other than the symmetric . 
distributions like Gauss’, Student’s, Cauchy, etc. This paper 

is concerned with a class of such distributions called ‘‘ Perks’ 

functions’. The original Perks’ function is a ratio of two 

exponential polynomials. The case when the density 

function is of the special type c/(e*+k-+e-*) where c and k 

are parameters, is discussed in a detailed manner. For 
k = o and 2 this becomes the density of a hyperbolic cosine 

and logistic respectively. The characteristic function of this 

density is found for all values of k and compared with the 
characteristic functions of Normal, double exponential 
and Cauchy distributions. The moments are found in terms 

of the Euler and Bernoulli numbers and the first four 
moments are calculated and listed in a table for the Normal, 

logistic, hyperbolic secant and double exponential densities. 

Finally, linear combinations of independent random vari- 

ables with a double exponential distribution are considered 

and it is shown that by propérly choosing the combination 

we can arrive at the Normal, degenerate and logistic distri- 

butions. Further, a counter-example for the centrai limit 
theorem is constructed by the same method. If instead of 
the double exponential distribution we start with a hyper- 

bolic secant, it is pointed out that similar relationships can be 

established. 


(K. R. Parthasarathy) 
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WHITTLE, P. (Applied Maths. Lab., D.S.I.R., New Zealand) 


2.9 (2.1) 


2.1 (0.2) 


On the smoothing of probability density functions—In English 


J. R. Statist. Soc. B (1958) 29, 334-343 (6 refs.) 


Extending the method of a previous paper [P. Whittle, 7. R. 
Statist. Soc. B (1957) 19, 38-47] on the smoothing of frequency 
spectra, the author considers the estimation of a probability 
density function by means of a linear smoothing of an 
observed (sample) density. Thus, the estimate (at each 
value of the variate x) is found by evaluating a certain weight 
function (depending on x) at each of the sample values and 
averaging over the sample. The weight functions are to be 
determined by some optimality criterion. 

Some sort of model must be postulated if there is to be 
any basis for estimation. ‘The author proceeds to consider 
the curve being estimated as one of population of curves 
and for which population the correlation between adjacent 
_ ordinates approaches unity. This does not prevent any one 
curve from having many kinds of discontinuity and irregu-_ 
larity. ; 

The weight functions are then chosen to give unbiased 
estimates with minimum variance (relative to this popu- 
lation). ‘The problem reduces to that of solving a certain 
integral equation. It is shown that the resulting estimates 
are invariant with respect to monotone transformations of 
the variate scale. 

In the case where, in the population, the correlation 
between different ordinates depends only on the difference 
between the abscissae, the integral equation can be solved 
by Fourier transforms : at least for the case of a distribution 
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on the whole real line. A particular example is solved 
explicitly: for this example, it is shown that the effect 
of considering distributions over only a finite interval is 
negligible for large sample except at the ends of the interval. 

The asymptotic behaviour of the mean square deviation 
is examined and is shown in general to depend on the 
behaviour of the autocovariance function for small displace- 
ments. The details for the general case depend on the 
density near zero of the eigenvalues of the kernel in the 
integral equation for the weight function. 


(C. L. Mallows) 


DARWIN, J. H. (App. Maths. Lab., D.S.I.R., New Zealand) 
On corrections to the Chi-squared distribution—In English 


3-1 (0.4) 


J. R. Statist. Soc. B. (1958) 20, 387-392 (4 refs., 1 table) 


The null hypothesis distribution of the Chi-squared good- 
ness-of-fit test function (here called y) is considered when 
the null hypothesis is fully specified for a multinomial 
distribution. When the number of observations (N) is large 
the distribution of x? is approximately that of x? with k—1 
degrees of freedom. (K. Pearson, Phil. Mag., 1900). 

An asymptotic series of correction terms has been found 
by Hoel (Ann. Math. Siatist., 1938) where these terms form 
a sequence of continuous varying functions whose co- 
efficients are of increasing power in 1/N. Ii is also known 
that when N is small the discrete nature of the variation 
of x? renders such asymptotic expansion inadequate. It is 
this effect which is treated in the present paper. For large 
N as k increases, the jumps in the cumulative probability 
distribution function of x? are shown to be of order 1/N*/4, 
and an expression is obtained for their magnitude. 

The method used is to express the multinomial distribu- 
tion in terms of the Fourier inversion of the joint char- 
acteristic function of the multinomial variables. This is 
used to obtain the Laplace transformation of the probability 
density function of y? which can be expressed in terms of 
theta functions. The properties of these functions are used 
to express the Laplace transformation in such a way that it 
can be split into two parts ; one involving Hoel’s corrections 
and the other, the major term, retaining the discrete nature 
of its variation. ‘This latter is expressed as a multiple series 
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of functions which are first modified by neglecting a set of 
terms of order 1/eN and then Laplace-inverted in terms of 
Bessel functions of large argument, O(N?) which are replaced 
by their asymptotic expansions. The result is an expression 
which is shown to be of order 1/N*!4. 

The analysis is checked and illustrated by considering 
the binomial case, k=2, where the probability distribution 
function has only a simple series of discontinuities. The 
principal term in the expression for these is evaluated. A 
numerical illustration is given for the case N=10, p1=2/7 
where the cumulative probability distribution function of 
x? is considered at the jump points together with the 
approximation given by the cumulative probability distribu- 
tion function of x? with one degree of freedom together with 
two corrected expressions (using the author’s results). The 
two forms of correction arise since the correction term is not 
only a step function taking different values either side of a 
jump point but also takes an intermediate value at the 
jump point. This allows the cumulative probability dis- 
tribution function of x” (a step function) to be approximated 
by the cumulative probability distribution function of xX? 
with one degree of freedom evaluated either at the lower 
edge of the step with the correction added, or at the upper 
edge with the correction subtracted. In the case considered 
the author’s corrections give a considerable improvement. 


(D. E. Barton) 


3.8 (1.6) 


A contribution to the theory of Extreme Values—In French 


Publ. Inst. Statist., Paris (1958) 7, 37-121 (3 figs.) 


This paper is the first part of the thesis presented by the 
author for his Doctor of Science. 

The problem which sets the foundation for all this 
research is the study of stochastic behaviour of the ‘“‘ middle ”’ 
of a sample. That is to say, half the sum of the greatest 
and smallest values, i.e. the mid-range. The author supposes 
all the variables to be of a single dimension. The first 
chapter is devoted to the over-estimation and under-esti- 
mation of extreme values according to the three topologies 
used in this work: convergence in probability, almost 
certain convergence, almost completely certain convergence. 

In the second chapter the stability of extremes in the 
sense of Kolmogoroff is approached. The author in this 
way rediscovers some theorems of which certain are due to 
Gnedenko: the examples use the concentration function 
of Paul Lévy. 

The last two chapters are used to give the conditions 
of almost completely certain convergence and almost certain 
convergence (this latter case being the more difficult). The 
author has here concerned himself with new results. This 
paper will be followed by another which will be published 
in a subsequent number of this same journal. 


(D. Dugué) 
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HOOPER, J. W. (Netherlands School of Economics, Rotterdam) 
The sampling variance of correlation coefficients under assumptions of 


fixed and mixed variates—In English 
Biometrika (1958) 45, 471-477 (4 refs., 2 tables) 


The writer considers T independent observations of the 
M + A-variate (y, x) where y = (y,... yyy) andx = (Caperae an) 
from the same population. He contrasts the models (i) 
where (y, x) is a multinormal variate with zero mean, (ii) 
the regressive model where x is a parameter (‘“ fixed variate ’’) 
and where Y = &(y) = xm (where zw is an MA matrix) 
and y—Y is an M-variate multinormal. He proposes a 
“mixed ”? model to include (i) and (ii) as special cases (iii) 
where x = €+w where € is a parameter and w is a A-variate 
multinormal with zero mean ; the vector y —xz being multi- 
normally distributed independently of w. Under (iii) the 
asymptotic sampling variance (as T->00) of a sample canonical 
correlation coefficient (r) is given as 2T var (r) = (1—p?)? 
(2—p*p?), where p is the corresponding parent canonical 
correlation and p = &(€é x é’), the variables being in canonical 
form so that &(xx’) = 1, o<p<1: proof being left to an 
appendix. The asymptotic sampling covariance of any pair 
of sample canonical correlations is shown to be zero. 

Three special cases are considered: M=1, A=1; 
A>M=1; w =o (i.e. p = 1) corresponding respectively 
to bivariate (zero order) correlation, and models (i) and (ii). 
Tables are given of T var (r) for p? = o(o-1)1, p? = o(0°1)1, 
and of T var (R?) = 2R*(1 —R*)*(2 —R*p?) for R? = o(0'1)1 
and p? = o(0'1)1 (where R is the sample multiple correla- 
tion coefficient, M = 1, and R its population value). The 
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3-5 (3.9) 


writer notes that the asymptotic sampling variances are 
continuous in p, o<p<1, and that the choice of models (i) 
(p = 0) and (ii) (p = 1) or (iii) makes a considerable differ- 
ence in the asymptotic sampling variances of r and R (for 
large p, R) and in that of the zero order correlation (but 
not the regression) coefficient. The appendix proves the 
stated results by means of elementary but lengthy algebra 
heavily dependent on the underlying assumption of nor- 
mality : matrix notation is not used. 


(D. E. Barton) 


3-7 (3.9) 


On moments and cumulants of systems of statistics—In English 


Sankhya@ (1958) 20, 1-30 (16 refs.) 


The author presents a direct proof of Fisher’s sampling rules 
for k-statistics which is both reasonably rigorous and easier 
to follow than the earlier ones and shows how the results 
can be modified to obtain the sampling properties of a general 
set of statistics 21, 2%, 23, ... (homogeneous polynomial 
symmetric functions of the sample values, of degrees 1, 2, 
3, ... respectively). The proofs are mostly based on the 
Fisher-Kaplan approach, univariate and less general results 
being obtained in special cases. Simple explanatory examples 
have been introduced in various places. The author also 
defines a series of J-statistics whose expectations are popula- 
tion moments about the mean and whose sampling pro- 
perties can be investigated by using those of 2;. Rules are 
provided for obtaining (1) the expression of (ordinary or 
mixed) moments of the z, in terms of population moments, 
(2) the expression of (ordinary or mixed) moments of the ee 
in terms of the cumulants, and (3) the expression of (ordinary 
or mixed) cumulants of the z, in terms of cumulants. 

The corresponding rules for the 2-statistics appropriate 
to samples from a p-variate population are also given. 
Proofs of all the rules stated are given. 


(P. K. Kumaran) 
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MERRINGTON, M. & PEARSON, E. S. (University College, London) 
An approximation to the distribution of non-central t—In English 
Biometrika (1958) 45, 484-491 (5 refs., 2 figs., 1 table) 


3.2 (2.7) 


The probability density distribution and first four moments investigation demonstrates the practical utility of tables 
of non-central ¢ are stated. The relationship between the which enable a simple system of non-normal curves to be 
moment-ratios, 6, and B,, and the parameters f and § is used in approximating to the percentage points of other 
explored. A diagram is given of contours of constant f untabled distributions through the use of four moments. 
and 6 in the f,, B, plane. The beta points for the non- 
central ¢ are found to be entirely within the area correspond- 
ing to Pearson Type IV curves. As 5 increases the limiting (M. Merrington) 
form of the distribution is that of the reciprocal of x. 
Comparisons are made between the true upper and 
lower 5, 1, and 0-5 per cent. points and those found by 
interpolation in the ‘‘ Percentage points of Pearson curves,” 
Table 42 in Biometrika Tables for Statisticians. ‘True values 
were found either from N. L. Johnson and B. L. Welch’s 
paper: ‘‘ Application of the non-central ¢-distribution,” 
[Biometrika, 31, 362-389] or from ‘“‘ Tables of «he non- 
central t-distribution ” by G. J. Resnikoff and G. J. Lieber- 
mann. At the six percentage points, the Pearson Type IV 
curve provides a very good approximation to the non- 
central ¢-distribution over a wide range of values of f and 8. 
The difference between the true and approximate percentage 
points being never more than one in a hundred. 
Inaccuracies were found in the Resnikoff and Lieber- 
mann tables of percentage points for high values of f. The 
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SRIVASTAVA, A. B. L. (Statist. Lab., Indian Inst. of Tech., Kharagpur) 3-2 (3.1) 
Effect of non-normality on the power function of t test—In English 
Biometrika (1958) 45, 421-430 (12 refs., 3 tables) 


The author uses techniques which have been made familiar 
by Gayen and Ghurye. Starting with a parent population 
which may be described by four terms of Edgeworth’s 
series to obtain expressions in terms of the Hh functions 
by means of which the power function can be calculated. 
Some calculations are given showing the effect of non- 
normality on the power function for samples of sizes 5, 10, 
and 20 and for the standardised difference of means 0, I, 2, 
3, 4, 5. The inverse problem of the determination of the 
critical region corresponding to a predetermined level of 
significance is also considered using the Cornish-Fisher 
expansion. The general conclusion is that unless the parent 
population represented by the Edgeworth series is markedly 
non-normal the effect of non-morality is not marked. 
When kurtosis in the parent population is not very different 
from that in a normal population the skewness effect becomes 
more pronounced. 

There is an appendix by E. S. Pearson in which the 
correspondence between the theoretical results of Mr. 
Srivastava and the early (1929) empirical sampling results 
of Pearson and Adyanthaya is pointed out. 


(F. N. David) 
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BAHADUR, R. R. (Indian Statistical Institute, Calcutta) 


Examples of inconsistency of maximum likelihood estimates—In English 


Sankhya (1958) 20, 207-210 (4 refs.) 


The maximum likelihood (M.L.) method provides an 
estimate of the entire population distribution even if one is 
interested only in a particular parameter. It is, therefore, 
of interest to enquire whether the maximum likelihood 
estimates of the population distribution converge, in some 
given intrinsic sense, to the actual distribution as the sample 
size tends to infinity. After remarking that Wald’s famous 
proof of the consistency of the M.L. estimates can be 
formulated so as to yield regularity conditions which ensure 
convergence in this sense, the author states that the need for 
regularity conditions is not clear from previously published 
examples of the failure of the M.L. method. This is because 
these examples are based either on convergence definitions 
which are far from intrinsic or on the discontinuous nature 
of the parameter to be estimated, considered as a functional 
on the set of alternative distributions. 

The author then gives two examples where maximum 
likelihood estimates of the population distribution function 
fail to converge to the actual distribution function: con- 
vergence being understood in the sense of P. Lévy. In the 
first example the requirement that the estimate maximise 
the likelihood is so inconsequential that some M.L. estimates 
do not converge at all. The set of alternative distributions 
in this example is the set of all probability distributions on 
[o, 1] which are absolutely continuous, having continuous 
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BASU, D. (Indian Statistical Institute, Calcutta) 
On sampling with and without replacement—In English 
SankhyG (1958) 20, 287-294 (3 refs.) 


The author considers a population of N identifiable units 
and discusses estimation procedures for the mean of some 
observable characteristic in this population under various 
types of sampling schemes—sampling with or without 
replacement and with equal or unequal probabilities. 

In sampling with replacement and with equal probability 
for a fixed sample size of m the author shows that the mean 
of the v distinct units in the sample (ignoring repetition of 
units that happen to occur more than once in the sample) 
is an unbiased estimate for the population mean and has 
uniformly a smaller variance than the sample mean. ‘The 
proof is similar to that of Rao-Blackwell theorem and 
utilises the fact that the observations corresponding to the 
distinct sample units constitute a sufficient statistic in this 
problem which, incidentally, is not a complete one. The 
same result is proved in a different way by Des Raj and 
Khamis [Ann. Math. Statist, 29, 550]. 

In sampling with replacement and with unequal pro- 
bability, the author shows that the observations corresponding 
to the v distinct units along with the associated probabilities 
of selection are sufficient. Hence it is possible to improve 
upon the sample mean by taking its conditional expectations 
given the values for the sufficient statistic. The actual 
computation of this conditional expectation is, however, 
very complicated, and even more troublesome is the com- 
putation of an estimate for its variance. 


4-1 (4.3) 


densities f(x) such that oX f(x)<2. In the second ex- 
ample the author constructs a countable family of admissible 
distributions on the set of integers, such that although M.L. 
estimates exist, they always converge to a fixed one from 
the admissible set, no matter what the actual distribution 
may be. 


(V. S. Varadarajan) 


4-9 (8.2) 


The author also considers sampling with replacement 
and with equal probability for a fixed number » of distinct 
units and shows that the mean of all the n sample obser- 
vations, though an unbiased estimate of the population mean, 
has uniformly a larger variance than the mean of the v distinct 
units. 

The author similarly considers sampling without replace- 
ments and with unequal probabilities and shows that the 
set of ordered sample observations along with their corres- 
ponding probabilities of selection are sufficient in this case. 
Hence any estimator which makes use of the order in which 
samples were drawn can be uniformly improved upon by 
taking its conditional expectation for given values of the 
sufficient statistic. 

The paper also gives the sampling distribution of », 
the number of distinct units in a sample of fixed size drawn 
with replacement and with equal probability and gives only 
the first two moments of this distribution for the unequal 
probability case. 


(S. K. Mitra) 
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BULMER, M. G. (Department of Social and Preventive Medicine, Manchester) 


4.4 (7-1) 


Confidence intervals for distance in the analysis of variance—In English 
Biometrika (1958) 45, 360-369 (9 refs., 6 tables, 1 fig.) 


A “ parametric’? model of the analysis of variance is con- 
sidered, E(Y) = XB where Y is a nX1 column vector of 
observations with a variance o”, X is an nx k matrix of known 
constant coefficients and 8 is a kx 1 column vector of para- 
meters. ‘The null hypothesis specifies r(<k) parameters 
(or linear independent functions of these), denoted by Vi 
and rx 1 column vector. 

Two types of confidence interval for the ‘‘ distance ”’, 
6, a function of (Y—Y,), of the true hypothesis, Y, from 
the null hypothesis, are constructed : (i) a “ simultaneous ”’ 
confidence interval, ensuring at least a preassigned confidence 
coefficient, for the “ distances ”’ from several null hypotheses 
—which can also be used in conjunction with Scheffé’s 
(listed in references) method of assigning simultaneous 
confidence intervals on individual parameters; (ii) an 
approximate confidence interval narrower than (i), for this 
“distance ”’ from a single null hypothesis, with a confidence 
coefficient as near as possible to some preassigned value. 
The methods discussed can also be used to find, in chi- 
Square test for goodness of fit, the confidence interval for the 
distance between the true and null probabilities of an 
observation falling in a particular class. 

The confidence interval (i), expressed as a simple function 
of the “testing” and “error”? sums of squares, is illus- 
trated by an example on analysis of variance. The idea is 
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extended to find the confidence intervals separately, for 
several components into which Y can possibly be split up 
—an example of chi-square test for goodness of fit has eluci- 
dated such an application. 

The confidence interval in (ii) is obtained by approxi- 
mating this distribution of the minimum value of the 
“testing? sum of squares, under a particular alternative 
Y = Yj, by a chi-square distribution with correct mean 
and variance. The actual computation of this approximate 
confidence interval, as illustrated by an example, involves 
mainly two quantities: (a) X = 8/0”, which can approxi- 
mately be obtained by an interative process from the percen- 
tage points of F-distribution with fractional d.f.’s one of which 
is a function of A, (6) the evaluation of the percentage points 
of an F-distribution for fractional d.f.’s f, and f,, a special 
table having been provided for 1 and 5 per cent. when 
fiS2°2 and fg =~. The adequacy of the approximation 
is investigated in the last section where it is found that 
the true confidence interval is slightly narrower than is 
asserted by the approximation. 


(B. K. Ghosh) 


4-4 (3-3) 


Statistical inference on the parameters of non-normal populations—In Spanish 


Trab. Estadist. (1958) 9, 117-140 (7 refs.) 


A discussion about statistical inference in cases where the 
parent population is continuous with finite mean and finite 
standard deviation, but non-normal. Having in mind the 
possible application of results to practical situations, the 
author sets as a postulate that finite populations, in the case of 
simple random sampling, are random samples of infinite popu- 
lations : ‘‘ Not a random sample in the restricted sense given 
in the standard texts of sampling theory and methods, but 
rather in the more general and intuitive sense which serves 
as a basis of all statistical considerations.” 

Taking a specially defined sample of means from con- 
tinuous distributions, conditioned with respect to a previous 
random sample, the author discusses the possibility of 
obtaining a sample of means of values, when these values 
are directly obtainable from the population. 

The case is then discussed of confidence intervals when 
its length, or the confidence coefficient, are random vari- 
ables ; the size of the sample being previously fixed. Substi- 
tuting an approximate confidence coefficient for the exact 
one, confidence limits are set for the mean of the sampling 
distribution of sample means: considering both the cases 
of infinite and finite parent populations—using, in the latter 
case, the author’s postulate mentioned above. Afterwards 
the author’s attention is directed towards stratified sampling 
and he proceeds to extend his first postulate considering each 


finite stratum or sub-population as a random sample of 
an infinite population, in a way similar to the previous one 
for a non-stratified situation. 

Next is discussed the case of confidence intervals when the 
number of the elements in the sample is a random variable, 
the confidence coefficient and the length of the intervals 
being fixed in advance, or at least an upper bound for the 
confidence interval. "The author derives his proof from a 
lemma which is the well-known property that a variate 
t may be generated by dividing a normal variate with zero 
mean and unit standard deviation into the square root 
of an independent chi-square variate divided by its number 
of degrees of freedom. The two cases of infinite and finite 
populations are also considered and confidence intervals 
are set for the mean of the sampling distributions of sample 
means. Rules are also given for the computation of the 
sample size, depending on the number of samples: Special 
consideration being given in this case to stratified popu- 
lations. 


(F. Azorin P.) 


DUNN, 0. J. (Iowa State College ; now at University of California, Los Angeles) 4.2 (4.8) 
Estimation of the means of dependent variables—In English ai 
Ann. Math. Statist. (1958) 29, 1095-1111 (7 refs., 2 tables) 


Several methods are given for deriving sets of simultaneous 
confidence intervals for the means of a multivariate normal 
distribution. 

Exact simultaneous confidence intervals are obtained 
(Section 3) for the case of known variances and for the case 
of equal, but unknown variances. Knowledge of the 
covariances is unnecessary. 

Several sets of intervals having bounded confidence level 
are obtained (Sections 4-7) for the two aforementioned 
cases. In addition, the method which employs Hotelling’s 
T-statistic is used (Section 4.2) to derive intervals for the 
general case where the variances are unknown and arbitrary ; 
here the sample size is assumed greater than the number of 
variables under study. 

Finally, the various methods presented are compared in 
terms of the expected values of the lengths of resulting 
intervals. This is done numerically, for confidence co- 
efficient 0:95 and several sample sizes. For general use, 
the winner in this respect is the method of Section 6 based 
on a Bonferroni inequality. 

The problem of dependent variables in which means of 
one or more variables depend on observations (not just 
means) of one or more other variables, though illustrated in 
examples in Section 2, is not treated. 


(Robert Hooke) 
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EL-BADRY, M. A. (Institute of Statistics, Cairo Univ.) 4.1 (4.3) 
The optimum estimate of moments—In English 
Egyptian Statist. F. (1957) 1, 33-38 (2 refs.) 


The problem of estimating p,, the kth central moment 
(k>1) of a distribution, of which a set of statistically inde- 
pendent observations is known, is considered. 

The author defines the optimum polynomial estimate as a 
polynomial in the observation such that its expectation is 
equal to p,, and which has the minimum variance in the class 
of unbiased estimates. 

He shows that such an estimate must be 

I. a symmetric polynomial in the observations 

2. homogeneous of degree k 

3. and that only one of the observations in each of the 

terms of the polynomial can have a power greater than 
one. 

For this he makes the assumption that the moment of 
degree 2k exists and gives an explicit expression for the 
estimate. 

He applies his result to give simple expressions for the 
estimate of the variance, and the estimates of the third and 
fourth central moments. 


(J. B. Simaika) 
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FERGUSON, T. (University of California, Los Angeles) 


4.2 (2.4) 


A method of generating best asymptotically normal estimates with application 
to the estimation of bacterial densities—In English 


Ann. Math. Statist. (1958) 29, 1046-1062 (18 refs.) 


The minimum x? method of estimation introduced originally 
by Karl Pearson and designed for use with the multinomial 
distribution has been extended and modified by Neyman 
and his co-workers in the past several years. These modi- 
fications have been introduced to simplify the calculation 
procedures while retaining the properties that these estimates 
have in common with maximum likelihood estimates, viz. 
asymptotic normality and minimum asymptotic variance 
(under suitable regularity conditions). Estimates with these 
properties are called best asymptotically normal estimates. 
A review of four modifications of the minimum x? method 
is given both for the multinomial distribution and for a more 
general situation. References are given to papers giving 
examples of the application of these methods. 

In all of these methods the estimation procedure involves 
the minimisation of a quadratic form in the observations. 
The author shows that it is also possible to start directly 
with a linear form and obtain best asymptotically normal 
estimates and, further, this procedure may be coupled with 
any of the usual modifications. 

The procedure developed is applied to the problem of 
estimating bacterial density by the dilution method. In 
this method, from a large volume of the liquid containing 
A bacteria per unit volume (A unknown), samples of size 
%; unit volumes are withdrawn and tested for presence or 
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(University of Aberdeen & A.R.C. Unit of Statistics) 


absence of bacteria. The probability that no bacterium appears 
in a single sample of volume «; is p; = e % The estimate 
of A is based on the number of sterile samples and since 
p; is an exponential function of A, the maximum likelihood 
estimate is the root of a complicated non-linear equation. 
References are given to papers where methods to solve the 
equation are discussed. ‘The author’s modification of mini- 
mum x” estimation still leads to an exponential equation 
but one claimed to be easier to solve than that yielding the 
maximum likelihood estimate. Numerical studies of the 
small sample properties of the new estimate are promised 
for a later paper. 


(D. G. Chapman) 


4-3 (10.7) 


The efficiencies of alternative estimators for an asymptotic regression equation—In English 


Biometrika (1958) 45, 370-388 (9 refs., 4 tables) 


This paper investigates the properties of various methods 
of estimating the parameter in the regression 
E(y)=a—B p?. 

Two models are considered for the generation of the devia- 
tions: (@) the y are independent and normally distributed 
about their expectations with constant variance, and (db) the 
y arise from a continuous autoregressive process in which 
the variance increases with x and successive values of y are 
correlated. 

The method of statistical differentials is used to. obtain 
general expressions for the bias and variance of a ratio. In 
the further investigation, detailed results are obtained only 
for the case (case 4, say) of four observations equally spaced 
in x. For this case, explicit expressions (in terms of p) are 
obtained for the following estimators, for both the models 
(a) and (6) above ; these are evaluated in four tables. 

(i) Patterson’s linear estimators. ‘These are the ratios 
of weighted averages of the y (omitting the first) to 
the same weighted average of the y (omitting the 
last). In case 4, there is one degree of freedom in 
the weights, and three values (1, 1-25, 2) of an 
arbitrary parameter are suggested as being reasonable 
for different ranges of p. This estimation procedure 
has high asymptotic efficiency (as the variance of the 
y approaches zero). 


(ii) The estimate obtained by performing the standard 
regression of y,,, on y,. For model (6), this gives 
the maximum likelihood estimator. In case 4, 
model (a), it has high efficiency, and the bias is of 
the same order as for (i), though of opposite sign. 
Two estimates proposed by Hartley, using more 
complicated autoregressions of the y. In case 4, in 
comparison with (i) and (ii), these have similar 
variances but rather larger biases. [Biases not 
evaluated for model (6).] 

A “general quadratic estimator” (ratio of two 
quadratic functions of the y) of which in case 4 a 
particular example is chosen for study. For model 
(a), this has high efficiency (bias not evaluated) but 
has no special merits. 

This study is continued by Patterson, Biometrika (1958) 


45, 389-400 (q.v.). 


(iii) 


(iv) 


(C. L. Mallows) 


yg en Pe 


» Bert nah epee vane 
ee Bags ‘te 


te x haga rt Nis aed ape tye wat 
peranierii( "i! sph eB, Pa ‘hae is 
RRR s ive) eat as 
F bide ed Be i af 


af Ati Siihe the oa) 


tavoeta vst Sa Ae i Vy st 3 | Bye 4 


Aa 5 


ase Sere yc o: 
Pe aati i Lira esi 
a he ut Bist 


asa a Na i il ali 


know a priori (from sources other than the sample the 


a simple approximate v: iar 


certain sums are | 


GUTTMAN, I. Ve eae 
MAN, 1 (ni solution of a 


Biometrika (1958) 45, 365-367 ‘refs.) 4 ie ee 
ese Bs hat eee ; oes i, eae. + 
Suppose ¢ is a complete sry auicisgt statistic, ssumi eee ee cree ne oa 


only non-negative integer values, having a probability a 
distribution of exponential type. depending on sarameter RMT Scie ts ah, 


@ which takes values in an interval basi eed «ah ane er ay Sar pee 
Then there exists an (essentially unique) unbiased estimate cet Ae eee ae em a to, 
with uniformly minimum variance of any chosen function 4 ee Sate ee oa N wos e era eee 
of @ if, and only if, a certain simple function is analytic at bees ate ds bl eee Ne eee 
one: and has a power series expansion whose rth coefficient 1 Rabe Se eara Rad 9 ce tek 
is zero whenever the probability that ¢ equals r is zero. The Re Se a ae eT ic Seda at 
estimate is obtained very simply in terms of these coefficients. | RR eid Rent Sy. Heli ret Le Ren te 
This result facilitates the construction of unbiased ae bb le fa ep DR See 
uniformly minimum variance estimates without recourse to Laeger ae Lee TCT 
the usual method involving the evaluation of the conditional oh a ey en Le sce ene 
expectation of an arbitrary unbiased estimate. It gives a weet nt “ Ptah e a eae 
result in cases where the method given in the 1946 paper by ah vies Cobh tm rsa et gant Se Rian has 
Girshick, Mosteller and Savage [Ann. Math. Statist. 17, 13] | bce rues 
fails. ‘The method is applied to the cases of negative binomial PRN ie Senay cae 
and Poisson distributions. ns ee Pree 
(C. L. Mallows) penta ae Cee aaa 
5 : ~ 7 : ; rs Lae 
‘ i Cito Rt oy, 
49 | o | ek (Lon Cope int er ae 


HAJEK, J. (Math. Inst., ois ctonald pre of soioees Prague) _ Has 
On the theory of ratio estimates—In English ast eee 
Aplik. Mat. (1958) 3, 384-398 (9 refs.) PEs, Papper bs Sri. 
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by simple random sampling, we wish to estimate the aggre- gested by the linear approximation, 
gate value of a variable, say y, under the supposition that we the ratio estimate is c 
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LEJEUNE, J. 


(Université de Paris) 


4.3 (2.8) 


On an a priori solution of Haldane’s a posteriori method—In French 


Biometrics (1958) 14, 513-520 (4 refs., 3 tables) 


The binomial distribution with the zero observations trunc- 
ated frequently occurs in human genetics. When data on 
families in which a defective (taré) trait is transmitted are 
available, frequently the families comprising the data are 
such that at least one of the offspring shows that trait. This 
is especially the case when the trait is governed by a recessive 
gene. Under these conditions estimation and tests of the 
proportion, p, of the offspring carrying the trait tend to be 
numerically tedious. 

The author rejects the procedure of Bernstein and 
Hogben for testing p = 1/4, 1/2 and is concerned to simplify 
the application of the method of maximum likelihood esti- 
mation of p as developed by Haldane. The estimation 
equation is equivalent to that obtained by the method of 
moments and is a linear function of p/{1—(1—p)"}; no 
explicit solution is available. 

This paper presents tables for facilitating the compu- 
tation of the maximum likelihood estimate and its estimated 
variance. ‘The first table tabulates the above function to 
four decimal places fo r m =2(1)15 and p = 0:60(0°05)o'05 
and 0-0275, 0:0225. 

The use of the tables is as follows: an arbitrary value, 
usually 0°25, 0°50, is chosen for p and the expected number 
of abnormal offspring is calculated according to the formula. 
This is repeated until the two neighbouring values p,, po, 
which give expected numbers straddling the observed 
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number, are discovered. At this stage linear interpolation 
1s accurate to within 1 per cent. 

Haldane has also given an expression for the inverse 
of the variance of the estimator. The author rearranges this 
into a linear function of a function of p, n which is tabulated 
for the above values of p, n. The inverse variance is calcu- 
lated for p,, p, and that for the interpolated value of p 
obtained by linear interpolation on these quantities. 
Although this underestimates the variance slightly, it is 
within 1 per cent. and so sufficiently accurate for practical 
purposes. 

The use of both tables is illustrated by application to 
some data on ataxia. The terms a priori, a posteriori are used 
by the author to refer to probability and relative likelihood 
calculations respectively. ‘Thus the title seems to arise 
from the previously noted equivalence of maximum likeli- 
hood and moment estimators. 


(C. W. Clunies-Ross) 


4-3 (2.5) 


Estimation of parameters of mixed exponentially distributed failure-time 
distribution from censored life test data—In English 
Biometrika (1958) 45, 504-520 (17 refs., 6 tables, 3 figs.) 


This paper is concerned with the estimation of population 
parameters for sub-populations where the measured char- 
acteristic of both the sub-populations and the whole popula- 
tion is distributed exponentially. The particular field of 
application is in connection with life test data derived from 
sample observations censored with respect to the time at 
which the test is terminated. After a brief introduction, a 
statement of the population model and the method of 
sampling, the authors deal with the estimation of the para- 
meters according to whether or not the relative magnitude 
of the sub-population parameter is known. ‘The method of 
maximum likelihood is used and estimating equations are 
given for the case of two sub-populations, i.e., the parameter 
for each together with the proportion that the size of one 
bears to the whole. These estimating equations involve an 
additional factor for which an estimating equation is also 
given together with an iterative process for estimation based 
upon a modification of the MLE obtained in 1956 by 
Deemer & Votaw [Ann. Math. Statist. 26, 498-504] for 
samples drawn from a singly truncated exponential distribu- 
tion. In the case where the relative magnitude of the sub- 
population parameters is known and the given procedure 
yields estimates which reverse the relationship, an adjusted 
estimation procedure is given: this condition is termed a 
“ crossover ”’ of the estimated sub-population parameters. 

The small sample properties of the estimates were 
obtained by empirical methods for a number of parameter 


points, i.e., combinations of sample size and value of the 
sub-populations’ parameters. ‘The results are shown in a 
full table and it is concluded that the estimation is poorest, 
both with respect to bias and variance, for those parameter 
points at which a large number of crossovers occurred: a 
substantial improvement is achieved by using the adjusted 
procedure. The large sample variances were obtained by 
inverting the symmetric information matrix and a short 
table is given comparing the asymptotic variances with the 
estimated variances for five parameter points. The effect 
of the level at which the sample is censored is considered 
in relation to the efficiency of the estimates—here in con- 
nection with the problem of whether the cost of an increase 
in the testing time will be offset by enough additional 
information. 

After showing that the results for the mixed exponentially 
distributed sub-populations are but a logical extension of 
those for a single exponential distribution, the authors 
generalise their model to a situation for which the exponential 
distribution is a particular case. The model chosen is the 
well-known family of distributions associated with the name 
of Weibull. ‘The authors state that it should not be too 
difficult to set up an iterative method for the solution of 
the ML equations provided the number of sub-populations 
is small: reference is made to the use of standard numerical 
methods for the solution of equations by iteration with or 
without the use of a computer. 


(Wm. R. Buckland) 
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PATTERSON, H. D. (Rothamsted Experimental Station, Harpenden, England) 


4-3 (10.7) 


The use of autoregression in fitting an exponential curve—In English 


Biometrika (1958) 45, 389-400 (5 refs., 4 tables) 


This paper complements the investigations of Finney, 
Biometrika (1958) 45, 370-388 (q.v.) into the properties of 
certain estimators of the “base”? parameter p in an ex- 
ponential regression. ‘The residuals are assumed to be 
independent, normal with zero means and constant variances. 


(i) The simplest estimator considered [denoted by (ii) 


in the summary of Finney’s paper] is that obtained 
by straightforward analysis of the regression of Void 


on yz. 
(ii) More generally, the regression of y, 41 on a weighted 
average ky, +ly,,, is considered. 
(iit) Finally, a “‘ general quadratic estimate’ (ratio of 


two quadratic functions of the y) is considered. It is 
not obvious whether this is more or less general than 
Finney’s estimator (iv). 

Using matrix notation, general expressions for the lead- 
ing terms in the bias and variance of (iii) are obtained. The 
variance is independent of kand/. It is shown how to choose 
the estimator so as to minimise the variance for any chosen 
value of p, and suggestions are made for computing the 
estimator. 

Taking p equal to zero [which gives (ii) above], the 


Finney’s paper]. It is found that the ‘“ regression ” 
estimator is substantially less efficient than the “linear ” 
estimator in the case of 7 observations, over the range of 
the most useful values of p. The bias of the “ regression” 
estimator is a function of k and J, is tabulated for case (i) 
above, and is compared with that of the “‘ linear ”’ estimator. 
The bias is large, and becomes more serious with increasing 
numbers of observations. It can be controlled by suitably 
choosing k and J, and is exhibited for some specimen values. 
The choice is difficult unless the number of observations is 
small. 

The conclusion is that in general the regression method 
is unsuitable. Preliminary results assuming the residuals 
to be autocorrelated [model (6) in the abstract of Finney’s 
paper] indicate that Steven’s method (Biometrics, 1951, 1, 
247-267) can safely be recommended for most biological 
applications. 


(C. L. Mallows) 


efficiencies of these estimators are tabulated for the cases 
of 4, 5, 6, and 7 observations, and are compared with those f 
of Patterson’s “‘linear”’ estimators [(i) in the abstract of i 
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RADHAKRISHNA RAO, C. (Indian Siatistical Institute, Calcutta) 4.1 (4.3) 


Maximum likelihood estimation for the multinomial distribution with 


infinite number of cells—-In English 
Sankhya (1958) 20, 211-218 (8 refs.) 


In this paper the author extends his earlier results [Sankhya 
18, 139] concerning maximum likelihood (M.L.) estimation 
to the case of a multinomial distribution with an infinite 
number of cells. 

The author shows, under the sole restriction that 27; log 
i> —©0 where the x; represent the population frequencies 
in the different cells, that the maximum likelihood estimate 
of the population distribution converges to the actual distri- 
bution in the strongest of intrinsic senses. This result 
can also be proved by using the arguments of Wald [Ann. 
Math. Statist., 20, 595] as pointed out by Kiefer and Wolfo- 
witz [Ann. Math. Statist., 27, 887] although the author’s 
proof is simpler and more direct. In the parametric situa- 
tion, the author proves, under the further restriction that 
the parameter to be estimated is a continuous functional of 
the distribution, that the approximate M.L. estimates 
corresponding to any fixed degree of approximation c 
(o<cX1), are uniformly consistent. 

In the next part of the paper, the author discusses 
maximum likelihood equation (M.L.E.) estimates, i.e. 
those which provide the maximum likelihood when the 
parameter is restricted to the roots of the likelihood equation. 
Under the additional assumption that there exists some 
neighbourhood of the true value of the parameter in which all 
the cell probabilities can be differentiated, the author demon- 


strates that, with probability one, for all sufficiently large n 
(n being the size of the sample), maximum likelihood 
estimates exist, that the class of M.L. estimates coincides 
with the class of maximum likelihood equation estimates, 
and that the M.L.E. estimates are uniformly consistent. 

The author remarks that the regularity conditions imposed 
by him in the derivation of the above results seem to be consid- 
erably weaker than those assumed by earlier authors, and that 
the methods developed in the present paper apply without 
modification to the situation when the number of parameters 
to be estimated is greater than one. 


(V. S. Varadarajan) 
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SARHAN, A. E. (High Institute of Public Health, Alexandria) 
Estimation of central tendency for asymmetric distributions—In English 


Egyptian Statist. ¥. (1957) 1, 7-14 (4 refs., 4 tables) 


The author considers three measures of central tendency 
(location): the mean, the median and the mode in asym- 
metric distributions. He estimates these measures by the 
unbiased linear minimum variance method and compares 
the variances of the different estimates for three particular 
asymmetric distributions. 

He shows that for 1. the two-parameter single exponential 
distribution and for 2. the one-parameter single exponential 
distribution, the variance of the sample median is approxi- 
mately one-half the variance of the sampie mean. While 
for 3. a skewed distribution with finite range (the mode 
dividing the range in the ratio 2:1), he shows that the 
variance of the sample mode is slightly smaller than the 
variance of the sample median and that this last variance is 
slightly smaller than the variance of the sample mean. 


(J. B. Simaika) 
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SHENTON, L. R. (College of Science and Technology, Manchester) 
Moment estimators and maximum likelihood—ZIn English 


Biometrika (1958) 45, 411-420 (16 refs., 2 figs.) 


A probability distribution function p(x ; 0,...4,) is dependent 
on A unknown parameters, is continuous in x, possesses all 
moments and is such that the range of x is independent of 
0 = (4,, ..., %,). The polynomials {g, (x)} dependent also on 
6, are defined by their orthogonality with respect to expecta- 
tions over p(x ; @). A series of these polynomials (to order r) 
is fitted (by moments) to the differential with respect to 
6; of log p(x; 6). The fitting polynomial is denoted 
A; («) and is sample mean value, denoted ./;'”, is then a 
linear sum of the first y sample moments (whose coefficients 
are functions of @) which has zero expectation. ‘The equa- 
{IOTIGN. Lj) On f= eT aces , h are then solved for @ to get 
estimates 6,, ...6,, which are functions of the first r sample 
moments. By what is in effect an application of Cramer’s 
theorem on functions of sample moments, the leading terms 
in the expansion of cov. (6;, 6;) as power series in 1/N are 
obtained and these are shown formally to approach the 
corresponding maximum likelihood results of Fisher as 
y->«. The approach to the limit is shown to be monotonic 
for h = 1 and the corresponding result for var. (0,)1 = 1, 
mt anes h is stated to hold generally (a method is given which is 
said to provide the outlines of a proof ). Three distributions 
are taken to exemplify the method. The first is the truncated 
Poisson where the method yields David & Johnson’s sufficient 
moment or maximum likelihood statistic. ‘The second is an 
artificial example (the product of e~* by a function linear 


4-1 (4.3) 


in x and the unknown parameter) where the estimator turns 
out to be the root of an rth order polynomial. ‘The third 
is the negative binomial of mean A and variance A(1+A/a). 
The estimator of X turns out to be the usual maximum like- 
lihood estimator (the sample mean) and & the root of a 
(2r—3)rd order polynomial. This, the author states, he is 
unable to show has a unique root but he shows that for 7 as 
small as 4 his estimator has high efficiency for « relatively 
small. The paper closes with a brief discussion of the 
analytic conditions on p(x; @) necessary to ensure that the 
efficiency of the estimators tends to unity as r—> » but no 
useful conclusion is reached. There is a final section 
referring to earlier papers of the author’s where analogous 
problems arise. ‘The author does not investigate the small 
sample behaviour of these estimators in this present paper. 


(D. E. Barton) 
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TANNER, J. G. (Road Research Lab., D.S.LR., England) 


A problem in the combination of accident frequencies—In English 


4.3 (5.2) 


Biometrika (1958) 45, 331-342 (2 refs., 3 tables, 2 figs.) 


The paper is concerned with the accident frequencies in a 
period before and after similar changes in road conditions 
have been made at a number of sites. Comparisons are made 
using a set of control areas for which it is assumed that the 
ratio of accidents in the ‘ before’ and ‘ after’ periods is free of 
error. This use of control sites permits that the period of 
observation on the changed sites may differ. The total 
number of accidents occurring at a given changed site 
is then assumed to divide binomially into the number of 
accidents in the before and after change periods, the prob- 
ability parameter of the binomial distribution being 
expressed as a function of x; which is a measure of the change 
at the 7th site and takes the value unity when no change 
occurs. 

In the first instance it is assumed that there is no varia- 
tion between the values x; at N sites and maximum likeli- 
hood applied to estimate the common value « and to test the 
hypothesis x = 1, i.e. that no change in accident frequency 
has been observed. 

Later an expression which measures the variation 
between the x; is constructed and which may be used to 
test the hypothesis that the «; have a common value, not 
necessarily unity. 

So that inference may be made with regard to all sites 
which undergo the similar change in road condition rather 
than to the N sites which are under direct investigation, 
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it is finally supposed that a parent population exists and 
that the values of x; form a sample from this population. 
Estimates of the first and second moments of the parent 
population are proposed together with a method for testing 
the hypothesis that the population first moment is unity. 

In an appendix the first four moments of the logarithm 
of the estimated population first moment are given to first 
order when all the x; are equal. 


(J. G. Saw) 


4.4 (3.6) 


Confidence limits for the mean of a normal population with 


known coefficient of variation—IJn English 


Aust. J. Appl. Sci. (1958) 9, 321-325 (1 ref., 2 figs.) 


Let « be the arithmetic mean of a random sample of size n 
from a normal population whose mean p and variance o* 
are unknown. If the coefficient of variation k = o/w of 
the distribution is known, the confidence limits for pw are 
determined by */(1+A) and x/(1—A), where A = Tkv/n, 
and the x confidence interval for the standardised normal 
variate is (—T, T). These limits depend only on k, ¥, and n. 


Let 


(i) 1p be the width 2Tky/1/n of the shortest confidence 
interval for », when the standard deviation o is known ; 

(ii) 1, be the corresponding width 2ts/./n of the con- 
fidence interval when o is unknown, ¢ being the value 
corresponding to 7 obtained from Students’ distribution 
for the sample of size n, and s the standard deviation of the 
sample ; 

(iii) 1 be the width 2 Ax /(1 — A®) of the confidence interval 
when k is known, using the above confidence limits. ‘Then 
by calculating the expected values Ao, A, and A of 1p, 1, and 1 
respectively, the efficiency of confidence limits based on 1 
is compared with those based on |, and l,. 

It is found that A, is always the shortest of the three. 
For small values of k, X and A, differ only slightly, while 
the difference between A, and A; is much larger. Only for k 
larger than 0:5 (approximately) does 4, turn out to be smaller 


than A. For large values of n, both A and A, approach rapidly. 
By plotting 100(A—Ag) /Ap and 1oo(A; —Ag)/Ag against different 
values of k, the expected values 4, A,, and A») are compared 
graphically for samples of size 4, 6, 8, 10 and 30, and different 
values of k, for 95 and 99 per cent. confidence intervals. 

An expression for the approximate expected value of s, 
the sample standard deviation is given in an Appendix, 
with an upper bound for the error involved. 


(B. R. Bhat) 
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ADKE, S. R. & DHARMADHIKARI, S. W. (University of Poona) 


5-7 (2.7) 


Gain due to sequential sampling from a Gamma population—In English 
Bull. Calcutta Statist. Ass. (1958) 8, 91-94 (4 refs., 1 table) 


It is well known that the use of sequential probability ratio 
test (SPR test) usually results in considerable saving in the 
average number of observations required to reach a decision, 
compared to that required by the classical Neyman-Pearson 
test of the same power. Wald has evaluated this gain, when 
observations are drawn from a Normal population with a 
known variance, to test the simple hypothesis JH, 0 
(mean @ = 69) against the simple alternative H, (@ = 0,). 
In this note the authors obtain the corresponding results 
for the Gamma population with mean Jo and variance Ic? 
for testing the hypothesis Hy (o = op) against H, (¢ = a) 
assuming that / is known. It is proved that the percentage 
gain in the number of observations is independent of 1. 
Hence, if 7 is taken as unity, the problem can be treated as 
one pertaining to the exponential population with mean oc. 
This simplified problem is treated by Epstein and Sobel 
[F. Amer. Statist. Assoc., 48, 486, and Ann. Math. Statist., 
26, 82]. The authors have reproduced the numerical results 
obtained by them for some values of K = o,/o,, « and B 
(the first and second kind of errors). The percentage gain, 
in these particular cases, is seen to range from 32 per cent. 
to 78 per cent. 


(K. Madhava Sarma) 
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5.2 (8.9) 


A significance test for fraction-defective by compressed limits—In Fapanese 
Statist. Qual. Cont. (1958) 19, 54-62 (8 refs., 1 table, 11 figs.) | 


When the probability that the variable under investigation 
exceeds the specification limits (in which case the item with 
this value is assumed to be defective) is very small, a large 
number of observations are needed to detect the change of 
the population mean of the variable by attribute; for 
example, by ‘‘ go-not-go ”’ device. It is well known that, in 
such a situation, compressed limits can be effectively used 
instead of the original specification limits in order to decrease 
the number of observations used in the control chart (cf. 
B.S.S. 1313 (1947): Fraction-defective charts for Quality 
Control). 

The author has applied this device to the significance 
test of the null hypothesis that the population mean has a 


given value against the one-sided alternative based on the 


-number of defectives and, given the method of obtaining 
the compressed limit which gives the minimum number of 
observations, to detect the given difference of two values of 
the population mean—under the null and alternative 
hypotheses for the given values of the probability of the 
first and second kinds of errors. The distribution of the 
basic variable is assumed to be normal with known variance. 

The main point of this method is to determine the 
compressed limit so that it maximises the difference of two 
values of the probability that the variable exceeds the 
compressed limit under the null and alternative hypothesis. 


This leads to the simple result that the desired compressed 
limit is half the value of the difference of the two values of 
the population mean mentioned above, divided by the 
population standard deviation. ‘The proof of the sample 
size being minimum, when the compressed limit is so deter- 
mined, is not rigorous but is based on the normal approxi- 
mation of the binomial distribution using the arcsine trans- 
formation. He also gives a graphical method of determining 
the number of observations and the critical number of 
defectives to test the hypothesis by using the compressed 
limit so obtained for fixed probabilities of the first and 
second kinds of errors. 

(This paper was written in Japanese, but the translated 
paper in English will be published by the Union of Japanese 
Scientists and Engineers in the near future.) 


(M. Siotani) 
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BEHRENS, W. U. (Kali-Chemie A.G., Hannover) 
A note on discriminatory analysis—In German 
Biom. Zeit. (1959) 1,-3-14 (3 refs., 4 tables, 2 figs.) 


This paper is concerned with the information supplied by 
inclusion of additional characters in problems of allocation 
and discrimination in the two-sample case. First the dis- 
criminatory power of one character is defined in terms of the 
standardised difference between sample means. Then it is 
shown that the additional discriminatory power supplied by 
a second character using a discriminant function is the 
greater, the higher is the absolute value of the coefficient 
of correlation between the two characters and the larger the 
difference between this correlation coefficient and the ratio 
of the second to the first standardised mean difference. 
Corresponding results follow for the addition of a third 
character. This statement, which is perplexing at first sight, 
is achieved by reformulating Mahalanobis’ D? as a function 
of the standardised sample mean differences and correlation 
coefficients between the characters. It seems to be valid 
only under the assumption (not explicitly stated) that there 
is no inter-character correlation between sample mean 
differences. The derived formulae together with a table 
make it possible to judge the usefulness of additional 
characters prior to an evaluation of the discriminant function. 


(R. Wette) 
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CLARINGBOLD, P. J. (University of Sydney, Australia) 
Multivariate Quantal Analysis—In English 
F. R. Statist. Soc. B (1958) 20, 398-405 (13 refs., 3 


In many biometrical investigations it is not possible to 
measure response on a continuous scale and recourse must 
be made to quantal scoring of signs exhibited by the animal. 
The author describes the canonical analysis of such multiple 
quantal variates. Considering g quantal variates (3, V2...Vq) 
the response vector is specified by o or 1 for each variate 
and the reaction pattern designated by the binary number it 
determines. ‘Transformations of the original variables are 
considered which include the possibility of a reduction in 
the number of variates. : 

Suppose an experiment is performed so that the three 
sources of variation may be examined: (i) experimental 
error, (ii) parametric or treatment effects and (iii) residual 
effects. ‘The dispersion matrices E and P are obtained by 
analysis of covariance of the variables into the first two 
sources of error. Determination of canonical variates, that 
is, estimation of the required matrix of transformation, 
enables the use of variables with unit standard deviation 
and which are mutually uncorrelated. ‘Successive rows of 
the transforming matrix obtained from the given equation 
correspond to the roots of the determinantal equation 
|P—®E |=o. Inthe absence of treatment effect in a canonical 
variate, the corresponding root should approach unity and 
may be tested by an approximation to the F-test. Thus 
the minimum number of variates required to describe 
significant treatment variation is determined and the treat- 
ment variation is maximised in the derived system. 


5.3 (6.4) 


5.8 (6.8) 


tables) 


One special case of importance is considered in which 
a score is obtained which gives the maximum dependence 
of response on dose. A numerical example of an oestrogen 
assay is analysed in a 34 factorial experiment in which there 
are three quantal variates. Two roots of the determinantal 
equation are significant and corresponding estimated 
standardised vectors ‘are obtained. In the experiment one 
of the factors was dose and the single variate which maximises 
dependence of response on dose is also obtained. 

The author mentions that some workers have reduced 
.a multivariate quantal problem to a univariate problem by 
pooling of classes and then applying standard methods of 
quantal analysis. He discusses the increase in precision 
resulting from the use of multivariate scores for response 
patterns, and stresses that, though the methods are somewhat 
laborious by manual computation, the availability of pro- 
grammes for automatic analysis greatly encourage their 
adoption. 


(N. W. Please) 


COHEN, L. (City College of New York) 
On mixed single sample experiments—In English 


5.2 (4.4) 


Ann. Math. Statist. (19538) 29, 947-971 (4 refs., 3 tables, 9 figs.) 


A mixed single sample (MSS) procedure is one which is 
derived from fixed sample (FS) procedures by choosing 
the sample size at random. For many problems use of 
MSS procedures will effect savings. For example, in testing 
two simple hypotheses the triple («, 8, 2), where « ( B) is the 
probability of Type I (II) error and is the expected number 
of observations, is used as the operating characteristic (OC) 
of an MSS test (MSST) and in several problems considered 
by the author (including that of testing two means of a 
normal distribution with variance known) for any FST 
whose OC is («, 8, n) there is an MSST whose («’, p’ n’) 
is such that «’ <a, B’<f, n'<n with strict inequality in at 
least one of these three. In the normal problem the good 
MSST?’s are randomisations of two FST’s, the first taking 
O observations and accepting one of the hypotheses always, 
the second being a Neyman-Pearson test based on taking a 
non-trivial number n,, of observations. In testing two 
rectangular distributions the first on (6,, 1 +0,), the second on 
(62, 1+.) the author shows that the FST’s are already the 
best one can do. 

The author characterises the admissible (in the OC 
sense) MSST’s for simple hypothesis testing by observing 
that an admissible MSST minimises for some positive 
a, b, c the quantity ax+bB-++cn from which it follows that 
an admissible MSST is a randomisation over v of FST’s 
of sample size v which are Bayes’ solutions within the class 
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of FST’s of size v and where the a priori distribution 
remains the same for all v. A prescription for computing 
the MSST is then given. Some examples are glven on how 
to obtain the best MSST for given «, f or for given («, 7) 


((B; 2). 


The author remarks upon certain peculiar features of 
these MSS procedures when extended to confidence interval 
estimation and composite hypotheses testing namely, that 
an MSS confidence procedure will cometimes choose a 
confidence interval of O length and O confidence coefficient 
and that an MSST in the composite testing problem will not 
have a power function which tends to 1 as the parameter 
moves far away from H,. This philosophical difficulty 
about choosing a sample of size O with positive probability 
also appears in the simple hypothesis testing problem. 


(J. Sacks) 


5.2 (2.5) 


On conformity of statistical data to a probability distribution—In Italian 
Ann. Fac. Econ. Com. Palermo (1958) 12, 209-219 (1 ref., 2 tables) 


In a statistical analysis the problem may arise of verifying 
the conformity of the observed data to the hypothesis that 
the empirical distribution of a particular phenomenon or 
character is according to a certain theoretical model. In 
this paper the author introduces some new statistics which 
aim at verifying the conformity of observed data to the 
hypothesis that the theoretical distribution is represented 
by the exponential Poisson law or by a new scheme of depend- 
ence recently constructed by Ferreri. ‘The proposed stat- 
istics are : 
ds as (x +2/x%+1)(1No+2/2? 241) 
and 
pa ra (%+3/x+1)(taNa+s/NatiNet2) 


for the Poisson law where n, = (NM*/r !)(e~™) and M is the 
arithmetic mean: N is the total frequency. For the Ferreri 
law we have : 


@, = [2—x(¢4—1)/1 —x(3—1)] [¢5—1/¢4—1]. 


These statistics use only some groups of the observa- 
tions, rather than the complete distribution, and therefore 
they can be used when the empirical distribution is trun- 
cated or censored. ‘The statistics have been applied to 
examine the degree of correspondence to the Poisson of a 
distribution of 400 squares of an “‘ emocitometer ” arranged 
according to the number of observed cells and the con- 


formity to the Ferreri distribution of the frequency distri- 
bution of 23,248 houses of Lione arranged according to the 
number of cancer deaths which occurred in each one during 
the period 1918-1938. 

In the paper are also indicated some simple criteria of 
estimation for the parameter M/ of the Poisson law and of 
the parameters A and d of Ferreri’s model. ‘These appear 
particularly interesting in the case of truncated or censored 
distributions as they, too, are based only on some of the 
groups of observations. 


(F. Galantino) 
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DOORNBOS, R. & PRINS, H. J. (Unilever N.V. Rotterdam & N. V. Philips’ Gloeilampenfabrieken, 


On slippage tests. III: Two distribution free slippage tests—In English 


Proc, Kon. Ned. Akad., Wetensch. (1958) 61 
Indag. Math. (1958) 20, 438-447 (6 refs., 2 tables) 


This paper is the last part of a series of four papers on tests 
for slippage of one of a group of variates, or one of a group 
of means of sets of variates, the first three of which were 
published in Indagationes Mathematicae (1956), 329-337 and 
(1958), 38-55. In these papers a number of general methods 
are applied to the cases of slippage tests for variates follow- 
ing various specified distributions with unknown parameter 
values and to some special cases. The foliowing methods 
were applied : 


(a) A simple and close approximation method to the true 
level of significance. 


(6) A method to obtain power functions in an approxi- 
mate way [derived from method (a)]. 


(c) A test for slippage of one of a set of means of groups 
of independent observations of unequal size. 


Method (a) was already introduced by Cochran. The 
present authors did not succeed in finding a general theory 
of slippage tests for one variate or one group mean. For this 
reason slippage tests for various cases are discussed in the 
papers. In the first papers, slippage tests for I’-variates (or 
variance estimates from normal samples), normal variates, 
Poisson, binomial and negative binomial varieties were 
discussed. ; 

The fourth and last paper gives first slippage tests for the 
largest and for the smallest of a number of column totals 
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5.5 (5.6) 
Eindhoven) 


obtained by the method of m-rankings. It contains also a 
table for the cases m=3(1)9 and the number of ranked 
objects R=2(1)10. 

The second test is a distribution-free k-sample slippage 
test for all sample sizes. The test statistic is obtained in 
the following way: all observations are ranked and each 
sample is tested with Wilcoxon’s two-sample test (one-sided) 
against the pooled other samples. The sample, which is 
associated with the smallest probability of exceedance of 
the values of Wilcoxon’s test statistic obtained, is said to 
have slipped if the probability concerned is smaller than the 
level of significance, chosen in advance, divided by the 
number of samples k. 

In a way similar to that used by D. van Dantzig [Indaga- 
tiones Math. (1951) 13, 1-8] to prove the consistency of 
Wilcoxon’s test against a certain class of alternative hypo- 
theses, the consistency of this test procedure against a related 
class. of alternatives is shown. 

Finally, the paper contains a table of critical values. for 
the test for Poisson variates described in the third part of 
this series of papers. ‘ 


(The authors) 


5.3 (5.6) 


On the Wilcoxon Test of comparing two samples—In Russian 
Bull. Acad. Sci.—Sér. Sci. Math. Astr. et Phys. (1958) 6, 611-614 


A statistical test is called consistent if it distinguishes two 
different probability distributions with the probability 
tending to 1, when the number of experiments increases 
to infinity. It is easy to verify that the Kolmogoroff- 
Smirnoff test is consistent. 

The paper contains an example proving that the Wilcoxon 
test is not consistent. 


(L. Kubik) 


Remark of the Regional Editor : 


From the following additional references it is known that 
the Wilcoxon test distinguishes between the distribution 


(oo) 
functions F(x) and. G(x) if and only if { F(x)dG(x) #3. 
: are AO. 

(a) D. van Dantzig (1951), On the consistency and the 
power of Wilcoxon’s two-sample test (Jndagationes 
Math., 13, 1-8). 

(b) E. L. Lehmann (1951), Consistency and Unbiased- 
ness of certain non-parametric tests (Ann. Math. 
Statist., 22, 165-180). 

(c) A. Rényi (1953), Neue Kriterien zum Vergleich 
zweier Stichproben (A Magyar Tudomdnyos Akad- 
émia Alk. Matematikai Int. Kézleményei, 2, 243-265). 


(A. Rényi) 


ISHII, G. (Ikuno Senior High School, Osaka) 
Kolmorgoroff-Smirnoff test in life test—In English 


5.6 (5.2) 


Ann. Inst. Statist. Math., Tokyo (1958) 10, 37-46 (10 refs.) 


The author considers certain distributions in connection with 
the Kolmorgoroff-Smirnoff distribution-free tests of good- 
ness of fit. Let there be two independent random variables 
with empiric distribution functions F,,(x) and G,,,(«) with a 
common continuous distribution function F(«). If the two 
sequences of observations from the empiric distribution 
functions are merged and ordered we have a sequence : 
Bas rent oe en Sey atta WE wan spOrtion. sof). this 
sequence is taken—as might occur if a life test was stopped 
at the mth failure—a theorem is proved on 


P{Gi(x) —F n(x) } <e, 


for —0O<x<zm. The method of proof is analagous to that 
used by Blackman in 1956 [Ann. Math. Statist., 27, 513-520]. 
The corresponding theorem relating to the opposite in- 
equality, 1.e., >—«’(Theorem 2) is followed by two theorems 
(3 and 4) dealing with P{F(x)—F,,(x)}<e for —o<x«<x, 
and P{F(x)—F,(x)}>-—e’ for the same restricted range of 
the variate x. These four theorems are concerned with a 
one-sided alternative to a test of a null hypothesis relating 
to the distributions. Theorems 3 and 4 have an application 
to the case where a life test is terminated at a prearranged 
time xp. 

The test of the null hypothesis for the two empiric 
distributions (Theorem 5) contemplates a two-sided alter- 
native and again the variable x is restricted to —O<x<2_. 
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In this paper only the case of k = 1 is treated, viz. : 
P{—’<G,(x) Ha) —<er. 


Two corollaries are given to Theorem 4 and another for 
Theorem 5. For all the theorems and the corollaries the 
expressions for the probability distributions considered are 
very complex. 


(Y. Suzaki) 


5-7 (4.7) 


Formulae for calculating the operating characteristic and the average 
sample number of some sequential tests—In English 
J. R. Statist. Soc. B (1958) 20, 379-386 (3 refs., 4 tables, 4 figs.) 


Work carried out by Page in 1954 (7. R. Statist. Soc. B. 16, 
36-139) in connection with a sequential sampling scheme 
enables approximations to be found to the operating charac- 
teristic and the average sample number. The present 
author takes up this work, uses the same method as Page 
but makes different assumptions, and obtains closer 
approximations. The method is applicable to a broad class 
of distribution functions for the random variable character- 
ising the unit of the sample, but calculations are performed 
for the normal case only. Comparisons are made between 
the proposed method of approximation, Page’s method and 
the true values. 


(F. N. David) 


LEHMANN, E. L. (University of California, Berkeley) 
Significance level and power—Jn English 


5.1 (5.2) 


Ann. Math. Statist. (1958) 29, 1167-1176 (13 refs., 2 figs.) 


The problems considered are (a) the general problem of 
choosing, for a significance test based on a given experiment, 
among possible significance levels « in the light of corres- 
ponding possible test powers 8 ; and (bd) the corresponding 
problem in the presence of a nuisance parameter v, par- 
ticularly in the case of multiparameter exponential distri- 
butions. As a basis for treating the latter case, a minimal 
essentially complete class of tests is characterised for the 
problem of a one-sided test of a hypothesis Hy: 0<0, 
(against #>6) on one of the parameters. Denoting 
by t the (complete) sufficient statistic for » when 6 = 0p, 
it is proposed to use conditional tests having conditional 
level a(t) at 0) and conditional power f(t) at 6,, where 6, 
is a chosen value of 6 for which it is important to reject Hp. 

Then the general problem (a) can be considered condition- 
ally on an observed value of t. Several possible approaches 
to the general problem (a) are described and suggested as 
more reasonable in many cases than use of a conventional 
level such as -o5 or -o1 for « or «(#). One of these is the 
choice of some constant k not very far from unity, and the 
choice of a test which satisfies « = k(1—8), or of a condi- 
tional test which satisfies «(t) = kA(1—f(t)). Several ex- 
amples are treated in detail. 


(A. Birnbaum) 
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Some non-parametric tests—In French 


5-6 (5.3) 


F. Soc. Statist. Paris (1958) 99, 254-284 (16 refs., 11 tables, 1 fig.) 


This paper summarises the work of recent years upon non- 
parametric two-sample tests which has been pubiished in 
statistical (English-language) journals. It is given in three 
parts: (1) Tests based upon the sum or average of ranks, 
(2) tests based upon the number and length of runs of 
alternatives, (3) tests based upon the signs of successive 
differences. ; 

In the first part the tests of Wilcoxon [Biometrics, 1945] 
Colin White [Biometrics, 1952] and Mann & Whitney 
[Ann. Math. Statist., 1947] are reviewed together with the 
k-sample test of Kruskal & Wallis [¥. Amer. Statist. 
Assoc., 1952]. Examples are given to illustrate these tests 
and the use of the special tables, extracts of which are 
included in the text. The second part summarises the tests 
proposed by Mood [Ann. Math. Statist., 1940], Mosteller 
[Ann. Math. Statist., 1941], Sved & Eisenhart [Ann. Math. 
Statist., 1943] and Wald & Wolfowitz [Ann. Math. 
Statist., 1940 and 1944]. Extracts of tables from the paper 
by Sved & Eisenhart are given and reference is made to 
the tables by Olmstead [Ann. Math. Statist., 1947]. 

In the third part, the author considers the tests proposed 
by Moore & Wallis [¥. Amer. Statist. Assoc., 1943] and 
Mann [Ann. Math. Statist., 1945]. In this connection the 


paper reviews not only tests concerning the number of 
positive and negative signs (of differences) but also the length 
of runs of differences of the same sign : the tables by Olm- 
stead are relevant in this latter instance. 


(D. Dugué) 


SALIB, S. R. (Faculty of Engineering, Ain Shams University) 


5.6 (5.1) 


Asymptotic power functions of some statistical rank order tests—JIn English 


Egyptian Statist. ¥. (1957) 1, 15-32 (5 refs.) 


Let X, Y, Z be three random variables and let F (t), Gd), 
and H(t) their distribution functions be such that they differ 
only in the value of the location parameter. The problem 
of testing the null hypothesis F (#) = G(t) = H(t) against 
the alternatives F(t)<G(t)<H(t) for all t¢ is considered. 

To test this hypothesis, the author considers three samples 
ayia Ty 2,0. 1, I) Jie 1,2, 6s.) Me and (RR SS oes 
mg and introduces the two statistics : 

1. The triplet statistic U,,,,,, defined by the relation 
mNMs Unnm, equals the number of possible triplets 
(x4, Vj, x) Such that xj,<9;<2,. 

2. The rank ordered statistic L,,,,,,, defined in the follow- 
ing way: Let (W,) t = 1, 2, ..., m-+ng-+n, be the ordered 
combined sample ; the statistic L,,,,.,, is the least number of 
nearest neighbour interchanges which turn (W,) into a 
sequence of three runs, with the first run formed by the 
elements drawn from the first population, the second run 
from the second population, and the third run from the third 
population. 

The author gives here the expectation and variance of 
the U,,n,,, under the null hypothesis and shows that its 
distribution is asymptotically normal. And, as he had 
previously obtained the expectation and the variance of 
Lynn, and had shown that its distribution is also asymp- 
totically normal, he compares the asymptotic efficiency of 
the two tests. 
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He then shows that, in the particular case where the 
distribution of the variables is rectangular, the relative 
efficiency of Lyn, with respect to Unnn is 2/15. And, in the 
case where the distribution of the variables is normal, he 
shows that the relative efficiency is 0°3. 

Lastly, the author notes that the number of comparisons 
used in the triplet statistic U has the order O(n*) while the 
number of pairs of comparisons in the L statistic is of the 
order O(n?). 


(J. B. Simaika) 


5.6 (1.3) 


A two-sample distribution free test for comparing variances—In English 


Biometrika (1958) 45, 544-549 (4 refs.) 


A two-sample distribution-free test is proposed for com- 
paring variances of two populations and a general formula 
is derived for its asymptotic efficiency with respect to the 
variance ratio (F-test). The proposed test presupposes 
knowledge about the relative locations of the two popula- 
tions. Where this information is not available it is suggested 
that deviations of the observations from the sample medians 
be considered and a modification of the test crition is pro- 
posed for this situation. It is shown that under certain 
regularity conditions the modified test is asymptotically 
distribution-free. 


(F. N. David) 


VINCZE, I. (Math. Inst., Hungarian Academy of Science, Budapest) 
Some two-dimensional distributions and limit-distributions in 


the theory of ordered sampling—In German 


5.6 (5.3) 


Mag. Tud. Akad. Kut. Mat. Intézet. Kézleményei (1957) 183-209 (15 refs.) 


In this paper the author determines the joint distributions 
and limiting distributions of some pairs of statistics. 

PEt fay. 10 on ON Hiya, s-, 7” be samples concerning 
the random variables € and y with common distribution 
function F(x). Let us denote further with F(x) and G(x) 
the empirical distribution functions of these samples. 

Let £9” and 79") denote in the following those smallest 


values of x, where the functions F,(x+0)—G,(x+0) and . 


| F,(« +0) —G,(« +0) |, respectively take their maxima. 

In theorem 1, the author determines with the aid of two 
combinatorical lemmas the joint distribution of the follow- 
ing two statistics : 


LF u(x) —Gr(x)], $F n( Eq” +0) +Ga( £,' -F0)] 


sup 
—o<x<0 

In theorem 2, the author considers the joint distribution 
of the statistics 

sup | F(x) —G,(x) |, 4LFn( no +0) + Gamo +0)]. 
— 0 <X%< 00 
This problem is reduced to a random walk problem. 

In theorems 3 and 4 the author determines the limiting 
distributions of the above distributions, in case noo. The 
marginal distributions of these limiting distributions are the 
well-known Smirnoff and Kolmogoroff-Smirnoff distri- 
butions, respectively relating to the one-sided and to the 
absolute deviation. 
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It is also shown that the limiting distribution of the 
variable £9") agrees with the common distribution of the 
variables € and 7. 

Finally in theorem 5 the limiting distribution of the 
quotient of the variables 


sup 
—ao<x<0a 
and 
LF n( Eo” +0) +Ga( 79 +0)] 
is determined. 
The author remarks that on the basis of these theorems, 
the tests for goodness of fit may be improved by using two 
statistics. 


(E. Csaiki) 


5.6 (5.3) 


Some simple unbiassed non-parametric tests for the two-sample problem 


with pairs of observations—-In German 
Metrika (1958) 1, 81-88 (9 refs.) 


A sample of size n consists of n pairs of numbers, each pair 
referring for example to the reaction of one person under 
two different conditions (treatments). For each pair, the 
difference of the two numbers is calculated. It is assumed 
that the differences are independent and that their distribu- 
tions are identical and continuous. The null-hypothesis is 
that these distributions are symmetric with regard to the 
value zero, corresponding to the assumption that both condi- 
tions (treatments) are equivalent. The alternative consists 
of all continuous distributions which are not symmetric with 
regard to the value zero. 

The test given below is uniformly unbiassed in the strict 
sense: the nu differences are ordered according to their 
absolute magnitude and then divided into r groups (r22), 
each group consisting of adjacent differences. Between 
every two adjacent groups one difference is deleted. ‘The 
deleted differences are neglected in the whole procedure : 
n, denotes the number of differences in the ith group, dy, 
the difference between the number of positive differences 
in the ith group and the number of negative differences in 
the ith group. Each test with the critical region Zg,|d;|p 2H; 
2;20, p= is uniformly unbiased in the strict sense. The 
proof of this theorem is given in a paper by E. Walter (so 
far unpublished). 


Subsequently, some special cases are considered : r=2, 


Ny =0, Nz=N—I, 213=82—p=1 is a test similar to the sign- | 


test, the only distinction consisting in the fact that the 
difference with the smallest absolute amount is neglected. 
T=2, m=n—Rk—-1, ny=k, 21=0, ¥2=p=1 is a so-called 
maximum-test: the null-hypothesis is rejected, if the k 
differences with the greatest absolute amount are all of the 
same sign. In this case, the size of the test is 2 *t, in- 
dependent of the sample-size. 

Some special features of uniformly unbiassed nonpara- 
metric tests are discussed. Special consideration is given 
to the maximum-test whose power is compared with the 
power of the Wilcoxon (Mann-Whitney)-test and the sign- 
test in the special case of differences resulting from normally 
distributed variables with different means and identical 
variances. 


(F. Ferschl) 


BARTON, D. E. & CASLEY, D. J. (University College, London) 
A quick estimate of the regression coefficient—In English 


Biometrika (1958) 45, 431-435 (12 refs., x fig.) 


In a bivariate sample (x,y;)i = 1, 2, ... , from a population 
in which the regression of y on x is linear the difference of 
the means of the & largest and the k smallest of the ordered x’s 
is divided into the difference of the means of the correspond- 
ing y’s to obtain an unbiased estimate of the slope of the 
linear regression. The efficiency of this estimate relative 
to the sample regression coefficient is found particularly for 
the normal bivariate distribution. In this latter case the 
estimate is shown to be the most efficient in large samples 
near the value k = n X 0:27 when the efficiency is 81 per 
cent. For small samples the first terms in an asymptotic 
expansion indicate an efficiency of between 70 per cent. 
and 80 per cent. for k between one-third and one-quarter of 
the sample size. 

The distribution of the estimate is shown to be sym- 
metrical and to tend to normality ; the leading terms of an 
asymptotic expansion for the second moment ratio, f., of 
this distribution are obtained. The method used is to find 
the movements of the estimate, first conditional on x, 
Xo, e254 0, and then over.x;, X2,...., X%, conditional on the 
upper and lower kth x-values. The expectation of these 
functions is obtained by the methods of David and Johnson, 
“Statistical Treatment of Censored Data,’ Biometrika 
(1954). 


(D. E. Barton) 


Ti. 


BEJAR, J. 
Contingency tables—In Spanish 
Trab. Estadist. (1958) 9, 85-101 


If quantitative observations of an aleatory variable are 
classified according to N distinct characteristics (quantitative 
or qualitative), the corresponding table of frequencies, which 
will be a double-entry table, is called a contingency table. 

This work deals with problems related to the hypothesis 
of independence of all the characteristics or groups of char- 
acteristics which form the contingency tables. It aims to 
give a general method in order to solve the independence 
problems of any number N of classifications or char- 
acteristics or p sub-groups of characteristics. The case of 
two characteristics is first considered and a general sketch is 
given of the method of testing the hypothesis of independ- 
ence of both classifications ; according to the “‘ Mathe- 
matical Methods of Statistics” of H. Cramér. This 
demonstration is generalised to the case of N classifications. 
The author then proceeds to consider the case of three char- 
acteristics. He puts the problem of testing the hypothesis 
that three of them should be independent according to the 
information given by a sample of n observations. The result 
is a generalisation of the former. It shall be coincident with 
it if one characteristic is taken without subdivision into 
intervals. The test ends with a distribution asymp- 
totically convergent to the x? distribution of Pearson with 
rst +r—s—t-+2 degrees of freedom; where 7, s, and ¢ are 
the number of intervals into which each characteristic is 
respectively subdivided. 


(Instituto de Investigaciones Estadisticas, Madrid) 


6.1 (4.3) 


6.9 (5.2) 


The next part of this paper considers the hypothesis of ~ 


the independence of one characteristic (7 intervals) as opposite 
to the other two (s and ¢ intervals). The test is performed 
by means of a x? test with rst —r—st+1 degrees of freedom : 
it coincides with the case of two characteristics if either 
t=1o0rs=1. This case of the independence of a char- 
acteristic as opposed to the other two jointly considered is 
solved (as stated by the author) in Introduction to the Theory 
of Statistics. The test in which every characteristic is 
independent from the other is considered for the general 
case of N characteristics. The variation of each is resolved 
into J,, Js, ...1, sets or groups. The test is completed with 
a x2 test with appropriate degrees of freedom. 

Finally, suppose that a group is arranged with the first 
h, characteristics, another with the second h, characteristics 
and so on until p groups are independent. ‘The test leads 
asymptotically to a x? distribution with J), J,, Js, ... ly tht 
lnytay eos dag 


- +h at Ih, yt2 + degrees of freedom. 


(S. Rios) 
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COX, D. R. (Birkbeck College, London) 


Two further applications of a model for binary regression—In English 


Biometrika (1958) 45, 562-565 (7 refs.) 


Two distinct problems are considered, both dealing with the 
analysis of quantal data, i.e. data in which each observation 
takes one of two possible forms, o and 1, say. (Thus o may 
be “non-defective ’’, 1 “‘ defective ’’.) 

The first problem is concerned with matched pairs. 
Let there be pairs of observations, the two observations 
in each pair corresponding say to two treatments A and Bs 
there being in general a correlation within pairs in the 
response. For each pair let the odds of the response 1 
for B be a constant, y, times the odds for A. The optimum 
test of the hypothesis ¢ = 1 (i.e. that A and B are equivalent) 
is shown to be one due to McNemar. Confidence intervals 
for % are obtained ; if two or more independent sets of data 
are available, the equality of their % values is easily tested. 

In the second problem, there are 2 independent trials 
and a sequence ;,..., p, is given. It is required to test 
whether the observations are consistent with the hypothesis 
that p,; is the probability of response 1 on the ith trial, 
71 = 1,...,n. Two statistics are introduced, one for examin- 
ing whether the p; are correctly spread, and the other for 
examining whether they are correctly centred. The two 
statistics can be combined to give a single x? test with two 
degrees of freedom. 


(D. R. Cox) 
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COX, D. R. (Birkbeck College, London) 


The regression analysis of binary sequences—In English 


6.8 (5.2) 


¥. R. Statist. Soc. B (1958), 20, 215-242 (discussion, 34 refs., 1 table) 


A sequence of o’s and 1’s is observed and it is suspected 
that the chance that a particular trial is a 1 depends on the 
value of one or more independent variables. Tests and 
estimates for such situations are considered. ‘The method 
is to assume a logistic dependence between the probability 
and the independent variable and, by considering an appro- 
priate conditional distribution, to obtain an optimum exact 
test and estimate for the slope parameter of the logistic 
relation. 'The problems considered differ from those of, 
for example, bioassay, in that it is assumed that the slope of 
the logistic line is the parameter of interest, the position of 
the line being a nuisance parameter. 

The following special cases are examined and related to 
earlier work : 


(a) 2X 2 contingency table ; 

(b) 2Xk contingency table with a quantitative variable 
associated with columns ; 

(c) test for serial order effect ; 

(d) test that the probability of occurrence is a specified 
value, against a trend-like alternative ; 

(e) test that the probability depends on the result of the 
previous trial ; 

(f) test that the probability depends on the total number 
of 1’s in ali previous trials. 


In the discussion, examples of applications calling for 
these and more difficult methods are given, together with 
general and particular theoretical comments on the methods 
used. 


(D. R. Cox) 
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6.8 (5.2) 


COX, D. R. (Report of Discussion on paper by) 
The regression analysis of binary sequences—In English 
J. R. Statist. Soc. B (1958) 20, 232-242. 


The discussion on this paper consists of eight verbal con- 
tributions, one of which was extended in writing, together 
with an additional contribution in writing received after the 
meeting and the author’s detailed reply. ‘Those who took 
part, in order of speaking, were :—C. A. B. Smith, G. H. 
Jowett, Miss V. R. Cane, D. V. Lindley (subsequently 
extended in writing), D. E. Barton, H. E. Daniels, G. M. 
Jenkins and G. A. Barnard. The written contribution was 
from A. Stuart. 

After commenting upon the possible applications of the 
paper, Dr Smith suggested that the difficulties of using 
the Gaussian or similar approximation to positively restricted 
variables might be overcome by using the distribution of 
the square of Gaussian variable with non-zero mean. 

Dr Jowett and Miss Cane also commented upon appli- 
cation of the methods of the paper and Mr Lindley con- 
tinued by offering a justification for the use of a conditional 
distribution in the argument of Section 4.1 of the paper. 
The justification rests upon the point that a desirable con- 
ditional test is a generally desirable test. 

Dr Barton drew attention to the surprisingly close 
parallel between the present paper and results in the theory 
of “ pure ranking ”’ as developments of paired comparisons. 
He also drew attention to earlier work on the Wald-Wolfo- 
witz runs test in connection with the binary Markoff chain 
by F. N. David [Biometrika (1947) 34, 335-339] as developed 
by Barton and David [Biometrika (1957) 44, 168-178] and 
Barton, David and Mallows [Biometrika (1958) 45, 166-177]. 
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DANIELS, H. E. & KENDALL, M. G. (Univ. of Birmingham and London School of Economics) 


6.8 (5.2) 


Professor Daniels showed that even in the case of random 
sampling from a finite population the logistic law is a natural 
approach, apart from its simplicity and mathematical con- 
venience, if the required approximate distribution is obtained 
by the method of steepest descents. Dr G. M. Jenkins 
contributed a simpler solution to the problem discussed 
in Section 6 of the paper—a solution which avoided any 
arbitrary assumption such as that of the logistic law. A 
criterion is presented which is similar to the non-circular 
serial covariance of lag one with a mean correction: raising 
the possibility that binary sequences, as considered in this 
paper, may be handled by a suitable modification of the 
criterion used in the autocorrelation analysis of continuous 
variables. 

Professor Barnard pointed out some cautionary points 
to be observed in using the logistic—after remarking upon 
the value and simplicity of the logit transformation. Mr 
Stuart’s written contribution, inter alia, draws attention 
to some additional points regarding the optimum character 
of the Wilcoxon Test for location problems. 

The author’s detailed reply offered further comments 
on points raised by the discussion. 


(Wm. R. Buckland) 


6.2 (2.2) 


Short proof of Dr. Harley’s theorem on the correlation coefficient—In English 


Biometrika (1958) 45, 571-572 (4 refs.) 


It has been shown by B. I. Harley that the mean value in 
repeated sampling of the inverse sine of the sample cor- 
relation coefficient is equal to the inverse sine of the popula- 
tion correlation coefficient. This result is proved by Daniels 
and Kendall using a combination of geometry and concord- 
ance. ‘Their method of proof is much shorter than that 
given by Harley and demonstrates effectively the basic 
reason for the existence of the relationship. 


(F. N. David) 
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HODGES, J. L., Jr. (University of California, Berkeley) 
18 itting the logistic by maximum likelihood—In English 
Biometrics (1958) 14, 453-461 (6 refs., 3 tables) 


Fitting a logistic response curve to quantal data by maximum 
likelihood usually entails an iterative process. The “ transfer 
method ” is presented as a means for eliminating or shorten- 
ing the iterative process, depending on the degree of accuracy 
required. Common practice in logistic curve fitting involves 
sight fitting a straight line to points plotted on a logistic 
scale or on logit paper. Estimates of the parameters of the 
logistic curve from such a line are subjective and not well 
defined, serving only as the first step in the iterative process. 
The transfer method is based on the maximum likelihood 
estimates of the minimal sufficient statistics for the para- 
meters, namely, the total response, estimated by the sum of 
the observed (positive) responses in several groups subject 
to different stimuli levels, and the total ‘‘ moment ’’, 
estimated by the sum of the products of the stimuli and their 
corresponding observed responses. 

The method is graphical and requires a straight line 
sight fitted to the original data as before. The aim is to 
determine a new set of responses, corresponding to the 
stimuli, which are collinear and which preserve the observed 
total response and total moment. This is accomplished by 
a series of transfers of responses from points above the line 
to points below the line, subject to the foregoing conditions. 
A new sight line is drawn and the process repeated until the 
points are collinear to graphical accuracy. Computers, 
tables or charts are not required. If graphical accuracy is 
sufficient, the parameters of the logistic curve may be 
determined from this final line. Otherwise these values may 
be used to initiate a single cycle of iteration which is claimed 
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6.7 (0.2) 


to result in accuracy usually attained only after several cycles. 
The transfer method assumes that the observations follow 
the logistic binomial model. The fitted line may be used to 
carry out a chi-square test (d.f. = no. levels—2) of this 
assumption using the original unadjusted data. The transfer 
method is easiest to apply in cases where an equal number 
of trials is made at each level of stimulus and where the 
levels of stimulus are equally spaced. Examples are given 
for this case and for one in which the number of trials is not 
the same for each level of stimulus. The method is still 
applicable, though less simple, when successive levels of 
stimulus are not equally spaced. A third example is given 
to show the ease with which this method handles the zero 
or 100 per cent. response. 

This paper deals only with the technique of maximum 


likelihood estimation and not with its desirability for fitting — 


the logistic. Where the maximum likelihood estimates do 
not exist for certain values of response, it is also impossible 
to find transfers which preserve the total moment. The 
transfer method requires the existence of simple sufficient 
statistics and thus is not directly applicable to the fitting of 
the integrated normal. Nevertheless it provides a good set 
of starting values for the iterative process in this case. It 
is noted that probit instead of logit paper may be used with 
somewhat less accuracy and that methods such as the 
Spearman-K4rber technique may be used in lieu of the 
initial sight fitted line. Reference is made to extant and 
forthcoming charts and tables by Berkson which in part 
make this an interim method. (M. F. Peck) 


6.6 (10.0) 


A restatement of the conditions for identifiability in complete systems 


of linear difference equations—In English 
Metroeconomica (1958) 10, 182-190 (3 references) 


In Cowles Monograph No. 10, Koopmans, Rubin and 
Leipnik stated and proved certain theorems relating to 
identifiability in complete systems of linear difference 
equations where those equations were subject to disturb- 
ances, i.e. had stochastic error terms. ‘The author, in using 
these theorems, has found that not only may the concept 
of identification as stated not be precisely the one required 
but that the theorems tend to lead to a certain amount of 
confusion. He suggests that some attention might usefully 
be given to conditions which make a set of coefficients either 
“just identifiable ” or ‘‘ unidentifiable ”’. 

After giving particulars of his notation the author states 
his problem as being concerned with methods of testing 
whether the specification of any given model is sufficiently 
narrow to allow identification of the particular (but not 
fully known) structure or the identification of a certain set of 
its characteristics. ‘The general sharpening of the concept 
of identification used in this paper rests on the definition 
of structures which are observationally equivalent. ‘Two 
structures are “‘observationally equivalent” if, for the 
observed variables of the set of initial variables and for the 
reigning joint distribution of exogenous variables, the 
structures give rise to the same (conditional) distribution 
of the observable variables. 


The author then states that the linearity of the model 
leads to the expectation that observational equivalence is 
closely related to invariance under non-singular linear 
transformations. After stating a lemma dealing with a 
necessary and sufficient condition, the author gives his own 
theorem—implied by this lemma and the definition of 
identifiability. He then states that the simplest application 
is the case in which the model specifies By in Ax; = Boys + 
T'(z,;) = u; to be the identity matrix, J, and proceeds to 
apply his theorem to recursive models. Further sections of 
the paper investigate the identifiability of rows of the 
matrix A and the effect of homogeneous restrictions to single 
rows. 


(Wm. R. Buckland) 
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LEVIS, V. (Faculty of Economics, University Institute, Venice) 6.1 (-.-) 
Mean Square Regression linear spaces—In Italian 
R. C. Mat. Univ. Roma (1958) 17, 423-438 


Let us consider an n-dimensional random variable, repre- 
sented on the Euclidean space S,,; the author studies the 
regression of an r-dimensional component variable (referred 
to a subspace S, of S,) on the complementary (n—r)- 
dimensional component variable (referred to the comple- 
mentary space S,_,), and particularly the mean-square 
regression linear space (‘‘spazio lineare dei minimi quad- 
rati’’) S,_,, defined as the best linear estimate, according 
to the principles of least squares, of this regression. 

The author calculates, by means of concise matrix 
algebra, the coefficients of the equations of si and then 
proves the following propositions : 


I. the space Se depends upon the space at infinity 
Se, of S, only; i.e. is invariant for the group of 
affine transformations under which S~, does not 
vary ; 

2. confining ourselves to spaces through the mean 
point or referred to mean points as origin, the relation 
between S, and Chee is reciprocal (i.e. the mean 


. . * . 
square regression space relative to S,_, is S,) and 
precisely is connected with the polar system of an 
ellipsoid of S;,. 


(G. Dall’ Aglio) 
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MOUSTAFA, M. D. (Institute of Statistics, Cairo University) 6.2 (6.1) 
A note on Partial Correlation—In English 


Egyptian Statist. F. (1957) 1, 1-5 (3 refs.) 


The author considers the case of three correlated variables 
and studies the partial correlation coefficient between two 
of them. 

He defines the partial correlation coefficient as the 
ordinary coefficient of correlation between the residuals 
obtained by subtracting from the two variables their true 
regression function on the third variable. 

The form of the regression function is generally unknown 
and has to be assumed. The author proves that: If at least 
one of the regression functions is rightly assumed except 
for a constant, then the coefficient of correlation obtained 
by using the assumed regression function has the same sign 
and is in modulus less or equal to the true partial correlation 
coefficient. 

In particular, if at least one of the regression functions 
is linear, he proves that the well known formula for the 
coefficient of partial correlation gives a value which has the 
same sign and is in absolute value less or equal to the true 
coefficient of partial correlation. 


(J. B. Simaika) 


RIZZI, A. (University of Roma) 6.0 (6.4) 
Cograduation tables in the three-dimensional case—In Italian 
Statistica (1959) 19, 72-78 (4 refs., 8 tables, 1 fig.) 


The purpose of this paper is to extend the cograduation 
tables, constructed by Salvemini, to the case of triple 
variates. ‘Tables of this type occur in the calculus of the 
indexes of homophily. 

In Italian usage, these indexes are special measures of 
concordance between statistical phenomena. 

To overcome difficulties in regard to multiple variates, 
the author makes use of the technique used by Salvemini 
in connection with cograduation tables. In addition, he 
shows a logical scheme relative to a particular cograduation 
table. The plan can be adapted to automatic computation ti) 
on electronic machines. } 


(V. Amato) 
87 
ROY, S. N. & GNANADESIKAN, R. (Institute of Statistics, North Carolina) 6.1 (4.8) 
A note on “ Further contributions to multivariate confidence bounds ’—In English 


Biometrika (1958) 45, 581 


The authors refer to a previous paper of theirs in Biometrika 
(1957) 44, 399-410. In this previous paper dealing with the 
general linear hypothesis, confidence statements were 
derived and similar statements were obtained involving the 
truncation of columns of a given matrix. 

The authors point out in this present note that there is 
another type of truncation which arises in discussing the 
uninvariate general linear hypothesis. They illustrate this 
by considering the problem of testing the equality of several 
treatment effects in ordinary analysis of variance. 


(F. N. David) 


SANKARAN, M. (Presidency College, Madras) 


On Nair’s transformation of the correlation coefficient—In English 


Biometrika (1958) 45, 567-571 (5 refs., 2 tables) 


Nair suggested a statistic x — (r—p)/(1— pr) where r is the 
sample, and p the population correlation coefficients. It is 
shown that x and the inverse sine of x are both asymptotically 
normally distributed for m large and that the inverse sine is 
more quickly normal than is x. It is also shown that the 
sine of x can be used to evaluate the probability integral of r. 
A confidence interval for p is also obtained. The author 
remarks (p. 568) “...it is not clear that the inverse sine 
transformation of x has any practical advantage over the 
2-transformation, indeed it is not appropriate in some 
situations where 2 is useful.” ’ 


(F. N. David) 
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6.2 (2.2) 


6.6 (4.8) 


The estimation of economic relationships using instrumental variables—In English 


Econometrica (1958) 26, 393-415 (16 refs.) 


This paper applies the use of instrumental variables, first 
suggested by Reiersel in 1941 [Econometrica 9, 1-23], to 
the more general case in which relationships are not exact : 
the ideal economic variables assumed generated by a set of 
dynamic stochastic relationships and the actual economic 
time series subject to random disturbances or measurement 
errors. The author considers the problem of determining 
the coefficients for a single relationship where the residuals 
are not autocorrelated. The main results are concerned 
with asymptotic properties of the estimates and to this end 
the author first considers the asymptotic covariance matrix. 
With this he proceeds to discuss the problem of choice 
between the possible instrumental variables in order to 
obtain the reduced set necessary to enable the method to be 
used. At the end of this section is a short discussion of the 
advantage of adding one new instrumental variable—stating 
that eliminating one from the reduced set cannot improve 
the asymptotic variance matrix. 

In the fifth and sixth sections of his paper, the author 
draws a parallel between his analysis and the canonical 
correlation approach of Hotelling [Biometrika (1936) 28, 
321-335] and of an alternative derivation in the form of a 
minimax approach similar to that suggested by Durbin 
[Rev. Int. Statist. Inst. (1954) 22, 23-54]. The paper then 
goes on to consider the asymptotic distribution of the 
smallest characteristic root and other matters prior to 
deriving two confidence regions for the unknown coefficients 
and a brief report of some lengthy (unpublished) calculations 


on the bias of the asymptotic approximation: the biases 
in the estimates tend to be large when the relationship is 
almost unidentified. The last two sections of the paper 
give a comparison with alternative methods of estimation : 
(i) Least Squares, and (ii) Limited-Information maximum- 
likelihood methods. 

The author hopes subsequently to publish some of the 
large amount of work that has been carried out to test in 
practice the method outlined in this paper. From previous 
results and work some tentative conclusions are offered: 

(i) The use of the instrumental variable method will 
produce convenient estimates of coefficients even 
when large measurement errors are apparent. 

(ii) Better results may appear to be obtained if an 
attempt is made to reduce measurement errors— 
without introducing autocorrelation in the residuals. 

(iii) Large measurement errors in predetermined vari- 
ables are likely to cause biased estimates if the Limited- 
Information maximum_-likelihood method is used. 
The Least Squares method also produces strongly 
biased estimates when there are large measurement 
errors, or when some of the independent variables 
are not predetermined and the relationship is 
confluent. 

(iv) Large numbers of instrumental variables may not 
improve the accuracy of the estimates. 


(Wm. R. Buckland) 


SHEPS, M. C. (Harvard University Medical School) 


6.9 (4.3) 


An examination of some methods of comparing several rates or proportions—In English 


Biometrics (1959) 15, 87-97 (10 refs., 4 tables) 


The author considers various statistics used for quantitative 
comparisons of relative frequencies in the fields of public 
health and medicine. Studies like those of the association 
between smoking and lung cancer, between congenital 
malformations and exposure to radiation, or between coronary 
disease and obesity require some kind of quantitative measure 
of the “change in risk”. Rates, such as mortality rates, 
have been compared by calculating the ratio of two mortality 
rates, the ratio of two survival rates, and the difference 
between the two mortality or survival rates. 

The author analyses these statistics in terms of several 
models. ‘The model defined as Model 1 applies to the study 
of different populations with mutually independent risks. 
He maintains that for this model the difference in rates has 
an advantage over a ratio because the numerical value of the 
ratio depends on whether mortality rate or survival rate (one 
minus the mortality rate) is selected for comparison, while 
the difference is unaffected by this selection. 

An alternative model (Model 2) is suggested for the 
common situation where the factor being investigated is 
assumed to modify the reaction of the population in one 
direction only. Under this model the expected rates in the 
several samples are not independent. As a result, the 
difference is a product of two probabilities, the ratio between 
one set of rates is compounded of several probabilities, while 
the ratio between the complementary set of rates is the 
estimate of a single probability. When a deleterious effect 
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is expected, the ratio of success rates is an estimate of the 
diminished chance of ‘‘ success’ or survival; and when a 
beneficial effect is expected, the diminished risk of failure 
is estimated by the ratio of ‘‘ failure ’’ or mortality rates. 

The author derives the maximum likelihood estimates of 
the parameters, and their variances, for Model 2. ‘The 
aforementioned statistics are then expressed in terms of these 
estimates. He proposes the use of a new statistic, an esti- 
mate of one of the parameters in Model 2. In terms of the 
smoking and lung cancer problem this parameter is the 
following probability of “risk”?: Given that a smoker 
belongs to a category of persons who would otherwise escape 
lung cancer, he now suffers the risk that he will fail to do so. 

This method is also applied to the poliomyelitis vaccina- 
tion problem where the Model 2 parameter is the increase 
in protection given by vaccination rather than an increase in 
risk. ‘The model is incapable of handling a situation that 
could occur in such studies, namely, of a vaccine actually 
producing the disease in some percentage of the cases. 

Several different modifications of the model are indicated 
for the case with more than two samples. 


(R. J. Taylor) 


6.1 (10.7) 


Covariances of least-squares estimates when residuals are correlated—In English 


Ann. Math. Statist. (1958) 29, 1251-1256 (4 refs.) 


The author is concerned here with certain problems which 
arise in linear regressions when the residuals are not in- 
dependent. Earlier work has given least squares estimates 
of the regression parameters and the residual variance, 
together with some asymptotic properties of these estimates ; 
references to this work are given. ‘The contribution of this 
paper is to obtain explicit formulae for the elements of the 
covariance matrix and the expected value of the residual 
variance when a specific variance matrix is assumed for the 
residuals ; analytically simple bounds for the elements of the 
covariance matrix are also obtained. The residual variance 
matrix which has been assumed is one in which all elements 
an equal distance off the main diagonal are equal. 

Two special cases have been worked out: (1) The true 
response is a linear combination of two non-stochastic 
variables, and an odd number of observations are made. 
(2) The true response is a linear combination of trigono- 
metric functions which are themselves known functions of 
non-stochastic variables. The general formulae obtained 
earlier for the covariance matrix and expected value of the 
residual variance are specialised for each of these two cases. 

This paper deals only with the situation where the 
residual variance matrix is known. It is pointed out that, 
while this is frequently not the case, the least squares 
estimates of the regression parameters and residual variance 


may still be obtained. Since this desirable property is not 
enjoyed by maximum-likelihood estimates, the author sug- 
gests the inadvisability of research in such a direction. 


(D. H. Shaffer) 
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TINTNER, G. (Iowa State College) 


6.7 (0.2) 


A new method for estimating the logistic function—In German 


Metrika (1958) 1, 154-157 (9 refs.) 


Time series concerning developments in the field of 
economics or population, often show a trend which can be 
described in a satisfactory manner only by an exponential- 
function of the time ¢ rather than by a polynomial in ¢. 
The fact that the polynomials become infinite for tco 
makes them inappropriate especially for the purpose of 
extrapolation. For many purposes, the logistic function is 
well suited for describing the trend: y, = k/1-+be-*. The 
methods commonly used for curve-fitting (method of least 
squares or maximum likelihood) lead to difficulties in this 
case, as the resulting equations for the unknown parameters 
are not linear. 

A way out consists in considering the differential equation 
of the logistic function or the difference equation of its 
reciprocal z; = 1/y;. These equations are linear relation- 
ships between magnitudes which can be calculated, in prin- 
ciple, from the data originally given. The differential equa- 
1 dys 
ye at’ 
which, however, can only approximately be determined 
from the data given. The difference equation of z; = 1/y;, 
on the other side, is linear in 2;,, and 2. Therefore, the 
common methods for estimating linear functions can be 
applied under the condition that the values of y; are given 
for equidistant times. With this method the development 
of the population in Sweden is treated as an example. In 


tion of the logistic function uses the rate of growth, 
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treating the differential or difference equation instead of the 
original function, one parameter—the constant of integra- 
tion—is lost. In the case of the logistic function, this 
parameter can be estimated by a formula of Rhodes. 


(F. Ferschl) 


6.8 (6.7) 


Multinomially grouped response times for the quantal response bioassay—In English 


Biometrics (1958) 14, 462-488 (15 refs., 3 tables) 


In work with quantal response data experimental approaches 
may be classified according to the method of observation 
employed. Suppose doses &,, &, ...,&m are administered 
tO 14, Np, ...)%m individuals respectively. Let the number 
of individuals responding to dose é; be observed at times 
Ti <Tig< ... <Tin, Which are pre-specified. The data are 
then said to be multinomially grouped according to response 
times. 

The aim of the authors’ analysis is to estimate a relation 
between 9, é, and 7, where p is the proportion of the popula- 
tion that responds to an amount of stimulus € in a period of 
time r. Particular attention is given to relations for which 
a transform y of p is a linear function (additive model) of 
transforms x and ¢ of é and 7, or a linear function plus an 
interaction term (non-additive model). A detailed discussion 
is given for these two models when y is the logit transform 
of p. Considerable attention is given to the problem of 
choosing metameters. ‘Two procedures for choosing these 
are presented for the above particular models, and an 
example is given in which the dose metameter 1s log dose 
and the time metameter is the reciprocal of the square root 
of time. The troublesome case when no response is obtained 
in some time intervals is discussed and an adjustment is 
recommended. 

Since the response numbers are accumulated for the 
successive times within a dose, the observed proportions 


are not independent within doses. This introduces some 
complexity into the analysis which is at least partly resolved 
by use of the minimum modified chi-square method of 
estimation of the coefficients of the relationship. Computa- 
tional details are given for this method of estimation. An 
appendix contains the derivation of the estimators of con- 
stants in the functional form by both the method of maximum 
likelihood and the method of minimum modified chi-square. 
A derivation of the variances of the estimates obtained from 
the estimated relation is given in the appendix for the 
additive and non-additive models. 


(C. P. Quesenberry) 
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COX, C. P. (National Institute for Research in Dairying, Univ. of Reading) 
Experiments with two treatments per experimental unit in the 
presence of an individual covariate—In English 


7.6 (9.7) 


Biometrics (1958) 14, 499-512 (13 refs., 9 tables) 


The methods available for eliminating temporal or spatial 
trends in experimental material from treatment comparisons 
include stratification in the design and covariance on the 
trend variable. They eliminate the mean trend over a set 
of experimental units, and are likely to be inadequate in 
some situations, particularly when within-block treatment 
comparisons are desired. The author has previously given 
a method for removing individual curvatures and isolating 
treatment comparisons in Latin square designs. In the 
present work this design restriction is removed, and a 
general solution obtained, giving within-individual treatment 
comparisons free from general curvature effects (assumed 
independent of treatment effects). The case considered is 
that in which there are two treatments per individual or 
block. Time trends are removed by the use of equally 
spaced observations. Upon making necessary changes, one 
may render the solution applicable to the more general 
problem of unequally spaced observations. 

For each individual the time is divided into three periods 
consisting respectively of t,, s, and t, equal and successive 
intervals. ‘The two treatments, of which the main effects 
are to be compared, are administered, one during each of the 
first and last periods. One independent observation is 
obtained in each interval, yielding t, observations in the 
first period and f, in the third. No observations are required 
in the intervening period of s intervals. In biological work 
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HALL, W. J. (University of North Carolina, Chapel Hill) 
Most economical decision rules—In English 


this may serve as a conditioning period ; in some cases it 
may be possible to take s = o. 

For the complete experiment, a mixed model is used with 
fixed treatment effects and random individual and inter- 
action effects. In application to an individual the treatment 
and interaction effects are confounded. The least squares 
technique is used to isolate a comparison involving the 
difference between the two treatment effects and the differ- 
ence of corresponding interaction effects. The sum of 
squares for treatments and interactions, as well as those for 
regressions up to the desired order are obtained. The 
establishment of orthogonality properties is indicated and 
the analysis of variance set up. 

Tests of significance in the complete experiment, are 
discussed and the precision of the method is examined. 
Reference is made to the question of allocating values to 
t, and t, so as to minimise the contribution of the un- 
explained variation to the variance affecting the treatment 
comparisons. A solution is given for the case of linear 
gradients, explicit solutions for other cases not being as yet 
available. Consideration is given to the testing of the 
adequacy of the mathematical model and the examination 
of the validity of the associated assumptions. The procedure 
is illustrated in a numerical example. 


(W. A. Glenn) 


7.7 (1.8) 


Ann. Math. Statist. (1958) 29, 1079-1094 (15 refs., 1 fig.) 


The problem considered is to determine optimum decision 
rules for what is often referred to as the discrimination 
problem, i.e. the problem of choosing one of m mutually 
exclusive alternative hypotheses A;. The decision procedure 
then is to subdivide the sample space into m mutually 
exclusive regions R; and if the sample point falls in R;, the 
decision A; is made (or hypothesis A; accepted). The 
author denotes by p;;(D) the probability that the decision 
procedure leads to decision A; when in fact A, is true. 
The matrix ||p,,(D)|| characterises a decision rule or dis- 
criminant function. 

The author then gives two definitions of a most eco- 
nomical decision rule. A decision rule Dy is most economical 
with respect to the vector a = (a4, ag, ... Gn) if 

pi(Dn) 24: iG = 1, 2,..- n) 
and WN is the least integer m for which this inequality is 
achieved by any decision rules based on samples of size n. 
In other words such a decision ensures that the probability 
of rejecting A; when A; is true has prescribed bounds and 
it does so for the smallest possible sample size. 

A decision rule Dy is most economical with respect to 
the matrix || 6|| if p4(D)S By @, 7 = 1 ... m) and N is the least 
integer n for which this inequality is achieved by any decision 
rule. ‘These two definitions are extensions of the Neyman- 
Pearson Type 1 and Type 2 Error concepts respectively. 


However, the author then sets up special types of loss 
functions, which lead to a solution of the problem of optimum 
discriminant functions from results of Wald’s minimax 
theory. Conditions are given that a most economical 
decision rule exists. "These are given both in fairly general 
terms and for families of distributions of exponential type. 

Some examples are given where the random variables 
are normal and the hypothesis A; specifies that their mean 
lies in some interval (6,1, @;). As is intuitively to be ex- 
pected, the decision regions are of the form accept Aj, if 
Ci_1 S%;<; (o known) or c’;_, St <c’ (¢ unknown) where ¢ is 
the usual Student ratio. The c’s are to be found by solving 
iteratively a series of equations in nm+1 unknowns involving 
the normal distribution function. The equations for the 
c’s are similar except that the normal distribution is replaced 
by the non-central ¢-distribution. Since a different set of 
these equations must be solved for each a (or|| f ||) and set of 
6’s, the exact optimum solution is not too well suited as 
yet for practical use. 


(D. G. Chapman) 
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MATHER, K. & JONES, R. MORLEY (University of Birmingham) 


73°C) 


Interaction of genotype and environment in continuous variation: II. Analysis—In English 


Biometrics (1958) 14, 489-498 (4 refs., 2 tables) 


This paper examines the requirements for successful separa- 
tion of the effects of genotype-environment interaction from 
other components of variation, and for the estimation of these 
components by comparisons and analyses of the variances 
and covariances which are obtainable in experiments. 

When polygenically determined continuous variation is 
involved, important comparisons can be achieved by raising 
true-breeding lines and their F, in the same range of environ- 
ment where these environments are so distinguished from 
one another that differences between them are large in their 
effects when compared with variation within them. Then 
genetic parameters can be estimated as well as the environ- 
ment and interaction parameters. 

When statements of genetical constitution must be based 
on knowledge of ancestry, on observations or on assumptions 
of linkage and selection, as with continuous variation, then 
families must be grown over a comparable range of environ- 
ment. If the range of environments is large enough so 
that the sampling variation is negligible in respect of geno- 
type and of environments then the differences between their 
mean phenotypes will reflect only their differences in average 
genetical constitution. Variation between individuals will 
reflect the differences amongst the environments of indi- 
viduals and also interactions. If the families are distributed 
at random over the environment, this paper furnishes the 
variances of the segregating generations. 
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If the individuals of all families are not distributed at 
random over a common range of environments, but if each 
family is raised in one or more compact groups, which are 
themselves distributed over a wider range of environments, 
the environment of differences affecting the variation of 
individuals within the plot are then different from those 
affecting the variation in plot and family means. In this 
case a method is given to separate the effects of environ- 
mental differences into those within plots and those between 
plot means. The interaction components of variation 
between plots may be isolated, but analysis of variation 
within plots is again possible only in terms of transforms. 

A previous paper, Part I—Description, was published 
by the authors [ Biometrics, 14, 343-359]. 


(C. Y. Kramer) 
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DARROCH, J. N. (Department of Mathematics, University of Cape Town) 
I. estimation of a closed population—In English 


The multiple recapture census. 
Biometrika (1958) 45, 343-359 (16 refs.) 


In this paper the author considered multiple-recapture 
census for which the number of samples s is usually con- 
sidered fixed. In a few cases s is allowed to be variable. 
The sample sizes are regarded both as constants, and as 
binomial variables. The investigation is restricted to a closed 
population in which there is neither departure due to death 
or emigration, nor augmentation due to birth or immigration. 
The author will take account of both departure and augmen- 
tation in a second paper. 

Two alternative models are constructed: Model (A), 
a multinomial density and claimed to be new, is applicable 
when the effort put into the catching of every sample is 
fixed before the experimenter starts ; Model (B), a generalised 
hypergeometric density which forms the basis of most of 
the previous investigations, is used when the experimenter is 
determined to catch a specified number of individuals at 
each sample. 

The two models are shown to lead to the same estimate 
of the population, n, and to the same asymptotic variance of 
this estimate. Model (A) is introduced as it is capable of 
generalisation which (B) is unable to accommodate. Con- 
fidence interval estimation under Model (B) and a special 
case of one individual per capture are investigated. 


99 


DE LUCIA, L. (University of Rene} 


8.3 (4.3) 


In the case where s is a random variable, a model (A,) 
is constructed and finally the author proves that for a certain 
class of estimates of n, those derived form (A) and (A,) 
by the method of maximum likelihood are asymptoueally 
the best. 


(C. Crouse) 


8.1 (4.4) 


Problems of sampling techniques in market research—estimation of proportions—IJn I eho 
Quad. Ass. Ital. Studi Mercato, no. 3, Giuffré, Milano (1958) (19 refs., 12 tables) 


In this report of 128 pages are collected the lessons held, 
during a course in Market Research, on the estimation of 
proportions in random sampling of attributes, which, as 
the author puts in evidence, is one of the most frequent 
cases in Marketing problems. The author exposes in an 
elementary form the basic notions of this field in the case of 
random sampling without replacement, emphasising the 
role of the problem of repeated trials and of the connected 
random variables. 

The first part is devoted to simple sampling. The author 
exposes first the classical methods and formulas for esti- 
mating the proportion and the difference between two 
proportions, and for evaluating their standard errors. ‘Then, 
reporting some results given in a previous paper [Riv. Ital. 
di Econ., Demog. e Stat. (1957) 11, 123-142], he deals with 
the estimate of proportions with regard to a part of the 
universe. In this study, as he shows, it is possible to include 
the case of a sample in which are present sampling units 
that cannot be classified. In the last section the methods 
of estimating the sampling variance by means of the sampling 
data are given. 

In the second part, devoted to stratified samples, the 
author repeats the topics of the first one. There is added 
the treatment of two most usual methods of dividing the 


total sample size among the strata: the partition being 
proportional to strata sizes, together with optimum 
allocation. Some fictitious examples, extensively worked out, 
accompany the development of the theory. 


(G. Dall’ Aglio) 
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ELTETO, 6. & SARKADI, K. (Math. Inst., Hungarian Academy of Sciences, Budapest) 
On the statistical problems arising in microbiological examination of 


tomato purées—In Hungarian 


8.8 (1.8) 


Mag. Tud. Akad. Kut. Mat. Intézet. Kézleményei (1958) 219-226 (2 refs., 2 tables) 


The authors deal with questions connected with the follow- 
ing problem : 

A lot of tomato purée will be accepted by the consumer 
if the Howard-count found in the sample taken out from the 
lot in question does not exceed a fixed limit. What limit 
is to be used by the producer on the occasion of a previous 
quality control for the purpose of retaining the bad lots ? 
It is assumed to be legitimate to regard the distribution of 
the Howard-count as a Normal one. The variance is 
assumed to be known from empirical data. The error arising 
from the heteroscedasticity can be eliminated by a homo- 
scedastic transformation. 


The problem is discussed according to the following | 


three principles :. 


1. It is prescribed that the joint probability of accepting 
in the first and rejecting in the second examination should 
not be greater than 1—f. The solution is given by formula : 


x = t+206-(4/1— 8). 


o”? denotes the variance of the Howard-count, ¢ and x 
the critical value in the final and in the preliminary sampling, 
respectively, and ®—1(u) denotes the inverse of the standard- 
ised normal distribution function 


u 
e/2 dz, 
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GARNER, N. R. (Aerojet-General Corporation, Azura, California) 


Curtailed sampling for variables—In English 


2. Let w, be the loss caused by the wrong decision in 
the first examination if we retain a lot which would be 
accepted, and wy the loss if we dispatch a lot which will be 
rejected. In this case the minimax solution is given by 
formula 


Ww 
x = t+200-} ( 4 : 
Wiy+We 
3. It is assumed that the expected value of the Howard- 
count of a lot is itself a random variable varying in the lots 
transported from the same producer according to the normal 
probability law. The problem can be solved numerically 


by the aid of the table of the two-dimensional normal 
distribution. 


(I. Vincze) 


8.8 (5.7) 


JF. Amer. Statist. Ass. (1958) 53, 862-867 (6 refs., 1 table) 


The usefulness of curtailed sampling is evident when testing 
is costly, when each test takes a long time, or when failure 
might constitute a hazard. Thus, if under a single sampling 
plan specifying sample size n the eventual rejection decision 
can be foretold with certainty after only r<n observations, 
it will be advantageous not to sample the remaining (n —7) 
items. Curtailed sampling has long been prescribed for 
attribute single-sampling ; for variables single-sampling, 
curtailed procedures have been suggested in the past but no 
optimum condition for this kind of sampling procedure has 
as yet been suggested. : 

This paper supplies such an optimum condition for 
curtailed sampling of variables in the one-sided case. 
condition is optimum in the sense that, if sampling is curtailed 
and the “‘ reject ” decision made after only r <n observations, 
then there exist (n—r) possible remaining observations such 
that rejection would “‘ just barely ” be indicated on the basis 
of the full sample of n. 

Suppose the lot is to be rejected if X+ks>L,, where 
X and s are the mean and standard deviation for the entire 
sample of n observations and k and L,, are specified by the 
sampling plan. Then curtail the sampling and reject if 
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The - 


X,+ks,>Ly, where X, and s, are the mean and standard 
deviation for the curtailed sample of the first r observations 


PN eaten) Tian (n—r)(r—-1) 
n—1I nrk® 
If the lot is to be rejected when X —ks<L,, curtail and. 
reject if X,—,s,<Ly. 


(H. T. David) 


GHOSH, S. P. (Calcutta University) 


A note on stratified random sampling with multiple characters—In English 


8.1 (8.7) 


Bull. Calcutta Statist. Ass. (1958) 8, 81-90 (3 refs., 5 tables) 


The optimum allocation ina case of stratified random sampling 
for multiple characters is defined in this note as that allo- 
cation which for a fixed total sample size minimises the 
generalised variance, namely, the determinant of the variance 
covariance matrix (;;) of the estimates of the population 
values. The equations providing the optimum values of 
strata allocations as functions of all the variances and co- 
variances are complicated. However, for obtaining the 
solutions numerically, an iterative procedure is suggested. 
As a first approximation, the individual minimum values of 
each o;; is obtained and substituted in a set of equations to 
give the values of n; for the different strata. The values of 
n;, So obtained, were substituted in another set of equations 
to obtain the second set of approximations of o;;. This 
process is repeated till stable values of n; are obtained. The 
condition under which the proposed optimum allocation 
will be identical with the average of Neyman’s allocations 
for the different characters is also considered. 

A numerical example with six strata and three characters 
is used to illustrate the proposed allocation with a view to 
demonstrate empirically the rapidity of the convergence 
of the iterative process. Optimum strata allocations have 
been obtained 

(i) assuming the correlation between two characters to 
be the same for all strata and 

(ii) assuming that there is no correlation between the 

characters. 
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HAMMOND, L. C., PRITCHETT, W. L. (Univ. of Florida, Gainesville) & CHEW, V. (North Carolina 


Soil Sampling in Relation to Soil Heterogeneity—In English 


Using this example, the variances of the estimates of the 
population means under the optimum allocation have been 
compared with those obtained by 


(i) proportional allocation, 
(ii) average of Neyman’s allocation, and 


(iii) allocation arrived at by minimising the sum of the 
ratios of the variances for the characters to their 
respective optimum variances. 


(M. N. Murthy) 


8.1 (8.7) 
State College, Raleigh) 


Proc. Soil Sct. Soc. Amer. (1958) 22, 548-552 (12 refs., 10 tables) 


After reviewing recent statistical studies that have been 
made on the problems of soil sampling, the paper presents 
analyses of soil samples taken from cultivated fields (tobacco 
and corn), pastures and orchards (oranges). ‘The seven 
locations ranged in size from 2 to 35 acres and included 6 
representative soil types in Florida. ‘The samples were 
analysed for moisture, pH, and extractable potassium, 
calcium, magnesium and phosphorus. 

A three-stage sampling scheme was adopted. The 
number of first-stage units in the field and in the sample 
varied with the size of the field, but the number of second- 
stage units was 9 in all cases, of which 2 were randomly 
selected. Two one-inch diameter cores were taken from 
each second-stage unit. There was no compositing of cores. 

Applying the usual analysis of variance techniques, 
but taking the finiteness of the populations of sampling 
units into consideration, variance components were calcu- 
lated for each of the fields and for each of the variables 
measured. The efficiency of the multi-stage sampling 
scheme relative to simple random sampling was found to be 
high (generally over 70 per cent.) for ‘“‘ macrouniform ” 
fields, i.e. for fields where the variation among the larger 
units is small compared to variation among the subunits. 
For such fields, the administratively more convenient 
multi-stage sampling is recommended over simple random 


sampling. 
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The paper concludes with some cost consideration and 
optimum allocation of the sampling units. For a simple 
assumed cost function, the larger number of samples required 
for multi-stage sampling scheme to achieve a given precision 
actually cost only about one-third that for simple random 
sampling. 


(V. Chew) 
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MOSHMAN, J. (Corp. for Econ. and Indust. Res., Arlington, U.S.A.) 


8.3. (-.-) 


A method for selecting the size of the initial porapicn in Stein’s two-sample procedure—ZJn English 
Ann. Math. Statist. (1958) 29, 1271-1275 (5 refs., 1 table) 


In 1945, Stein published a paper [Ann. Math. Statist., 16, 
243] giving a two-sample procedure whose power was 
independent of the variance. No specification was given for 
the choice of the size of the initial sample but in 1953 a 
paper by Seelbinder [Ann. Math. Statist., 24, 640] suggested 
that it should be selected to minimise the expectation of the 
total sample. An efficient choice will hold the expectation 
of the sample small-and will further reduce the probability 
of an extremely large total sample. This note by the present 
author explores the problem and shows that an upper per- 
centage point of the distribution of the total sample size, 
when used in conjunction. with the expectation, is a guide 
to the efficient choice of size for the initial sample. The 
procedure outlined is a rapidly calculable function based 


upon ordinary x? (distribution) tables and is a function of 


the population variance. 

_ If an initial sample (Vo) is taken and an estimate s? 
(based upon N,—1 degrees of freedom) of the population 
variance is computed, then the total sample size, N, is: 
max [(s?/z)-++1, No] where z>0. It follows that P(N = No) 
= P(ns?/o?) = P[x?/n<nNpo2/o7] where n = Ny—1. Similarly 
PIN = m),=P[n(m—1)2/0? < x?/n<nmz/o?]. These two 
expressions facilitate the calculation of a given percentage 
point, N,, of the distribution of N : in fact, 

N, =-max{No, A~}[100 pth percentage point of x°(m)/n] +1} 
where A = 2/07. 
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STANGE, K. (Technische Universitit Berlin) 


An objective rule is given for the choice of No and it is 
suggested that p would be chosen independently from non- 
statistical grounds. An example is given (A = 0-1) together 
with a brief discussion of some subsidiary questions. 


(Wm. R. Buckland) 


8.8 (11.2) 


The graphic treatment of plans for quantitative inspection—In German 


Metrika (1958) 1, 111-129 (3 refs., 10 figs.) 


The paper deals with sampling procedures for variables. 
The use of variables inspection (measurement) instead of 
attribute inspection (classification by gauges) is of advantage, 
if the cost of measurement is nearly equal to the cost of 
classification, if the cost of inspection is very high anyway 
(for example due to the destruction of the inspected items) 
or if the fraction defective is very small. 

The quality characteristic x is assumed to be normally 
distributed with parameters », o. An item is regarded to 
be defective if x exceeds some given value JT. The quality 
of a lot is judged by calculating the test statistic z = X+ks 
and accepting the lot if z<T, rejecting if z>7T. 2 has a 
mean of = »+ko and variance 

2 2 

n a 
(n = sample size). It is approximately normally distributed 
for n>10. In order to obtain the acceptance probability, 
one has to consider u,= T—é /o?, which is a linear function 
of the magnitude u = T—p/o which is relevant for the 
fraction of “‘ defective’ items (items with x>T7T). The 
coefficients of this linear function depend only on n and k. 
Therefore, the operating characteristic (OC) curve, giving 
the acceptance probability as a function of the fraction 
defective, is a straight line if plotted in double-probability 


paper (i.e., a probability paper, which has probability scales 
on both ey Therefore, if two points of the OC curve 


are given (for example producer’s safe point, producer’s 
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risk and consumer’s safe point and consumer’s risk), it 
can be plotted easily on this double-probability paper. 
Starting from the OC curve on the double-probability 
paper, the sample size m and the constant k can easily be 
determined by graphical methods. The use of this graphical 
method makes the use of any table superfluous. 

In order to condense the information, given by a sequence 
of samples (for instance in order to decide whether the 
inspection level should be reduced or tightened), usually 
the statistic 


T—x' /. sg n—I 
t’ => 7 (# SS Dik, = Zs) 
Tye mn—I j= 
is used. It is proposed to use instead of ¢’ the statistic 
Bees RG TP i3 
Et, (= ie): 
™ j= 


Both, #’ and #, are asymptotically enlaaee distributed with 
the same variance and therefore equally well suited for this 
purpose. 7 can, however, be easily determined by graphical 
methods, as the magnitudes #, are already given. 

Similar results concerning attribute inspection will be 
given ina future paper. In this paper, also a simple equation 
will-be given for the relation 1,,/n, which compares the 
sample sizes in the case of measurement (,,) and classifica- 


tion (”;). 


(J. Pfanzagl) 


STANGE, K. (Technische Universitit, Berlin) 
Belt Sampling—In German 


Metrika (1958) 1, 177-223 (5 refs., 20 figs.) 


Many production processes, especially in the chemical 
industry, are of a continuous nature. For supply and delivery 
a belt system is used. If samples are taken from the belt, 
in many cases the results of subsequent samples are corre- 
lated. The correlation is a decreasing function of the 
distance between the samples. The shortest distance for 
which adjacent samples are non-correlated is called the 
region of correlation. 

The accuracy of the sample mean & depends on the 
variance of the individual values, the sample size n, the way 
in which the sample is taken (random sample, stratified 
sample or systematic sample), the length of the belt L, the 
distance between two adjacent samples / and the region of 
correlation A. If /> the correlation is of minor importance 


and the variance calculated from the sample can be used for, 


estimating the variance of @. If /<A the correlation may be 
of great importance. A correlation integral is introduced, 
by means of which a critical distance of adjacent samples 
and a critical sample size n, is defined. 

If the quality of a given belt-section of length L is to be 
inspected by systematic sampling (i.e. sampling at equal 
distances), the variance decreases very rapidly (in some 
instances to the order of 1/n?) if the sample size is greater 
thann,. ‘The variance within the belt section can practically 
be eliminated by taking a sample of size 27,. 

If the purpose of the inspection is to judge the process 
itself, samples of size L are taken from different belt sections. 
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STEVENS, W. L. (University of Sao Paulo, Brazil) 


8.9 (10.7) 


In this case the variance converges for nto the value of 
the correlation integral. For the correlation functions 
usually observed this limiting value is reached for n = 2n,. 
Therefore the accuracy cannot be increased by taking 
samples larger than n = 2n,,. Thus in both cases there is no 
essential gain in accuracy by taking samples of a size greater 
than 2n,. 

These results also hold true if the subject under investi- 
gation are mean values for longer periods. If on the other 
hand the purpose is to regulate the process, values for specific 
time points are needed. In this case the accuracy can be 
increased by taking samples at shorter distances. 

The results of this inquiry are of practical importance 
only if the correlation function is stable during time. Whether 
this is the case can not yet be decided and many empirical 
investigations will be necessary to clarify this question. 


(A. Adam) 


8.1 (4.3) 


Sampling Without Replacement with Probability Proportional to Size—In English 


J. R. Statist. Soc. B (1958) 26, 393-397 (1 ref.) 


This paper puts forward an alternative to Yates & Grundy’s 
{%. R. Statist. Soc. B., 1953] extension of Horvitz and 
Thompson’s [¥. Amer. Statist. Assoc., 1952] result and it is 
not self-explanatory without reference to both of the earlier 
papers. A finite population of N units, each of which pos- 
sesses a measurable characteristic y, is supposed. Itis desired 
to estimate the mean value () of this characteristic in the 
finite population using the y-values of a sample of N units. 
In drawing the sample the probability of selecting a particular 
unit is proportional to a second measure denoted by « which 
measures the size of the unit. The x-measures of each and 
every member of the finite population are supposed known. 
Stevens considers the case where the elements of the 
population fall into groups having the same size, there being 
N, elements of common size x; in the ith group and N;2n 
for alli. He proposes that n groups be selected with replace- 
ment, the probability of selection of the ith group being 
proportional to the total sizes, i.e., N; x, in the group and 
that if the ith group is selected ¢; times then ¢; units are 
selected from that group without replacement. On this 
basis he obtains an unbiased estimator called 7, of the total 
Nu which is formally the same as that for the case of sampling 
with replacement. He also obtains the variance of this 
estimator and an unbiased estimator of this variance. 
Stevens considers that Yates & Grundy’s method is 
“hardly practicable ” for n>2 but that his method can be 
easily applied for an n of (say) ten with, if necessary, certain 
modifications (e.g., pooling the size groups) ‘which he 
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discussed and claims will his 
results. 

If p,; is the probability of selecting any given element 
in the ith group and if r denotes the result of dividing the 
y-value of an element by its probability of selection (viz. 
p; if its size-group is the ith) then 7 is the mean of the n 
sample values of r. Since the probability of selecting t 
elements in the ith group is the binomial probability of 
selecting ¢ in a sample of m independent events with constant 
probability N.p;, and since the conditional probability of 
selecting a given element from amongst those in the 7th 
group is t/N;, then the overall probability of selecting that 
element is pj. 


n7, called o* (R), is obtained in a similar but more lengthy 


not materially affect 


‘manner by summing the conditional probabilities over the 


multinomial probabilities of selecting given numbers of units 


Hence, E(n7) = nNyu. The occurrence of | 


from a specified set of groups, thus obtaining the probabilities — 


of selecting a given pair of elements. The unbiased 
estimator of this variance is simply obtained by observing 
that the mean value of the sum of squares of the 2 r-values 


about 7 differs from (n—1)o?(R)/n by a simple multiple of | 
XN ,o% (where o% is the variance of the y-values within the | 


ith group). 


An unbiased estimator of this is found by means | 


of a weighted sum of the within-group variances of the 


r-values. 


(D. E. Barton) 
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WEILER, H. (Comm. Sci. Ind. Res. Organisation, Prospect, N.S.Wales) 
Mean charts of constant false alarm risk for non-normal distributions—Jn English 
Aust. F. Appl. Sci. (1958) 9, 326-331 (5 refs., 4 figs.) 


In quality control methods, the mean p of a population is 
usually controlled by entering the means @ of random samples 
of fixed size m into a control chart with limits »+Bo/4/n, 
where o is the standard deviation of the population and 
B = 3 or 3:09. If the population is normal, the probability 
that % falls outside these limits is 0-001 when B = 3:09, 
so that such an event mostly indicates a change in population 
mean. However, many populations occurring in practice 
are non-normal, so that false alarms (% outside the limits for 
unchanged ,) may occur more often than expected on the 
basis of normality. : 

In this paper, the author establishes a relation between n 
and B such that control charts will provide the same pro- 
tection against false alarm when non-normal distributions 
are encountered. It is shown how the width of the control 
limits can be much reduced as the size n of the sample 
increases. 

By using an expansion in terms of the standard normal 
distribution ¢ (x) and its derivatives for the probability p 
that @>u-+ fo/1/n, an expression is obtained for the false 
alarm risk py; = py t+A¢"(B)/6./n, an upper barrier to the 
values of p for a given set of populations whose skewness 
does not too often exceed A. Here py is the false alarm risk for 
a standard normal population, and ¢”(B) the second derivative 
of the standard normal distribution at x = B. 
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8.9 (2.7) 


Curves of constant false alarm risk per sample are given, 
in which B is plotted against n for fixed values of A = 0, 
0°5, I, 1°5, 2. Further curves of constant false alarm risk 
Pi = 0°00I, 0:002, 0-003, 0004 per sample, with B plotted 
against n are provided for A = 1, 2. 

Lastly, curves of constant false alarm risk per article 
tested are given, B being plotted against m for A = 0, 0°5, 
1, 2. The average number of articles tested before raising 
an alarm is found to be about 4000, 1877, 1176, 690 respect- 
ively in each of these cases. 


(J. Gani) 


8.7 (1.8) 


On the choice of the best amongst three normal populations with known variances—In English 
Biometrika (1958) 45, 436-446 (9 refs., 4 tables, 2 figs.) 


This is a problem in statistical decision functions. It is 
given that there are three normal populations with known 
variances. Decisions have to be taken in (i) the problem 
of selecting the population with the largest mean given a 
probability of making a wrong decision, and in (ii) the 
problem of determining the smallest sample sizes required 
to select the population with the largest mean—for given 
values of differences between population means and for 
given values of decision probabilities. The statistic which 
is used is the standardised difference between the two 
largest sample means. Tables whereby the significance of 
the statistic under (i) can be judged are not given but can 
be obtained from the author if desired. An alternative 
procedure based on the difference between the largest sample 
mean and the grand mean is compared with the procedure 
proposed under (i) and it is suggested that this alternative 
procedure is not so good. 

Under (ii), the planning of an experiment, some assump- 
tion is necessary about differences between the population 
means. An abacus is given whereby, for fixed values of the 
first and second kind of error, the equations of condition 
for the sample sizes can be solved. This solution was pro- 
posed by Bechhofer: an optimum solution is achieved by 
repeated application of Bechhofer’s method. 


(F. N. David) 
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ATIQULLAH, M. (University of Dacca, Pakistan) 


9.1 (-.-) 


On configurations and non-isomorphism of some incomplete block designs—In English 


Sankhya (1958) 20, 227-248 (20 refs.) 


This paper discusses the existence and isomorphism of, and 
relations between, certain classes of Balanced Incomplete 
Block (BIB) and Partially Balanced Incomplete Block 
(PBIB) designs. Let W, denote the class of balanced 
incomplete block designs with parameters v = n(n—1)/2, 
b = n(n+1)/2, Rk = n—1, ry = n-+1, A= 2 and 7, the class 
of two associate partially balanced: incomplete block de- 
signs with parameters v = b = n(n—1)/2, r= k = n—1, 
Ay = 2, A. = I, ty = (n—2)(n—3)/2, Nz — 2n—4, 
Pi, = (n—4)(n—5)/2. The following theorems are proved : 


Theorem 1:A necessary and sufficient condition for 
W,, to exist for n= 4 is that 7, exists. 

Theorem 2: If n24 and n(n—3)/2 odd, a necessary 
condition for 7, to exist is that (n—1) is a perfect 
square. 


The non-existence of designs W, for n = 6, 9, 13, 14° 
and n = 18 follows from these two theorems. 

The problem of obtaining all non-isomorphic solutions 
of dn for various values of n is then taken up. It is shown 
that for m = 4, just one and for m = 5, just five solutions 
exist : for n = 8 two new solutions are obtained. 

It is further shown that a Triangular Singly Linked Block 
design I, which is a two-associate PBIB design with para- 
EGA. Ul = Gamay 1) == Pe (G30) Pa 
Ag = OF ny =-2n—4, nN, = (n—2)(n—3)/2 and pi, = n—2 
has a unique solution for any n. 
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(East Malling Research Station, Maidstone) 


Lastly, two methods of generating two-associate PBIB 
designs D, with parameters v = n(n+1)/2, 6 = n(n—1)/2, 
kR=n+1, 7 =n—r, Ay = 2, Ay = 1, m = (n—1)(n—2)/2, 
Ny = 2n—2, pi, = (n—3)(n—4)/2 from a symmetric BIB 
design with r= +1 and \ = 2 are described and it is 
shown that the two methods may lead to non-isomorphic 
solutions for D,. 


(J. Roy) 


9.1 (9.7) 


Families of designs for two successive experiments—In English 
Ann. Math. Statist. (1958) 29, 1063-1078 (9 refs., 7 tables) 


The author discusses all known families of O: PP designs 
with two associate classes, where a design of type O: PP 
is one in which the first set of treatments is arranged in 
randomised blocks (a type O, for orthogonal, design) and 
the second set of treatments is arranged as a partially balanced 
incomplete block design with two associate classes (a type P 
design) with respect to both the blocks and the first set of 
treatments (referred to as rows and columns, respectively). 
The O:PP designs are classified according to the 
association scheme of the P-type design since, for practical 
purposes, the O : PP designs have been constrained to have 
the same associate classes with respect to rows and columns. 
The majority of the O : PP designs are group divisible with 
the remainder having triangular and Latin square association 
schemes. ‘The author states that it appears unlikely that 
any other association schemes will provide O : PP designs. 
The author shows that not every pair of group divisible 
P-type designs can give rise to an O: PP design. However, 
there are 56 pairs that do and the parameters for the resulting 
56 families of group divisible O: PP designs are listed. 
Also listed are the parameters for one family of triangular 
designs and one family of Latin square designs. The 
numerical values of the parameters are tabulated for all the 
“useful”? designs (designs having more than 2 replicates or 
treatments, not more than 30 replicates, treatments, rows or 
columns and not more than 150 plots in all). Most of the 


useful group divisible designs occur in 3 of the 56 families, 
namely, certain subsets of the designs which are either 
singular with respect to both rows and columns, or singular 
with respect to one and either regular or semi-regular with 
respect to the other. The remainder of the useful group 
divisible designs are confined to 10 of the remaining 53 
families. For each of the 43 families of group divisible 
designs that have no useful designs the author has tabulated 
the smallest designs that are theoretically possible. In all, 
the author has listed the parameters of 127 useful group 
divisible designs, eight triangular designs and two Latin 
square designs. 

A catalogue of the useful group divisible and triangular 
designs has been prepared and is available at East Malling 
Research Station. General methods of analysis of the 
O: PP designs have been given by the same author in one 
of the references. 


(John W. Wilkinson) 
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GOMES, F. P. & GUIMARAES, R. F. (S40 Paulo Univ. & Comp. Paulista de Estradas de Ferro, Brazil) 
Joint analysis of experiments in complete randomised blocks with 


some common treatments—In English 
Biometrics (1958) 14, 521-526 (6 refs., 4 tables) 


The paper deals with the intra-block analysis of a group of 
experiments in complete randomised blocks, where some 
treatments, but not all, are common to the whole group. 
The common treatments are called by that name and the 
others are called ‘“‘ regular treatments”. The “ regular 
treatments ” belonging to one experiment are not common 
to any other experiment. 

The design studied for the analysis is a special case of 
intra- and inter-group balanced incomplete blocks, proposed 
earlier by Rao in a paper listed in the references. The 
whole set of experiments is regarded as one trial with incom- 
plete blocks, and under the assumption of equal variances a 
joint analysis is carried out. 

Least squares estimates of the treatment means on 
Model I are obtained. For the common treatments these are 
simply the arithmetic means ; however, the estimate of the 
mean of a regular treatment in any experiment has to be 
adjusted by subtracting from the arithmetic mean the 
difference between the mean of the common treatments in 
that experiment and the overall mean of the common treat- 
ments. 

The estimates of the variances of various possible con- 
trasts among the estimated treatment means are also given. 
These variance estimates are 

(1) V(Y;) of the estimated contrast Y, between the 

means of two regular treatments in the same experi- 
ment. 
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(2) V(¥Y2) of the estimated contrast Y, between the 
means of two regular treatments in different experi- 
ments. 

(3) V(Y3) of the estimated contrast Y, between the 


means of two common treatments based on all 
observations. 


(4) V(Y.) of the estimated contrast Y, between the 
means of a common and a regular treatment. 


_ In the case whan only one treatment is common 
V(Y,) = V(Y,) and the contrasts Y, and Y, are similar. 
As a test of these contrasts Tukey’s studentised range test 
is proposed with k = no. of regular treatments+1. The 
other type of contrast which may be considered is that of (2). 
Here again Tukey’s test may be used with k = number of 
regular treatments. When only contrasts between the 
common treatment and a regular one are to be considered, 
Dunnett’s procedure with least significant difference = 
t./V(Y.), where ¢ refers to the multivariate Student t-distri- 
bution, is to be preferred. 

The analysis of variance of the design is given. 
illustrate the various procedures, a group of two experi- 
ments on different species of Eucalypts, each with five 
replications and ten treatments, one of which was common to 
both the experiments, is analysed. 


(R. L. Chaddha) 


9-5 (6.8) 


Application of quantal-response theory to cross-comparison of 


taste-stimuli intensities—IJn English 


Biometrics (1948) 14, 548-557 (8 refs., 3 tables, 2 figs.) 


This article concerns some taste testing experiments which 
were carried out using the method of paired comparisons. 
The purpose of these experiments was to illustrate the use 
of this technique when applied to a quantal response bio- 
assay problem. In particular, various concentrations of 
sucrose were compared with various concentrations of sodium 
chloride, citric acid and quinine sulphate. These represent 
the four primary tastes. 

The underlying model is based on Thurstone’s “ phi- 
log-gamma ”’ hypothesis in which the use of the cumulative 
normal for the sigmoid is replaced by the logistic function. 
This is more in keeping with the assumption of the Fechner 
law that the strength of the sensory process is proportional 
to the logarithm of the stimulus. If C,; and C; represent 
two different concentrations of a stimulus, then the model is : 


Y = In[P/(1—P)] = a +8 log (Ci/Ci) 
where P is the proportion of the time that C, is judged 
stronger, say, than C;. E(a) = 0 and £ is a parameter of 
sensitivity. ‘The use of the Weber fraction in measuring 
acuity is also discussed and, in fact, is equal to 
F = antilog [(In, 3)/8]—1- 


For a given substance, E(Y) = 0 occurs at Ge GF lt; 
however, valid comparisons are being made between two 
qualitatively different stimuli, the concentration ratio at 


Y =o will be the taste-intensity ratio of the pure stimuli, 
weight for weight. This ratio is the main objective of this 
type of experiment. The appropriateness of the model can 
be studied by goodness-of-fit tests. 

For the case where different substances are being com- 
pared, the model is then generalised to 
Y = In[P/(a—P)] = a+8; log C:+ 

+B; log C;+yiz (log C;) (log C;) 
where C; and C; are now concentrations of different sub- 
stances. If the response to both substances is the same, 
then y = 0, 8; = f; and the model would become the same 
as the original one. 

The experimental results are tabulated and discussed. 
These include estimation of the relative intensity for the 
four substances which are then compared with earlier, less 
statistically rigorous, methods. 


(J. E. Jackson) 
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MAURICE, RITA J. (University College, London) 


Selection of the population with the largest mean when comparisons can 


be made only in pairs—In English 


9.5 (1.8) 


Biometrika (1958) 45, 581-586 (7 refs., 4 tables, 1 fig.) 


Where it is required to select the best (indicated by the 
largest mean) of a number of normal populations of known 
equal variance and these populations can be compared only 

in pairs, the usual procedure is to use the standard experi- 
mental design for incomplete blocks with two plots per 
block. A possible alternative procedure would be to use a 
procedure of successive elimination, following the procedure 
familiar in cup-ties and tournaments. The simplest case 
in which the latter procedure may be used is for the selection 
of the best one out of four and its performance in this situa- 
tion is examined. 

The effect of the blocks on the experimental results may 
be constant or subject to random variation. If the block 
effects are constant, the cup-tie procedure appears to attain 
a given probability of correct selection more cheaply than 
the incomplete blocks procedure. This held for a number 
of probability levels, but the general case could not be 
proved. Assuming that the loss resulting from an incorrect 
choice is proportional to the difference between the largest 
and the chosen mean, the minimax principle may be used 
to determine the sample sizes of both procedures. The 
expected loss of the cup-tie procedure is then less than that 
of the incomplete blocks procedure whatever the size of the 
difference between the largest and the other means. 

If the block effect is a random variable, the incomplete 
blocks procedure achieves a given probability of a correct 
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SCHEFFE, H. (University of California, Berkeley) 
Experiments with mixtures—In English 


choice more cheaply than if the block effects are constant. 
The analysis of the cup-tie procedure is not affected by 
the change in the model. If the variance of the block effect 
(o;) is sufficiently small the incomplete blocks procedure 
will be cheaper than the cup-tie procedure. For the prob- 
ability levels selected, oj must be less than or equal to 0-80? 
(where o? is the variance of the populations) for this to hold 
in all cases. When the minimax principle is used to deter- 
mine the sample sizes the incomplete blocks procedure is 
always better if oj = 0:50”, and always worse if o; = 207. 
For oj; = o? and oj = 1°50" which procedure results in the 
smaller expected loss depends upon the value of the differ- 
ence between the largest and the other means. The com- 
parisons are made assuming o;/o0? is known. In practice 
this ratio may have to be estimated and the incomplete 
blocks procedure will then compare less favourably with the 
cup-tie procedure. 


(R. J. Maurice) 


9-7 (9.3) 


Y. R. Statist. Soc. B (1958) 20, 344-360 (2 refs., 4 figs.) 


The author discusses the design of experiments in which 
the response depends on the properties of the components 
present in a mixture and not on the amount of the mixture. 
The factor space for a g-component mixture with the experi- 
ment arranged to have m-+1 equally spaced proportions 
(Qs Sufiiy AiiGencca , 1) of the components, is a simplex lattice 
denoted by {g, m}. All the possible mixtures correspond to 
points on the lattice. 

He considers a {g, n} polynomial of degree n in q variables 
having the same number of coefficients, of the various 
products of the variables, as the number of points on a 
{q, n} lattice. The author conjectures, proving it for the 
case of n = 1, 2 and 3, that the values of the polynomial 
can be assigned arbitrarily on the lattice and that its values 
on the simplex are then uniquely determined. Combina- 
tions of the coefficients of the polynomial may be interpreted 
in terms of responses to the mixtures ; linear coefficients 
to pure components, quadratic coefficients to synergism 
between two components and so on. The author further 
shows that for n = 1, 2 and 3, the relation of the coefficients 
to the responses at points on the lattice is related to ortho- 
gonal polynomials of the responses and considers significance 
tests of the fit of the polynomials. Suggestions for the allo- 
cation of numbers of observation at each point are made 
from consideration of the variances of predicted responses 


based on the order of the polynomial selected. Modifications 
of simplex lattice designs are suggested when it is necessary 
to keep the properties of one or more of the components 
below certain limits. 

In a mathematical appendix the following points are 
examined : 


(i) The number of points on a {g, m} lattice. 
(ii) 
(iii) 


(iv) 


The number of coefficients in a {g, n} polynomial. 
The canonical form of a {q, 2} polynomial for n = 1, 
2 and 3. 

The calculation of certain orthogonal polynomials 
and the maxima of their squares. 

(v) A statistic for testing several correlated differences 
between observed and predicted responses. 


Estimation of the regression coefficients when 
m>n and one observation is made at each point. 


(vi) 


(N. W. Please) 


SHAH, B. V. (University of Bombay) 


A note on orthogonality of experimental designs—In English 


Bull. Calcutta Statist. Ass. (1958) 8, 73-80 (no refs.) 


This paper adopts the definition of an orthogonal design as 
one which provides mutually uncorrelated estimates of any 
two estimable parameters belonging to two different groups. 
BY? For designs eliminating heterogeneity in one direction, 

the necessary and sufficient condition for orthogonality 
of estimable block contrasts and estimable treatment con- 
trasts is found to be that ni;/r, Should be constant for all 


Zand nijj z,; constant for all 7, 7,; being the number of times 


the 7th treatment occurs in the jth block, 7, the replication 
of the ith treatment and k; the number of plots in the jth 
block. Hence it is shown that orthogonality cannot be 
achieved in a connected incomplete block design. 

A similar set of necessary and sufficient conditions for 
orthogonality is also derived for general designs involving 
three groups of parameters, which are used for elimination 
of heterogeneity in two directions. 

A measure of non-orthogonality is defined in terms 
squares of extreme values of correlation coefficients between 
estimable block and treatment contrasts and its value 
computed for Balanced incomplete block and Group divis- 
ible partially balanced incomplete block designs. 


(I. M. Chakravarti) 
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9.7 (6.8) 


Design of a biological assay method for the detection and estimation of estrogenic 
residues in the edible tissues of domestic animals treated with estrogens—In English 
Endocrinology (1958) 63, 806-815 (11 refs., 2 tables, 4 figs.) 


Experiments are described to estimate the residual estrogenic 
activity arising from the administration of stilbestrol to the diet 
of laboratory mice. The logarithm of the ratio of uterine 
weight to final body weight was generally found to be linearly 
related to dose of stilbestrol. 

The treated mice were those fed tissues taken from other 
animals which had been fed diets containing stilbestrol. 
The control mice were fed tissues taken from animals which 
had not been fed stilbestrol. Four experimental designs 
are described to compare the tissues of these two groups of 
mice. 

A. Two-curve method. Stilbestrol was added to both 
control and treated tissue diets at concentrations 
of 0, 2, 4 and 6 parts per billion. Dosage response 
lines were obtained for both groups of mice. The 
amount of residual activity, M, in parts per billion 
of stilbestrol, was estimated from the horizontal 
distance between the two lines. This method could 
not be used if M>2. 


B. Range-finder method. If M>2, control diets were 


used with o, 3 and 6 parts per billion of stilbestrol 
added; five other groups were fed treated tissue 
diets diluted with 0, 50, 75, 87°5 and 93°75 per cent. 
of control tissue diet. ; 


C. Confirmatory dilution assay. Based on the range- 
finder experiment, the treated tissue diet was diluted 
with control tissue diet so that the expected con- 
centration of residual stilbestrol would be between 
1 and 2 parts per billion. A new two-curve experi- 
ment was then conducted. 

D. Single point assays. These are to be used only in 
case there is insufficient tissue to perform a two- 
curve experiment. 

Preliminary experiments had indicated that rats required 

3 to 4 times more tissue than mice in detecting given residual 
estrogenic activity. Statistical procedures are those of Bliss 
(The Statistics of Bioassay, 1952, Academic Press, Inc., 
New York). 


(M. Carter) 


9.1 (7.4) 


1] 


ADAM, A. (Institut fiir Statistik der Universitit Wien) 
Entropy and variance—In German 
Metrika (1958) 1, 99-110 (10 refs. 2 figs.) 


The paper is intended to direct the interest of the statistician 
towards information theory. Some relations are shown 
between entropy, in the sense of Shannon and variance, as 
commonly used in statistics. 

The author treats information theory based on axioms 
as formulated by Chintschin and Faddeyew ; he introduces 
entropy in Chintschin’s definition and then gives Faddeyew’s 
axiom 3’, which concerns the increase of entropy due to the 
addition of one more item (p, is divided into q,.p, and 
Y2Pn, +g = 1). The author points out that an analogue 
relation holds true for the variance. 

According to the author, the relationship between the 
calculus of variance and of entropy, respectively, would 
allow the formulation of several statistical methods, based on 
the calculus of variance, in an analogous way for cases where 
only the concept of entropy is applicable. 

A second example is concerned with regression analysis. 
Again the formalisms in the two “ languages ’”’ are given. 
It is shown that at least in special cases the coefficient of 
determination based on entropy may be better than a measure 
of contingency. 

Finally, it is shown that for a given variance the normal 
distribution has maximal indeterminacy. 


(H. Zemanek) 
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The after-history of pulmonary tuberculosis : 


10.5 (7.0) 


10.9 (10.2) 


a stochastic model—In English 


Biometrics (1958) 14, 527-547 (13 refs., 6 tables, 3 figs.) 


. Stochastic models of chronic illnesses have resulted from 
the realisation that the course of a disease exhibits regularity 
in groups of patients although it is unpredictable in indi- 
viduals. Validity of use depends upon whether clinical 
events are genuinely random in nature. Application to 
chronic disease requires a clearly defined starting point in 
time (diagnosis, completion of treatment, relapse) and 
definition of mutually exclusive clinical states. If instan- 
taneous transition between states is possible, the model 
may be made time continuous. Sometimes in chronic 
disease such instantaneous transition is precluded by the 
choice of clinical states. The model discussed in this paper 
contains two such states, arrested pulmonary tuberculosis 
and active pulmonary tuberculosis. The definition of 
arrested pulmonary tuberculosis requires that certain 
conditions be present during a six-month period before 
entrance into the state is possible; otherwise, the disease 
continues to be classified as active. The conditional pro- 
bability function for ‘‘ patient has arrested disease at time 
t,’ given that ‘“‘ patient had active disease at time ¢, ’, may 
lack a derivative and may be discontinuous. 

If inquiry as to state is made at specified times, a dis- 
crete stochastic process will result. The transition pro- 
bability function will be discrete and questions of differ- 
entiability and continuity will not arise. ‘The process may 
be stationary even if the continuous process was not. ‘The 
present model is such a case. Discrete processes permit 
more freedom of definition and since clinical data are usually 


summarised at discrete intervals of time, no information is 
lost. To avoid too great complexity stationary Markov 
processes are a necessity in any but the simplest situations. 
Data studied have been truncated at a time period before 
death from other than tuberculosis to avoid non-stationarity. 
Accuracy of classification of cause of death is examined to 
ensure validity of truncation. A study is made of some 
follow-up data of pulmonary tuberculosis. ‘The following 
six states were considered in the model (1) active tuberculosis 
that will remain active indefinitely, (2) arrested tuberculosis 
that will remain arrested indefinitely, (3) tuberculosis 
that has proved fatal, (4) arrested tuberculosis that is certain 
sometime to become active, (5) active tuberculosis that is 
certain sometime to become arrested, and (6) active tuber- 
culosis that is certain sometime to be fatal without ever 
becoming arrested. Transition probabilities are assumed 
constant with some “ impossible ”’ transitions having zero 
probabilities. The problem is reduced to one of estimating 
six unknown probabilities. Maximum likelihood estimators 
are derived with suitable iterants suggested. "The method 
is applied to the follow-up data and reasonably good agree- 
ment between observation and prediction is obtained. It is 
suggested that better agreement would result from a finer 
classification of ‘‘ active ’’ tuberculosis. 


(M. E. Turner) 
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BARKER, J. S. F. (Department of Animal Husbandry, Sydney University) 
Simulation of genetic systems by automatic digital computers. Parts III and IV—In English 


10.9 (11.7) 


Aust. F. Biol. Sci. (1958) 11, 603-612 & 613-625 (10 refs. 2 tables, 2 figs., & 22 refs., 3 tables, 3 figs.) 


III. Selection between alleles at an autosomal locus. 


‘This paper describes a Monte Carlo approach to analysis 
of the effect of selection on gene frequencies. An electronic 
digital computer, the SILLIAC at Sydney University, is used 
to simulate the selection processes operating in populations. 

The programme simulates selection between two 
autosomal alleles, allowing for selection at four stages of 
the life cycle: (i) Zygote selection: ability to survive 
from fertilisation to sexual maturity; (ii) Reproductive 
selection: differential reproductive ability of different 
genotypes ; (iii) Selection between gametes at meiosis in 
heterozygotes ; (iv) Selection between gametes on their 
ability to take part in fertilisation. ‘These, and the popu- 
lation size can be varied to investigate different genetic 
situations. 

To test the programme, two experiments reported in the 
literature have been simulated: (i) Selection between ST 
and CH chromosomal arrangements in the third chromosome 
of Drosophila pseudoobscura [T. Dobzhansky and O. Pavy- 
lovsky, Evolution (1953) 7, 198-210]; and (ii) Selection 
between glass and wild type in Drosophila melanogaster 
[D. J. Merrell and J. C. Underhill, Genetics (1956) 41, 
469-485]. 

Close agreement was obtained in (i): the importance 
of an adequate estimate of generation length in these com- 
parisons is discussed. Agreement with the results of (ii) 
was not so close: possible reasons for this in terms of 
selective mating are discussed. 
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IV. Selection between alleles at a sex-linked locus. 


A programme, simulating selection between two alleles 
at a sex-linked locus, and developed for SILLIAC introduces 
selection and chance effects at four stages of the life cycle. 

‘Two experiments on competition between sex-linked 
alleles are simulated, selection between (i) the ‘‘ sex ratio ” 
and standard X-chromosomes of Drosophila pseudoobscura 
[B. Wallace, Evolution (1948) 2, 189-217], and (ii) yellow 
and its wild-type allele in D. melanogaster [D. J. Merrell 
and J. C. Underhill, Genetics (1956) 41, 469-485]. 

(i) Wallace’s data includes estimated overall selective 
values, and selective values for various components of fitness 
from which specific selective values for the programme’s 
four stages are calculated. Two sets of simulated popu- 
lations, one using the overall, the other the specific selective 
values, gave similar results. Though completely unlike the 
experimental population, they were very similar to the 
theoretical curve of change in frequency of sex ratio. 
Difficulties of estimating selective values from studies of the 
components of fitness are discussed, the need for an overall 
measure of fitness being emphasised. 

(11) ‘The simulated populations for competition between 
yellow and wild type agreed closely with the experimental. 
The importance of selective mating in this competition is 
considered. 

' The author concludes that it is possible to simulate 
selection between two alleles both at an autosomal and a 
sex-linked locus by using automatic digital computers. 


(J. Gani) 


10.9 (0.0) 


The existence and stability of selectively balanced polymorphism at a sex-linked locus 
Aust. $. Biol. Sct. (1958) 11, 598-602 (9 refs., 2 tables) 


In a large random mating population, and in the absence of 
selective differences, the genotypes at a completely sex-linked 
locus with two allelomorphs A and a have the equilibrium 
frequencies p2AA, 2pqAa, and g’aa in the homogametic 
sex (female say), and pA and ga in the heterogametic sex 
(male), p and q being the population frequencies of the 
genes A and a respectively. 

When there are selective differences in the genotypes, 
the corresponding equilibrium values of the gene ratios in 
the gametic outputs are given by uy = (1+ta—2Saa)/ 
(1 +ta—2taSa) for the female and taus for the male, where 
ta is the ratio of the selective values of A to a in the male, 
and Sq, Sa are the selective ratios of aa and AA to Aa in 
the female. Hence it is shown that the necessary and 
sufficient condition for the existence and stability of a 
balanced polymorphism at a sex-linked locus is that 1+t4 be 
greater than 2,S,, and 2taSaa. In particular when ta = 1, 
the necessary and sufficient condition for the existence and 
stability of an equilibrium is that the heterozygote should be 
at a selective advantage over both homozygotes; but 
when tg has any value other than unity, it is neither necessary 
nor sufficient for the heterozygote to be at a selective advantage 
with respect to both the homozygotes in order that a stable 
balanced polymorphism may exist. 


These results are considered in relation to Wallace’s 
[Evolution (1948) 2, 189-217] data for populations of Droso- 
phila pseudoobscura in which selection is acting on the sex- 
linked condition ‘‘sex-ratio”’. It is also verified that the 
occurrence of heterozygotes with a frequency greater than 
that of the homozygotes does not necessarily indicate that 
the heterozygote has a selective advantage over both homo- 
zygotes. 


(B. R. Bhat) 
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BLACKWELL, D., BREIMAN, L. & THOMASIAN, A. J. (University of California, Berkeley) 
Proof of Shannon’s transmission theorem for finite-state indecomposable channels—IJn English 


Ann. Math. Statist. (1958) 29, 1209-1220 (6 refs.) 


The information source considered here has a finite number 
D of states, transitions between which are governed by an 
indecomposable D by D Markov matrix M. At the end of 
each transition one of an alphabet of A letters is emitted 
into an information channel. The channel has R states 
between which a transition then takes place, governed by 
one of AR by R Markov matrices C(1), ..., C(A) selected 
by the input letter. A letter is emitted at the output of the 
channel, depending on the final state of this transition. 
The combination of source and channel is then a finite-state 
Markov process governed by a DR by DR matrix L. The 
channel is said to be indecomposable if the matrix L is in- 
decomposable for every possible source. A definition of the 
capacity of such a channel is given. As a criterion for 
indecomposability it is shown that the channel is indecompos- 
able if, and only if, every possible finite product of the channel 
matrices C(1), ..., C(A) is an indecomposable matrix. It is 
shown that every channel with finite memory is indecompos- 
able, so that the results of this paper include those of earlier 
work on finite-memory channels. The authors then prove 
a modification of McMillan’s theorem, that there exists an 
entropy h for any Markoy process governed by an indecom- 
posable matrix, and that for almost all sequences of length 
N and probability p,, ee (—log p,)/N = h. Shannon’s 
co 


theorem is proved in two parts, showing that through an 
indecomposable channel (a) it is possible to transmit at any 
rate less than capacity and (6) it is impossible to transmit 
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16.5 (1.5) 


at any greater rate. Finally, in another form of Shannon’s 
theorem it is proved that an indecomposable channel of 
capacity H is “‘ adequate ”’ to transmit information produced 
by an ergodic source at rate h provided h <H, and that such 
a channel is not adequate if h>H; a precise definition of 
adequacy is given. 


(C. W. Helstrom) 


10.1 (10.2) 


An unfamiliar aspect of the stochastic theory of ‘‘ dead ’’-time counters—In French 


Ann. Soc. Sci. Bruxelles (1958) 72, 81-96 


This paper deals with the theory of ‘ paralysis ’’-counters 
[W. Smith, Proc. Camb. Phil. Soc. (1957) 53, 176] with a 
constant dead time 7 and Poisson input. Denoting by to 
an instant on which a particle has been registered, it is well 
known that the time-axis from ft) onwards is divided into 
contiguous intervals, alternatively ‘‘ covered” (with length 
7), and “ open ” having a length which is a random variable 
with an exponential distribution. The situation one meets 
in going backwards along the time axis from fy) is not sym- 
metrical with the one going along the time axis from fo 
onwards, due to the way in which a counter behaves. 
According to the author it is because of this asymmetrical 
situation not being evident that the division in intervals 
as indicated does hold for the backward case. He then 
proves it to be correct: the proof proceeds by analysing 
the succession of events preceding fp. 

The paper ends with an application of these results to 
the problem of estimating the spectrum of the flight-time 
of particles. 


(H. Breny) 
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FERRERI, C. (School of Statistics, Palermo) 


On the generalisation of Vianelli’s simultaneous conditional means—In I talian 


Statistica (1958) 18, 693-707 (7 refs., 1 table) 


The author develops the statistical problem of conditional 
quasi-arithmetic means defined by Vianelli. A convenient 
application of this model to simultaneous power means is 
discussed in this paper. The opportunity of such an ex- 
tension has also been considered in connection with other 
mean values. 

Cases of this character occur in the analysis of statistical 
time series. For a set. of » observed values in a time series 
Xi (1 = 1, 2,..., 2), each at different points of time t, the 
above-mentioned scheme is based on a minimum problem 
as follows 

od: 
O(X,) = 2 2 (xi4—Xi)? = minimum 
1=1t=1 
under the general condition 
«BX, ..., Xn3 b) = 0 
where the quasi-arithmetic means X;, in the Vianelli sense, 
are unknown and 6 is known. 

In order to determine X; by the well-known Lagrange 

method, it is necessary to differentiate the linear relation 
D(X) +AF(X; 5 b) 
with respect to each of the n+1 unknowns Xj, ..., Xp, A. 

By analogy with the above-mentioned problem, the 
author suggests a general expression from which the quasi- 
arithmetic means follow as special case. This generalised 
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scheme, based on conditional simultaneous power means, 
includes quasi-harmonic, quasi-geometric and quasi-quad- 
ratic means. 

At the end of the paper the author gives some numerical 
examples principally concerning the Cobb-Douglas function 
In connection with the theory of marginal productivity and 
for which he uses the original statistical data. 


(V. Amato) 


10.3 (10.1) 


Asymptotic distributions of “* psi-squared ”’ goodness of fit criteria for 


mth order Markov Chains—IJn English 
Ann. Math. Statist. (1958) 29, 1123-1133 (16 refs.) 


A probabilistic model for the random behaviour of an object 
which moves at discrete time points from one to another of 
a finite number s of states is said to be a Markov chain of 
order k (an integer 2o) if the conditional (transition) 
probabilities of the object entering a certain state at time f, 
given the whereabouts of the object at the previous k instants 
of time, suffice (together with certain initial probabilities) 
to determine all other probabilities in the model. The case 
k = 1 is the usual Markov chain. If k = 0, the behaviour 
of the object at any time ¢ is independent of its past. 
Considerable effort has been spent during the past five 
years on problems of inference in Markov chains. ‘These 
inference problems are of two types: (i) the problem of 
estimating the transition probabilities, and, (ii) the problem 
of testing the hypothesis that a given set of data has come 
from a Markov chain of a specified order. The paper being 
reviewed is one of a series of studies by the author into such 
problems. The present paper, however, is more of a con- 
solidating digression from, rather than a continuation of the 
main trend of his research, since it is concerned with deter- 
mining the validity of several conjectures made by I. J. Good. 
In order to make inferences, the experimenter will 
observe the behaviour of the process in question, the data 
being therefore a sequence of integer random variables, 
(X,, Xo, ..., X7) say, where the value of X; denotes the 
state the object is in at time j. .Let u = (u1, Uz, ..., Um) be a 
set of integers and let f,, be the frequency of times the object 
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passes through the states u,, uo, ..., Um consecutively during 
the first T units of time. ‘These statistics are sufficient for 
most inference problems. Let f,,,, denote the expected 
value of f,, under the assumption that the process considered 
is a particular Markoy chain of order n, and set f m n equal 


to the expectation of 7, under the assumption that the 
process is the Markov chain of order 7 having for its transi- 
tion probabilities their maximum-likelihood estimates com- 
puted from the observed data. For testing the hypothesis 
that a process is a Markov chain of order n on the basis 
of the f,, the following chi-square-like goodness-of-fit 
criteria have been proposed ; 


2G tifa, and $2.2 (hf, 


For these tests to be practical, their limiting distribution 
(T-+>-+0) should be known. Goodman proves the follow- 
ing conjecture of Good. Under the composite hypothesis 
that the process is a Markov chain of order n, the limiting 
m—n—l 
distribution of pees is that of 22 Yg; where Yg; is a chi- 
; t=1 
square random variable with g; = (i—1)%7™"*71 degrees of 
freedom, all random variables being independent. A second 
conjecture by Good concerning the asymptotic distribution 
of pe n 38 not proved but a modified version is stated 


and proved. 
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GUPTA, H. C. (University of Delhi) 
Diffusion by discrete movements—In English 
Sankhya (1958) 20, 295-308 (5 refs., 2 figs.) 


The problem of diffusion considered in this paper is that of 
the motion of a particle along a straight line with one of two 
possible velocities u,(>0) and u,(<o). It is governed by a 
specified initial distribution and a matrix of transitional 
probabilities. This can then be considered as a homo- 
geneous Markov chain with two possible states. Let »,(v.) 
be the time spent during the first 7 intervals in the state u,(u,) 
when it ends up in the state 1,(u.). The first section of the 
paper is devoted to finding (i) the characteristic function, 
(ii) the probability distribution, (iii) the moments of the 
random variables », and v,. In the limiting case, when the 
discrete movements become continuous, the limiting values 
of the mean and mean-square are evaluated. The modi- 
fications necessary, in case immigration or emigration of 
particles is present, is touched upon. 

In the next section differential equations governing 
continuous diffusion are obtained. It is a second-order 
partial differential equation which contains as particular 
cases, (i) the equation of heat conduction, (ii) the equation 
of motion of a sound wave, (iii) the equation of diffusion 
by continuous movements with a constant rate of emigration. 
The last section is devoted to the investigation of the pro- 
bability distribution of v,(v,) when one or more barriers 
(reflecting, partly reflecting and absorbing) are present. 
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Two queues in parallel—In English 
Biometrika (1958) 45, 401-410 (1 ref.) 


Each new arrival to a system consisting of two queues A and 
B is assumed to join the shorter queue, or, if they are of the 
same length, one particular queue, say queue A. The 
arrivals are assumed to comprise a homogeneous Poisson 
process, and the service is also assumed to be Poissonnian 
(though not necessarily with the same parameters for the 
two queues). Differential-difference equations are set up 
and, assuming equilibrium, various relations are found 
between the state probabilities, the individual (marginal) 
queue-length distributions together with some of the lower 
moments of the conditional queue-length distributions. If 
the length of queue B is independent of time, the problem 
reduces to one previously studied by the same author in 
“ Queueing with balking ” [Biometrika (1957) 44, 360-369]. 

Considerable simplification occurs when it is supposed 
that any member of a queue will instantaneously change 
queues whenever it appears advantageous to do so, 1.e., 
when the other queue is the shorter. Thus the lengths of 
the queues can never differ by more than one unit. Explicit 
expressions for the equilibrium state probabilities are given ; 
each of the three possible configurations (equal queues, 
queue A longer, queue B longer) has (conditionally) a 
geometric distribution. 


(C. L. Mallows) 
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The method of images is employed to reduce it to the case 
when barriers are absent. Some of the results of this section 
are generalisations of certain results of Chandrasekhar. 


(R. Ranga Rao) 
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10.5 (1.1) 


A property of ¥-divergences of marginal probability distributions—In English 


Czech. Math. 7. (1958) 8, 460-463 (5 refs.) 


f-divergence is an entropy-like measure introduced by 
H. Jeffreys to measure the difference of two probability 
distributions P and Q. 7 is the sum of relative entropies of 
P with respect to Q and of Q with respect to P. The paper 
is based on a rather general definition of j-divergence, 
holding for probability measures on arbitrary fields of sets. 
The marginal distribution of a p-measure means the measure 
restricted to a subfield of the original domain of definition. 
So, for the marginal distribution, the ¥-divergence is also 
well defined. 

At first, certain properties of $-divergence are listed. 
It is a symmetrical function of P and Q, it is always non- 
negative and equals zero if, and only if, both P and Q are 
identical. It is shown by a counter-example that #-diver- 
gence cannot be used as a distance in spaces of measures 
because it does not possess the triangle property. Another 
property is a kind of monotony with respect to the Borel- 
field. ‘That is, the #-divergence of marginal distributions 
is always less or equal to the divergence of the original field. 
The equality holds if the marginal 7-divergence is infinite or, 
when the subfield is sufficient to distinguish between P and 
Q. Let us note that the divergence can be finite only if 
P and Q are equivalent measures. 

The property announced in the title is the following: 
let us have an increasing sequence of Borel-fields and let 
F be the smallest Borel-field, containing every field of the 
sequence, then the #-divergence of measures P and Q 
with respect to F is the limit of #-divergences on the fields 
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of the sequence. In the proof, the case of an infinite limit 
is first excluded as trivial ; then it is shown that the measures 
P and Q are equivalent on F if the limit is finite. For the 
last step Doob’s theorem on martingales is needed. By an 
analogous method similar results concerning continuity 
of relative entropies can be obtained. Also if the limit of 
relative entropies of P with respect to Q is finite, probability 
P is absolutely continuous with respect to O on the field F. 

The theorem, when applied to stochastic processes, gives 
the result that the }-divergence between two distributions 
of a random function is the least upper bound of the diver- 
gences of finite-dimensional distributions. In this case the 
field F is the smallest Borel-field, generated by finite- 
dimensional cylinder-sets. It is shown that divergence of F 
equals that on a Borel-subfield generated by the values of 
the random function on a countable set of time-instants 
and the assertion easily fellows. In an article in the next 
number of the same journal the author gives an application 
of this theory to the problem of deciding, with probability 
zero of error, between two Gaussian distributions of a 
random function, when one realisation of the function has 
been observed. 


(P. Mandl) 


10.0 (—.—) 


A property of normal distributions of any stochastic process—In Russian 


Czech. Math. #. (1958) 8, 610-618 (6 refs.) 


In a previous paper in the same journal (Czech. Math. 7., 
8, 460-463) the author proved that 7-divergence (H. Jeffreys) 
of two probability distributions P and Q of a stochastic 
process is equal to the least upper bound of 7-divergences 
of finite-dimensional distributions of that process. Now, 
he uses this property to show that two normal distributions 
of a stochastic process can be only equivalent or singular. 
The distribution of a stochastic process is called normal 
when its finite-dimensional distributions are multivariate 
normal. Equivalency of P and Q means that a set has 
probability zero under P if, and only if, it is of Q-measure 
zero. P and Q are called singular if there exists a set which 
has P-measure one and Q-measure zero. In the case of 
singular distributions we can decide, with probability zero 
of error, between P and Q from one realisation of the process ; 
so there are no statistical problems. The author’s result 
can be restated in the form that, in any set of normal dis- 
tributions of a stochastic process from which no two are 
singular, we have probability densities (Radon-Nikodym 
derivatives) of all measures with respect to any member of 
the class. The probability densities are of interest, for 


instance, when we want to form most powerful critical 
regions. 
The proof is based on the following property of 7-diver- 
gence: to an arbitrary positive number e there exists a 
constant K such that we can decide from one observation, 
with probability of error less than e, between two arbitrary 
regular n-dimensional normal distributions of a random 


vector whose divergence is greater than K. When the 
divergence of two normal distributions is finite, they are 
equivalent. When it is infinite, because it is the least upper 
bound of finite-dimensional divergences and because of the 
property mentioned, a sequence of critical regions can be 
constructed, with probability of incorrect decision (both 
kinds of error) tending quickly enough to zero. From that 
sequence the region with probability zero of error is con- 
structed. 

The criterion of equivalence (finiteness of 7-divergence) 
is illustrated in two examples : 


(1) Change of the scale parameter. In this case P and Q 
have the same mean-value function, but the covariance 
function of Q is that of P, multiplied by a positive constant 
C. The result is that when C is different from one, the 
measures are equivalent if, and only if, they degenerate in 
finite-dimensional distributions. 

(2) Translation of a Gauss-Markoff process. In this case 
a necessary and sufficient condition of equivalency of two 
normal distributions, which differ only in their mean-value 
functions, is the finiteness of a Hellinger-integral. 


(P. Mandl) 
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HANNAN, E. J. (Australian National University, Canberra) 
The estimates of the spectral density after trend removal—In English 


10.6 (10.2) 


F. R. Statist. Soc. B (1958) 20, 323-333 (7 refs., 2 tables) 


Consider a discrete time-series which can be represented 
as the sum of (a) a multilinear regression on several known 
variables, with unknown coefficients, and (b) an independent 
stationary random process with an unknown spectrum. 
It is supposed that the regressor variables are such as to 
allow the (asymptotic) definition of autocovariance functions. 
Thus polynomial and trigonometric regressions are included. 
A further “ orthogonality” condition is imposed on the 
regressor variables. 

The problem is then to estimate the unknown spectrum 
when the regression coefficients are estimated by least 
squares—which is asymptotically equivalent to using the 
best linear unbiased estimates. This spectrum is needed to 
obtain the covariance matrix of the regression estimates. 

The estimation of the regression coefficients introduces a 
first-order bias into the straightforward (unsmoothed) 
estimate of the spectrum from the residuals. This bias is 
evaluated and is shown to be a simple function of the spectra 
of the regression variables. It corresponds (asymptotically) 
to multiplying the straightforward estimate by a factor which 
can be calculated a priori. ‘The sample autocovariances 
are also biased and these biases are similarly evaluated. 

To obtain a consistent estimate of the spectrum it is 
necessary to smooth the straightforward estimate. This 
introduces a further bias which will be relatively negligible 
if the population autocorrelations decrease sufficiently 
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The estimation of distributed lags—IJn English 
Econometrica (1958) 26, 553-565 (9 refs.) 


The purpose of this paper is to advance the understanding 


ice) 
of the lag scheme y, = « 2 A‘’x;_i:+u,; proposed by L. M. 
i=0 
Koyck [Distributed Lags and Investment Analysis (1954)]. In 
particular, the problem of estimating the parameters « and 
dis considered. A derivative of the above equation is used 
which involves only three observable variables, yz, yz and 
x, Ordinary least squares procedures, however, leads to 
biased estimates of the coefficient parameters since the com- 
posite error term in the derived form of the model is not 
independent of y;_, and also because the composite error 
terms are serially correlated. 

Koyck’s method of estimation, when the original series 
of error terms are not autocorrelated, is shown to be equival- 
ent to a type of maximum likelihood estimate. An alter- 
native computing method is suggested which enables the 
estimates of « and A to be obtained in one step and not two. 
In addition a generalised statement of the problem, with 
zero autocorrelation in the original error terms, is given which 
enables the method of solution to be generalised to any 
number of variables. 

The case of autocorrelation in the original series of error 
terms is then considered. When the autocorrelation para- 
meter is equal to A, it is shown that the least squares regression 
of y, on x; and y;_, would give consistent and efficient 
estimates. In other cases it is necessary to transform the 


rapidly. Both biases will be asymptotically negligibly 
relative to the standard error of the estimate. However, 
since the samples generally encountered in practice will be 
finite, it is thought worthwhile to investigate methods of 
eliminating them. Several suggestions are made, the 
simplest corresponds to multiplying the smoothed estimate 
by a factor (a smoothed version of that mentioned above) 
which again can be calculated a priori. This will not affect 
the order of the asymptotic variance of the estimate, but will 
reduce the order of the overall bias to that of the bias due 
to smoothing. 

A numerical example is given; it is admitted to be 
favourable to the use of the proposed factor since the bias 
due to smoothing is relatively small. 


(C. L. Mallows) 


Foundation, Yale) 10.2 (4.2) 


variables in order to obtain equations not involving auto- 
correlated errors. An equation system is developed which 
can be solved for «, X and the autoregressive parameter by 
an iterative procedure. 

The paper concludes with a general discussion of the 
suitability of the general model and its derivative for the 
investigation of particular problems in economics. It is sug- 
gested that with cross-section, as opposed to time-series 
data, the techniques for zero autocorrelation would be 
appropriate. 


(D. J. Day) 
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KOZNIEWSKA, I. (Planning & Statistics High School, Warsaw) 


1O.r (1.5) 


Ergodicity of non-homogeneous Markoff chains with two states—In Polish 


Collog. Math. (1958) 5, 208-215 (6 refs.) 


The scope of the note is to compare different notions of 
ergodicity for non-homogeneous Markoff chains with two 
states. The notions of ergodicity of a Markoff chain with r 
states, considered by Mrs Kozniewska, can be formulated 
in terms of transition probabilities h,, (m, n) from a state 
E, at moment m to a state H, at moment 7 as follows : 


I. Ergodicity in Kolmogoroff’s sense: for every 
m = 0, 1, ... and any fixed a, b and c the difference between 
Nig» (m, n) and hy (m, n) tends to zero if n increases to infinity. 

II. Weak ergodicity in Hajnal’s sense: for any a, b, ¢ 
the difference between hyp (0, n) and h,, (0, n) tends to zero 
with “ increasing to infinity. 


III. Strong ergodicity : for any a, h,, (0, n) tends with n 
increasing to infinity to a limit depending only on b. 

IV. Strongest ergodicity: for every m =o, I, ... and 
any @ the probabilities ha, (m, n) tend with increasing n to a 
limit depending only on 0b. 

Obviously strongest ergodicity implies the strong one 
and Kolmogoroff’s ergodicity implies that of Hajnal. Mrs 
KoZniewska gives examples of two-state non-homogeneous 
Markoff chains which are ergodic (a) in all senses, (6) only 
in the sense of II and III, (c) only in the sense of II, (d) 
only in the sense of I and II, (e) in no sense. Then, using 
the necessary and sufficient conditions for a two-state 
Markoff chain to be ergodic in the four senses in question, 
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formulated in terms of characteristic roots of matrices of 
transition probabilities from a state at moment 7 to a state 
at moment +1, she gives conditions (f) for the ergodicity 
I to be equivalent to ergodicity II, and (g) for ergodicity III 
to be equivalent to ergodicity IV. Conditions are also given 
for a two-state non-homogeneous Markoff chain to be 
ergodic in the strong sense without being ergodic in the 
strongest sense; to be ergodic in Kolmogoroft’s sense with- 
out being ergodic in the strong sense; to be ergodic in the 
weak sense without being ergodic in any other sense. Finally, 
the equivalence of the four kinds of ergodicity for two-state 
homogeneous Markoff chains is proved. 


(S. Zubrzycki) 


10.7 (0.4) 


Note on the characteristic function of a serial-correlation distribution—In English 


Biometrika (1958) 45, 559-561 (6 refs.) 


The moments of a two-parameter distribution appearing in 
several studies of serial correlation have recently been dis- 
cussed both by M. G. Kendall and by J. H. White. 

The present author finds a Neumann-type series for the 
characteristic function of this distribution for a generalised 
parameter range. The derivation of the characteristic func- 
tion is given and an expression is obtained for the pth moment 
about zero of the distribution. Agreement with Kendall 
and White is noted. The consequence of the standardised 
distribution to the normal, as one of the parameters tends to 
infinity, is demonstrated through the characteristic function. 

The author remarks that ‘‘ the technique of asymptotic 
differential equations employed may be useful in other such 
problems. It is familiar in physics under the name W.K.B. 
method.” 


(F. N. David) 
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LESLIE, P. H. & GOWER, J. G. (Bur. of Animal Pop. Oxford & Rothamstead Exp. Station) 
The properties of a stochastic model for two competing species—In English 
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Biometrika (1958) 45, 316-330 (11 refs., 5 tables, 4 figs.) 


Conditions are found for the stationary values of the system 
in terms of the parameters of a discreet time, interacting, 
deterministic model. Inequalities which arise lead to three 
types of stationary state, (i) stable, in which both species 
_ persist, (ii) unstable, in which either one or the other species 
persists according to the initial state of the system and (iii) 
in which one species persists and the other disappears 
independently of the initial state. 

The variety of stochastic models having the same deter- 
ministic equivalent may be considered to lie between two 


extremes : Model I in which the birth rate of each species | 


is constant, and Model II in which the death rate of each 
species is constant. A stochastic process is constructed by 
attaching random normal deviates to the number of each 
species in existence at each stage of the deterministic 
development. As a consequence of this the outcome of a 
system with an unstable stationary state (ii) depends not 
only on the initial state but on the early development of 
the system. Using sixty replicates on each of a range of 
initial states, contour lines for the probability that a given 
species will survive (are estimated) for both Models I and II. 

In the case of the stable stationary state (i) an analysis 
is carried out on the life histories of the two systems with 
the same initial states and deterministic model in which the 
birth rate of each species is constant (Model I). It is found 
that the fluctuation of the number of each species in existence 
about the stationary state follows a bivariate distribution 
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LEVY, P. (Ecdle Polytechnique, Paris) 
On various classes of random functions—In French 


$. Math. Pures Appl. (1959) 38, 1-23 (8 refs.) 


The author gives in this paper a complement to a theory of 
Laplacean random functions presented by him at the Third 
Berkeley Symposium and the more general study of linear 
correlation functions. The paper begins with work on 
random functions in time 4(t) which are, for all fixed (2), 
certain linear functions of values in (0, t) of an additive 
random function. 

The author represents these in the two forms of Stieltjes 
integrals ; then he shows that they are identical. ‘The con- 
fines of this set of functions contain all the linear correlation 
functions. 

The proofs use theorems of Frechet and Damois general- 
ising the results of S. Bernstein. The paper contains eight 
theorems and the author points out the interest of this ques- 
tion for the problem of classification of random functions. 


(D. Dugué) 


which is approximately normal but that the variation about 
the stationary state is greater than that which may be 
predicted from theoretical considerations. 

By using sixty replicates for each of a variety of initial 
states, an estimate is made of the probability that the 
stochastic outcome is the same as the outcome of a deter- 
ministic model with a type (iii) stationary state using Model I. 

Life histories are obtained by reading into a computor 
the parameters of a deterministic model together with the 
initial conditions and a pair of random numbers. These 
generate the pseudo-random deviates required at each stage 
of the stochastic development. 


(J. G. Saw) 
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MORAN, P. A. P. & WATTERSON, G. A. (Australian National University, Canberra) 
The genetic effects of family structure in natural populations—In English 


Aust. 7. Biol. Sci. (1959) 12, 1-15 (7 refs.) 


This paper deals with the rate of approach to homozygosity 
of dioecious diploid populations whose generations do not 
overlap. The probability that a finite genetic population 
without mutations is not totally homozygous at the t-th 
generation (f= 0, 1, 2, ...) is asymptotically equal to 
Cn, where A is called the rate of approach to homozygosity, 
and C is a constant. 

For certain models of Wright [Genetics (1931) 16, 97- 
159; Proc. Nat. Acad. Sci., Wash. (1933) 19, 411-419], 
Moran [Proc. Camb. Phil. Soc. (1958) 54, 60-71], and 
Watterson (unpublished data) with specific distributions of 
offspring per individual of the parent generation, the known 
values of A are given. ‘These results are generalised to cases 
with arbitrary offspring distributions, and with permanent 
marriage. 

First, the family structure of populations of fixed size 
N is considered, with N, males and N, females where 
N,+N,.=N. The authors obtain the generating functions 
P,(z), P.(z) for offspring distributions of male and female 
parents respectively, in terms of the birth and death rates 
for each individual of specified age ; the first two moments of 
these distributions prove to be of importance in subsequent 
theory. As examples, the generating functions applying 
to the previous population models are given. 

The effect of family structure on the genetic behaviour 
of a population after a large number of generations is now 
studied ; the population is assumed to be diploid, with a 
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diallelic genetic factor. "The autosomal case in a non- 
overlapping generation model with independent random 
mating between adults in each generation is considered first. 
It is shown how A, the largest non-unit root of the stochastic 
matrix for the Markoff chain of genetic states, may be 
evaluated from a set of difference equations as approxi- 
mately 
A= 1—N*NyP," (1) +N2P2" (1)]/8. 


This value remains unaffected when permanent marriage 
occurs, with N, = N, = N/2 and identical variances for 
both offspring distributions. Further results are stated for 
the sex-linked case. 

Finally, though the general solution for overlapping 
generation models is not known, it is shown that A is the 
same for a pair of models of Moran’s type whether genera- 
tions overlap or not. That'this result does not hold in general 
has been shown by Watterson, using a second pair of equiva- 
lent models. The conditions required to extend results in 
non-overlapping to overlapping generation models require 
further investigation. 


(J. Gani) 


10.6 (10.2) 


On asymptotically efficient consistent estimates of the spectral density 


function of a stationary time series—In English 


F. R. Statist. Soc. B (1958) 20, 303-322 (16 refs., 2 tables) 


The author considers continuous and discrete time series 
which are. stationary (in the wide sense) to order four, and 
constructs several general classes of estimates of their 
spectral density functions. Each of these estimates can be 
regarded as a weighted average over the periodogram, or as a 
weighted sum of the sample autocorrelations. ‘There are 
two important types of estimate, called exponential and 
algebraic, corresponding to weight functions which satisfy 
certain fairly weak inequalities. 

Considering the case where the period of observation 
tends to infinity, definitions are given according to which 
an estimate may be pointwise consistent or boundedly con- 
sistent; uniformly consistent; uniformly-boundedly _con- 
sistent ; functionally-uniformly consistent ; and integratedly 
consistent—depending on the integrated mean square error 
(i.m.s.e.). Each of these estimates being of a certain order, 
it is possible to define the minimum integrated mean square 
error and the asymptotic efficiency of an estimate. 

The highest order of consistency with which the spectrum 
can be estimated depends on the smoothness of the spectrum, 
as expressed by the rapidity of the decrease of the popu- 
lation autocorrelations as the lag tends to infinity. ‘The 
results are very different according to whether the decrease 
is (a) exponential, or (6) algebraic. The asymptotic 
im.s.e.’s for estimates of the exponential type [for (a)] 


and of the algebraic type [for (6)] can be evaluated and give 
simple results. In the first case, the actual choice of a weight 
function is not critical, but in the second case it is very 
critical. The behaviour of the autocorrelation function 
need not be known at all completely in the first case, but in 
the second case it must be known fairly well before it is 
possible to construct an efficient estimate. The class of 
spectra which are consistently estimated (in an integral 
sense) by a given estimate can in each case be simply 
characterised. 

Similar results are obtained for asymptotic mean square 
errors corresponding to the other definitions of consistency. 
The minimum i.m.s.e.’s are evaluated numerically for each 
of eight suggested estimates in the case of a simple Markoff 
process and for three sample sizes. 


(C. L. Mallows) 
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PFANZAGL, J. 
On the decomposition of statistical series—In German 


Metrtka (1958) 1, 130-147 (11 refs., 2 tables, 7 figs.) 


The paper deals with the analysis of statistical time series 
in components, starting from a type of series which can be 
represented by an additive composition of a smooth com- 
ponent (trend and business-cycles), a periodical component 
(season) and an irregular component. 

The author stresses the fact that the methods for the 
analysis of time series are frequently based on the assump- 
tion that the values of the seasonal component at correspond- 
ing points (e.g. in January, February, etc.) of different 
periods (e.g. years) are either constant or change only in a 
special way, depending on the corresponding values of the 
smooth component. In the latter case the amplitude of the 
seasonal component is usually assumed to be proportional 
to the smooth component. There exist, however, series 
which do not comply with these assumptions. The method 
of A. Wald takes this fact into account by starting with a 
more general model: the periodical component is assumed 
to have a constant shape, its amplitude is assumed to change 
only gradually, but in an arbitrary way, i.e. independent 
of the development of the smooth component. The author 
varies this model in such a way that the amplitude may be 
subjected to sudden changes from one period to the other. 
In accordance with this, the season is defined as the product 
of a periodical function, and—differently from Wald—of a 
magnitude (amplitude) which can assume a different arbitrary 
value every year. Consequently, the values of the seasonal 
component have to be computed separately for every period. 
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Under the usual assumptions on the behaviour of the 
irregular component this is possible if the smooth com- 
ponent can be approximated in every period by a linear 
function of time. By subtracting the seasonal component 
from the. original series a composite series of the smooth 
plus the irregular component is obtained. These two com- 
ponents can be separated with the variate difference method. 
The accuracy of the method is examined on the basis of an 
artificial series and its applicability demonstrated by analys- 
ing the series of scarlet-fever deaths in the years 1949 to 1957. 

The procedure developed by the author is a generalisa- 
tion of Wald’s method, as the assumption of an arbitrary 
change of the seasonal amplitude is less restrictive than the 
assumption that there is a continuous change in time. If the 
assumptions underlying Wald’s model are fulfilled both 
methods lead to practically the same results. 

The new method can be used especially for the analysis 
of series showing an unexpected change of seasonal fluctua- 
tions. This will mostly be the case with agricultural 
production, infectious diseases, etc. Finally, attention is 
drawn to the possibilities of applying such methods 
to the analysis of agricultural and industrial experiments, 
as in this field formally equivalent relations occur. 

A correction of several misprints is subsequently given 
in this volume of Metrika (1958). 


(O. Anderson, jun.) 


10.8 (1.4) 


On secondary processes generated by a random point distribution of Poisson type—In Hungarian 
Ann. Univ. Sci. Budapest—Sect. Math. (1958) 153-170 (13 refs.) 


In this paper an exact mathematical model of secondary 
stochastic processes and a condition under which the 
secondary process will be of “‘ Poisson’s type ” is given. 

Let T be an abstract space and Yr a o-algebra of 
certain subsets of 7. Suppose that in T a random point 
distribution of Poisson type is given, i.e. a random selection 
of a finite number of points such that if é(A) denotes the 
number of random points lying in Ae Pr then 


PLEA) = A] = (2/R YP), 


where (A) is a finite measure on the elements of Sr. To 
every random point in T there corresponds a secondary 
“‘ happening ”’ i.e. a random selection of exactly one element 
of another abstract space Y in which a o-algebra Sy is 
given. If the random points dre ty, f,, ..., fy and the corres- 
ponding ‘‘ happenings ”’ are 1, Vo, ---) Yn respectively, then 
a sample element of the whole secondary process ISIE siVa)5 ee 
(tn, Yn)» ‘Thus the secondary process is not else than a 


random point distribution in the product space Tx Y. 
The author proves that under a suitable condition on the 
space T and under the condition of independence of second- 
ary happenings corresponding to different random points, 
the random point distribution in TX Y is also of Poisson 


type. If 1(C, 2) is the distribution of the random happening 
y having the starting point teT and AeA then the - 
parameter of the probability distribution of random points 
lying in AX C is 


| u(A, t)\dt). 

Cc 
The parameter corresponding to other sets of ASrX Sy can 
be obtained by extension of the preceding measure. Appli- 
cations are mentioned, e.g. to a telephone problem, to the 
transparentness of a film darked by emulsion spheres and 
to wandering random points in a Euclidean space. The 
latter contains as a special case the result of Doob. 
(Stochastic processes, New York-London, 1953, pp. 404-407.) 


(P. Révész) 


RAMAKRISHNA, B. S. (Indian Institute of Science, Bangalore) 
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Information Theory and some of its applications—Jn English 


Current Science (1958) 27, 376-385 (x1 refs., 1 fig.) 


This is an expository paper dealing with Information Theory, 
as was developed by Claude Shannon, and some of its appli- 
cations to the transmission of messages over channels with 
and without noise together with problems of language and 
human communication. The model of a general communi- 
cation system is as that adopted by Shannon, viz. “ the 
information source”? from which message passes to the 
transmitter which releases the transmitted signal with 
“noise ”’ added, and the received signal passes through the 
receiver which gives rise to message again. This message 
comes to the “destination”’. Various properties of the 
quantity of information—2p, log, p,, where p,, p.... are 
the probabilities of ” symbols, are enunciated to justify 
this particular definition of information. The capacity of 
a noisy channel is defined in the usual manner by maximis- 
ing “‘ the amount of information sent less the uncertainty 
that remains of the transmitted message after the message 
is received at the destination”. ‘To use the channel to 
transmit information at any rate near to its capacity, one 
should assign probability p(;) to the input symbols such 
that 2pG,;, log pu,n/Pu- Du is equal to the required 
rate where pi,;) = p(j | 2). Pa, Pw =P, and 
p(t|j) is the channel probability: i denoting the input 
symbol and j the output symbol. 

Another important branch to which the theory of infor- 
mation can be applied is that of language and human com- 
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SMITH, W. L. 
Renewal theory and its applications—In English 
F. R. Statist. Soc. B (1958) 20, 243-284 (98 refs.) 


This is a long review paper with nearly one hundred refer- 
ences: it is followed by 18 pages of discussion for which 
a separate abstract is given. The topic reviewed is that of 
“renewal theory ”’ (alias the theory of recurrent events). 
A succession of ‘‘ lifetimes” X,, ... X, ... 1S envisaged, 
each an independent sample from the same law, and interest 
is focused on the distribution of the number N of such 
events occurring in a time t (measured from the start of the 
first lifetime). ‘The review is divided into four chapters. 
The first treats the fundamental theory based on the 
“renewal function”? H(t) = &(N;,) and here the emphasis 
is on rigorous proof and generality (which is found to exist 
chiefly in the limit as t>0). ‘The second chapter discusses 
“ regenerative stochastic processes” on the basis of Smith’s 


‘result that the embedded renewal process provides a proof 


of the existence of an equilibrium limiting distribution for 
the regenerative process. This is also done for semi-Markoff 
processes. The chapter closes with a brief treatment of 
cumulative processes (partial sums of lifetimes), their 
normal and degenerate limits. 

The third chapter treats four recent problems of a more 
particular and compound nature. The first is one of Takacs’ 
(arising from Geiger counter theory) where the expectation 
of an infinite sum of a function of two independent renewal 
processes is required together with its auto correlation 
function. Asymptotic results are given together with a 


(North Carolina State College, late Statist. 


munication. Information theory has been developed without 
attention to the meaning of the messages. But when human 
communication is considered it is remarked that the semantic 
problems involved cannot be ignored. But nevertheless 
the problem of semantic information is not entirely separable 
from the “ information of the telecommunication problem ”’, 
sometimes called the ‘selective information’’. For if we 
have the choice between two languages in which to transmit 
the same semantic information, it is suggested that language 
which requires the smaller number of bits of selective 
information should be preferred. By adopting this method 
the author states that he has succeeded in studying the 
relative efficiencies of different Indian languages for com- 
munication of the same semantic content. A note on the 
relationship between Informational and Theromodynamical 
entropies is provided at the end. 


i 


(K. R. Parthasarathy) 
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discussion of the corresponding spectral density function. 
The second problem is of the expectation P(¢) of an infinite 
product related to Takacs’ sum (which ‘has a death and 
immigration process application) and an integral equation 
is found for this. A differential equation is found for the 
corresponding Laplace Transform and its numerical solution 
is discussed together with explicit solution for the case of an 
underlying Poisson process. A series expansion for the 
asymptotic value of P(t) is found in another special case. 
The third problem is that of the application of “ paralysis ”’ 
and “‘ guarantee’ censors to a renewal process and a critical 
catalogue of some thirty relevant papers is made. A variety 
of more or less complete solutions is found for various special 
underlying distributions. "The fourth problem is Harris’s 
treatment of branching renewal processes where, at the end 
of each lifetime, two new lives start and an integral equation 
for the mean population size at time ¢ together with an 
asymptotic expression is given for this. The fourth chapter 
covers : Chang’s extension of renewal theory to vector X's ; 
Cox and Smith’s treatment of the case where the {X;} have 
different p.d.f.s. ; Wold and Doub’s work on the case where 
the p.d.f. of X,; depends, in a Markovian fashion, on the 
previous values of X. 


(D. E. Barton) 
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SMITH, W. L. (Report of a discussion on paper by) 
Renewal Theory and its Ramifications—In English 
J. R. Statist. Soc. B (1958) 20, 284-302 


The discussion on this review paper consists of nine verbal 
contributions, two written contributions received after the 
meeting and the author’s detailed reply. Those who took 
part, in order of speaking, were :—M. S. Bartlett, D. R. 
Cox, J. M. Hammersley, J. G. Skellam, L. R. Shenton, 
A. Mercer, R. Syski and D. G. Kendall. he written con- 
tributions were from L. Takdcs and S. Karlin. 

Professor Bartlett brought forward a “ portmanteau ”’ 
formula concerning renewal theory and also referred to the 
forthcoming publication of work relevant to m-dependant 
renewal process. Dr Cox discussed three generalisations 
including multidimensional renewal theory and a practical 
application involving renewal theory and the theory of 
order statistics. Mr Hammersley’s contribution was 
extensive and dealt in some detail with the method of 
treating an integral equation (section 3.9 of the paper) and 
of the approach of renewal and associated processes through 
integral equations for cumulative distribution functions. Mr 
Skellam discussed a class of immigration-death processes 
where the arrivals are at random. Dr Shenton commented 
upon a renewal process which arose in polymer chemistry 
and Mr Mercer a practical application to conveyor belting 
damage. In this latter case explicit formulae are not 
available and attention must be turned to non-parametric 
results. 

Mr Syski pointed out that the concept of regeneration 
points was introduced by Palm in 1943 [Ericsson Technics 
44, 1] and gave certain additional references to work in 
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connection with equations for infinite sums. (See section 
3.1 of the paper.) Mr D. G. Kendall’s contribution noted 
that there was a good deal of common ground between 
renewal theory and the theory of Markov chains. He had 
carried out work on the rate-of-convergence problem 
associated with a limit theorem due to Kolmogorov and was 
interested to learn whether anything had been done in this 
connection with the analogous theorem of Blackwell. He 
also pointed out that Palm’s term for what is now known as 
“ regeneration point ” was “‘ point of equilibrium ”’. 

Dr Takacs and Professor Karlin provided written con- 
tributions to the discussion. Dr Takacs showed how 
formulae (1.1) of the paper could be obtained without 
recourse to the Strong Law of large numbers and emphasised 
the key importance of the theorem (1.3) in the theory of 
stochastic processes. He gave an application of renewal 
theory to telephone traffic and a generalisation of renewal 
theory in connection with industrial processes. 

Professor Karlin’s written contribution drew attention 
to a further source of stochastic models whose analysis must 
rely upon renewal theory—inventory and_ production 
models. 

The author’s detailed reply summarised the discussion 
and offered comments further upon points raised by the 
participants. 


(Wm. R. Buckland) 


10.4 (—.-) 


Ann. Univ. Sci. Budapest—Sect. Math. (1958) 73-82 (8 refs.) 


Suppose that at a telephone centre calls are arriving in the 
instants 71, T2,.--) Tn)... Lhere are m full-availability channels 
and there is a waiting facility. If all channels are busy 
and the number of waiting calls is less than w (w is the maxi- 
mum number of waiting calls), then each new call joins 
the waiting line and waits until a channel is free. If all m 
channels are busy and w calls are already waiting then any 
new call is lost. 

We say that the system is in state H,(k = 0, I, 2,.... 
m, ..., m+w) if the number of calls being served and waiting 
is k. Define the random variable 7(¢) : 

Ok (BO, Kye. Ml, MALT, MW) 
if at the instant ¢ the system isin state H,. Put 7(t,—0) = %n- 
The problem is to determine the limiting distributions 
* 

lim P{y, = k} = Py, lim P{n(t) = k} = Pp. 

no to 
The author determines these limiting distributions if the 
following assumptions are made : 

(a) The instants 7, form a sequence of recurrent events, 
i.e. the random variables tn41—Tn (m = I, 2, ...) are identi- 
cally distributed and independent. Their common distri- 
bution function has a finite first moment. 

(b) Denote by x» the holding time of the call arriving at 
the instant t,. Suppose that x, are independent random 


TAA TnL 


variables which are independent also of the process {7,}. 
We suppose that the conditional distribution function 
P{Xn Sx | an = k} (Rk = 0, 1, ..., m-+w—1r) is an exponential 
one with the same parameter » for each R. 

The basis of the determination of the limiting distri- 
bution P; is that the sequence », forms a Markoff chain. 
Further, if the distribution function of the random variables 
T41--T is not a lattice distribution then the limiting 
distribution lex is also determined. It is shown in the 
paper that P;, and P;, (Ai=t0; ni, 
of the initial state. 

Finally the distribution function of the waiting time 
under the condition that a call is not lost is determined. 


..., m+w) are independent 


(J. Mogyorddi) 
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On the probability distribution of the measure of the union of 


random sets in a Euclidean space—In Hungarian 


Ann. Univ. Sci. Budapest—Sect. Math. (1958) 89-95 (8 refs.) 


Consider certain random points {Pi} distributed in the 
Euclidean space of m-dimensions. To each random point 
let there be a corresponding Borel-measurable random set 
Hi in this space. Let this random set depend on a random 
parameter «;. It is supposed that these random parameters 
are identically distributed independent random variables. 
Further suppose that the indicator functions of the random 
sets can be expressed by the same Baire-function f (r, a), 
where 7 = (11, 72, ...; 7m) denotes the points of the space 
and a is a positive real number. Namely, if «; = a; and 1; 
denotes the radius-vector of Pi, then we have 


1 if r €Hi 
re Co a i if r ¢Hi. 
Let ¢(x) = 1, if x = (x1, x, .... Xm) is a point of at least one 
random set and ¢(x) = 0 otherwise. If H is a fixed set with 
finite positive Lebesgue-measure in the Euclidean space of 
m-dimensions, then the expression 


@(H) = i hea 


gives the measure of those points of H, which are covered 
by at least one random set. @(H) generally is a random 
variable and the problem is to determine the distribution 
of (A) if different assumptions are made about the dis- 
tribution of the random points. The problem is investigated 
in the following four cases : 

(a) The random points are distributed only on the set H 
according to a Poisson process with density A, i.e. if v(H) 
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denotes the number of random points in H, then 
Prob {o(E) = n} = (1/n!) {Aw(ED}" exp [—Au(HD)] 
where (7) is the Lebesgue-measure of the set H. 

(b) The random points are distributed on the whole 
space R,, according to a Poisson process with density A. 

(c) The random points are distributed only on H. It is 
supposed that »(H)=n and the random points are distri- 
buted independently of each other and uniformly on H. 

(d) The random points are distributed only on H. It is 
supposed that »(H) is an arbitrary random variable which 
assumes integer values and under the condition o(H) =n 
the n random points are distributed independently of each 
other and uniformly on H. 

In the four cases mentioned above the author determines 
the moments of 0(H): 


MLO} =| |, Prob {Eea) = t - Dehn ae 


where 
Prob {£(x1) = 1, «5 (s) = 1} = 2 (—1)! 
k=0 
2 D(X;,, X 559 


; : b Beery %;,) 
Deis fasess SS 

Here the probabilities p(x,, x2, ..., x,) are expressed by 

the indicator functions of random sets and the distribution 

function of the random parameters. ‘The problem is thus 

solved, since the distribution function of O(H) is uniquely 


determined by its moments. (J. Mogyorddi) 


10.3 (10.9) 


Probabilistic treatment of meteoropathological phenomena—ZIJn Hungarian 
Mag. Tud. Akad. Kut. Mat. Intézet. Kézleményei (1958) 129-140 (7 refs., 3 tables) 


Let us suppose that certain climatic phenomena (e.g. 
meteorological frontal passages) occur at the instants 4, 
to, ...» tn, ... and certain biological phenomena occur at the 
instants U,, Ug, ..-) Un, ...- Itis supposed that every meteoro- 
logical front gives rise to an effect extended over an interval 
with length « (the instants ft, are the centres of the corres- 
ponding intervals). One can observe the instants of the 
meteorological fronts and of the biological phenomena in the 
time interval (0, J) and the problem is to decide in what 
degree the biological phenomena can be considered as an 
effect of the meteorological fronts. 

The mathematical model of the problem is the following : 

Consider the case when the time differences t,41—t 
are identically distributed, independent, positive random 
variables with distribution function F(x). We suppose that 
the climatic phenomena have no effect. In this case it is 
reasonable to suppose that the sequence of the biological 
phenomena {un} forms a Poisson process with density » 
and {u,} is independent of the process {tn}. 

Denote by vz the number of those biological phenomena 
in the time interval (0, 7) which fall into the intervals 
corresponding to the meteorological phenomena. Further 
denote by Ve the actually observed number of the biological 


phenomena. Then the probability P{v,<v;} may be con- 
sidered as a measure of the effect of the meteorological 
phenomena. 


rAh 


The author proves that if T is sufficiently large then 
Ply, Svy}~$, —MT/DT?) 


x 
edt and lim M{v;}/T = M 


ioe} T= 0 


where (x) = 1/(27)? | 


and (Dea cree OLN he 


T- 0 


M{v,} and D*{v;} may be determined with the aid of the 
distribution function F(x). 

Finally the problem is treated whether the biological 
phenomena can be considered as an effect of some climatical 
phenomena acting continuously with intensity A(t), 
(o<t<T). 


(J. Mogyorddi) 
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TANNER, J. C. 


: . (Road Research Laboratory, D.S.I.R., England) 
A simplified model for delays in overtaking on a two-lane road—In English 


10.4 (2.5) 


J. R. Statist. Soc. B. (1958) 20, 408-414 (2 refs., 1 table) 


The model considered is of vehicles travelling in both direc- 
tions along a long straight two-lane road. The vehicles 
travel at constant speeds of v in one direction and V in the 
other and are assumed to be spaced at random in their 
streams ; the density of vehicles in the two streams being 
qand Q. The problem investigated is the delay experienced 
by a single additional vehicle during uninterrupted travel in 
the first direction at a constant speed u (uv), but which 
may only pass a v-vehicle according to the following rules : 


(i) If the distance between the v-vehicle and the next 
approaching V-vehicle is d or more then it passes 
without slowing down. 

(ii) If this distance is less than d it slows down instantane- 
ously to v and waits until there is a distance D to the 
next approaching V-vehicle before overtaking. 

(iii) It overtakes v-vehicles one at a time. 


A function is obtained for the average speed u achieved 
by the u-vehicle in the course of a long journey along the 
road. ‘This function is tabulated for a few selected values 
of the parameters. For specified speeds of the vehicles, ex- 
pressions are also obtained for d and D—the minimum gaps 
in the opposing traffic which will allow passing to be safely 
performed. 

The author points out that the model does not allow 
for the bunching which occurs in practice, owing to the 
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URBANIK, K. (Math. Inst., Polish Academy of Science, Warsaw) 


variations in the speeds of vehicles in one direction, but 
suggests that it gives a fair approximation for roads with 
fairly light traffic. He mentions possible extensions to the 
model, to allow for bunching, by considering each stream of 
vehicles to be fed into the road as the output from a fixed— 
service—time-queueing process with random input. He 
points out that, though this type of bunching is different 
from that occurring in practice, the theoretical analysis of 
more realistic bunching is beyond the scope of present 
techniques. 


(N. W. Please) 


10.1 (—.-) 


Remarks on invariant functions in Markoff processes—In Polish 


Collog. Math. (1958) 5, 223-230 (3 refs.) 


The note referred to is one of a series of Urbanik’s papers 
devoted to a study of limiting properties of Markoff processes. 

There are considered homogeneous Markoff processes 
whose sample functions are step functions of continuous 
time ¢ assuming their values in a finite or denumerable set. 
Transition probabilities are assumed to depend continuously 
on time. A real valued function f, defined on the sample 
space of the process, is called invariant if it is measurable 
with respect to the class of sample function sets invariant 
under the shift transformations. ‘The purpose of the note 
is to examine the power of a set of values which can be 
essentially assumed by an invariant function. Let us call 
the function Q(x), defined for values x of the process as a 
limit of a probability that the process takes on x as its value 
at moment f if t tends to infinity, the limiting distribution 
of the process. ‘Two theorems are proved. ‘The first states 
that if the sum of all limits Q(x) extended over all values of 
the process equals unity, then the invariant functions of the 
process assume essentially at most as many values as there 
are elements in the set of values of the process. In other 
words, if c is the power of the set of values of the process, 
then for each invariant function there exists a set of power 
c of its values such that it assumes a value from that set with 
probability one. 

The second theorem states that for each sequence 4q;, 
qa, ... of non-negative real numbers with sum q less than 


unity, there exists a Markoff process with properties de- 
scribed above for which the limits Q(x) assume exactly all 
values of the sequence q;, gq, .... Thus their sum, being 
equal to q, is less than unity; moreover there exists an 
invariant function of this process which assumes essentially 
non-denumerably many values, i.e. for each denumerable 
set of values of that function the probability that it assumes 
a value from this set is less than unity. 


(S. Zubrzycki) 
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WETTE, R. (Zool. Inst., Biomath Dept., University of Heidelberg) 


10.9 (2.5) 


On the biomathematical basis of the logarithmic series distribution of elements % 
of taxonomic units in the natural system—Jn German 


Biom. Zeit. (1959) 1, 44-50 (7 refs.) 


Following a short discussion of Kendall’s and Skellam’s 
models for the generation of the empirically observed loga- 
rithmic series distribution of species per genus by stochastic 
processes in species evolution, a different model is suggested : 
the allocation of species to taxonomic genera by taxonomists 
is a stochastic process continuous in time with the following 
properties : 
(1) The probability of allocating an additional species 
to an existing taxonomic genus is proportional to the 
size of the latter (probability of average likeliness), 


(2) The probability of withdrawal of one species from a 
genus is proportional to the size of the latter (number 
of possible choices), with the exception that 


(3) Unspecific genera (of size unity) tend to persist or 
when being abolished change their name with finite 
probability. 


This is a modified Feller-Arley ‘“‘ birth-and-death ”’ 
process, the modification consisting of a reflecting barrier 
between zero and unity. The author avoids integration of 
the derived partial differential-functional equation for the 
probability generating function by considering only the 
limiting stationary distribution. This stationary distribution 
is a logarithmic series with a zero term, the latter being 
censored in the empirical distribution. 


(R. Wette) 
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ZITEK, F. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 


Random Functions of an interval—In French 
Czech. Math. fF. (1958) 8, 583-609 


The author considers random functions of an interval 
and their Burkill integral. In the paper it is supposed that 
random variables associated with any two adjacent intervals 
(i.e. two disjunct intervals the union of which is an interval) 
are stochastically independent. A random function of 
an interval is called additive if it associates with the union 
of two adjacent intervals the sum of random variables 
associated with separate intervals. The author introduces 
a notion of continuity in the following sense: the random 
variables associated with a sequence of intervals the length 
of which converges to zero converge to zero in probability. 
The Burkill integral is defined as follows: we divide a 
given interval into a finite number of subintervals, add the 
corresponding random variables and let the maximum 
length of subintervals tend to zero. If the corresponding 
sums tend in distribution or in probability to a certain limit, 
we say the random function has a (BB)-integral or a (pB)- 
integral, respectively, over the interval considered. The 
function is called integrable, if it has an integral over all 
subintervals of the given interval. Then the indefinite 
Burkill integral is an additive random function. Subse- 
quently, the author gives some criteria for the existence of 
Burkill integrals of random functions, among them a criterion 
for the case of continuous random functions which is based 
on the theory of limit laws of probability. 


10.0 (1.6) 


The notion of (BB)-integral is used to make more clear 
some aspects of the theory of stochastic differential equations 
in the interpretation of P. Lévy. The author starts with the 
fact that a differential equation determines principally the 
distribution law of increments of a random function to be 
obtained. These increments can be written in the form of a 
random function of interval. By integration of this function 
we get an additive random function which is the solution 
of the differential equation considered. 


(M. Sova) 


Remark of the Regional Editor : 


The results of the paper are partly overlapping with 
recent results of: A. Prékopa (1957), On stochastic set 
functions, III [Acta Math. Acad. Scient. Hung. 8, 375-400]. 

(A. Rényi) 


BARNARD, G. A. (Imperial College, London) 
Thomas Bayes.—A Biographical Note—In English 
Biometrika (1958), 45, 293-295 


Biographical information regarding Thomas Bayes, F.R.S. 
(1702-1761), the father of inverse probability, is very scanty. 
The author has collected together such information as there 
is and presented it as a short three-page preface to the 
reproduction of Bayes’ essay presented to the Royal Society 
by Richard Price in 1763. Bayes, the son of a Presbyterian 
minister, was educated privately and became himself a 
Presbyterian minister in Tunbridge Wells. It is not known 
when he went to Tunbridge Wells, he was there in 17.20 
and he died there in 1761, having retired from his ministry 
in 1752. He was buried with his parents and brothers and 
sisters in Bunhill Fields. 


(F. N. David) 
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FOSTER, F. G. (London School of Economics) 


11.9 (1.0) 


Ir.t (5.8) 


Upper percentage points of the generalised Beta Distribution, II[I—IJn English 


Biometrika (1958) 45, 492-503 ( refs., 1 table) 


The problem discussed in this and two previous papers by 
the author is that of testing a null hypothesis for k-variate 
normal populations by considering two independent estimates 
A and B, of what, on a null hypothesis, is the same dispersion 
matrix. The roots X of the determinantal equation of 
B—)A are invariant under all linear transformations of the 
Rk variates, and it has been argued that some functions of 
the roots might be a useful criterion for testing the null 
hypothesis. 

The author has been concerned with the tabulation of 
the greatest root A, or, what is shown to be equivalent, of 
the greatest root 6 of the determinantal equation of the 
difference between the sums of products matrices of B and 
A. 06, a function of X and the degrees of freedom of A and 
B, increases monotonically with 4, and all @ lie between o 
and 1. The distribution function of the greatest 0 depends 
on k, p and g (where p and q are respectively linear functions 
of the degrees of freedom of A and B) and is given the name 


of ‘‘ generalised Beta Distribution ”’ since it reduces to the . 


Beta Distribution when & is 1. 

The 80, 85, 90, 95 and 99 per cent. points were tabulated 
with k having value 2 in the first paper (published jointly 
with D. H. Rees). The value of k in the second was 3, and 
4 in the tabulations of the present paper. The degree of 
freedom of A and B take values 5-195, 4-11 respectively. 


Computations were carried out on the English Electric 
Company’s DEUCE, as they were for the second paper, to 
which reference is made for details of the programme and 
interpolation requirements. Modifications of the computing 
procedure for k equal to 4 are given in this present paper. 
The existing programme could now be used to extend any 
of these three tables to higher values of the degree of freedom 
and hence of p and qg. However, it is probable that this 
method would cease to be useful when k& exceeds, or even 
equals 5, and a new theoretical approach appears necessary. 


(D. Allibone) 
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HOTELLING, H. (University of North Carolina) 


11.9 (1.0) 


The statistical method and the philosophy of science—In English 


Amer. Statistician (1958) 12, 9-15. 


This is a philosophical discussion of probability and statistics 
which was presented as a presidential address to the Uni- 
versity of North Carolina chapter of Sigma Xi. Accordingly 
it is directed at an audience of general scientific (rather than 
statistical) background. The aim of the paper is to indicate 
the impact that recent statistical theories are likely to have on 
the philosophical disciplines of metaphysics, epistemology 
and inductive logic. 

The distinction between subjective and objective ideas of 
probability is discussed with particular reference to the 
theories of Keynes and Carnap. The axiomatic approach 
of Kolmogoroff is recalled. ‘The author suggests that it 
may be appropriate to take ‘“‘ probability’ itself as an 
undefined term, since every formal theory must necessarily 
contain certain undefined concepts. It is pointed out that 
the representation of ignorance by uniform probability 
distribution over a finite number of points may sometimes 
have justification via ergodic theory ; as an example, in the 
shuffling of a pack of cards a uniform distribution of permu- 
tations is approached as the shuffling time is increased. 

The philosophy of the usual techniques of point and 
interval estimation and of testing hypotheses is discussed. 
The impact of Fisher’s Design of Experiments and of recent 
developments in sequential experimentation are recalled. 
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Finally it is pointed out that there is a parallel between the 
pragmatic philosophy of William James and the theory of 
decision functions : One’s stated objective is not to discover 
a true value, but rather to minimise a loss. 


(R. J. Buehler) 


11.6 (0.1) 


LivavaTI—A new analogue computer for solving linear simultaneous equations and related problems. 


Part II. 


Design of Model II and its application to the solution of secular equations—In English 


Proc. Ind. Acad. Sct. (1958) 48, 269-283 (3 refs., 8 figs.) 


The general principles involved in the design of LILavaTI 
were discussed in a previous paper (in the same volume of 
the same journal) by the author jointly with G. N. Rama- 
chandran. This second paper deals with an improved 
model of the first computer. The new design employs 
potentiometers instead of the meters in the Model I and so 
has the following advantages over it: (a) considerable 
economy is possible in material and in number of operations 
(b) the circuit can easily be modified to solve secular equa- 
tions (c) operations can be made with a.c. and d.c. supply 
without changing the circuit appreciably. 

Another advantage is that in this model, values of currents 
need not be read at intermediate stages but can be measured 
when the iterative process has converged. The technique 
of measurement is based on a null setting so that it is capable 
of great accuracy in principle. The potentiometric arrange- 
ments and the setting of the circuit to solve n linear simultane- 
ous equations in m unknowns are explained. Details of the 
model are given with the help of diagrams. 

Modifications to be made in the circuit to extract the 
eigen values and eigen vectors of a matrix A, when the 
iterative method of Frost and T'amres is used, are considered. 
The computer described in this article can handle only a 
third order secular equations; but the author points out 


that there is no essential difficulty in constructing a larger 
computer based on these principles by making use of switches 
having longer capacity. Improvements in the model and 
the method of solving the equation AX=ABX will be dis- 
cussed in Part III (to be published). 


(P. K. Kumaran) 
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MENDENHALL, W. 


(London School of Economics and North Carolina State Coll.) 


11.9 (2.0) 


A bibliography on life testing and related topics—In English 


Biometrika (1958) 45, 521-543 (approx. 620 refs.) 


This select bibliography covers papers concerned with 
statistical theory and methods applicable to the study of 
the life characteristic of some biological or physical body. 

Although many of the statistical methods developed for 
life tests of industrial products are applicable to studies of 
life in other fields, and vice versa, major emphasis is given 
to the industrial applications. 

The special conditions associated with the life test and 
analysis of its observations are briefly considered: in 
particular the occurrence of ordered observations from 
censored samples with or without truncation of the popula- 
tion—and the accelerated test. The introductory matter to 
this select bibliography notes that very little has been 
published in connection with this last topic. Other topics 
of considerable interest in the study of the life characteristic 
are: extreme value theory, fatigue and wear tests, system 
reliability and the treatment of sensitivity data. The refer- 
ences have been classified as belonging to one or more of 
nine groups, as follows : 

C. Censored sampling, or sampling from a truncated 

distribution—both univariate and multivariate. 


R. Sampling with renewal. 

O. Order statistics. (Excluding papers on extreme value 
theory or censored sampling.) 

E. Extreme value theory. 

K. Papers concerned with failure rates and conditional 
failure density. Also included are tests of random- 
ness for a sequence of events occurring in time. 
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F. Fatigue testing and wear problems. 

M. Machine productivity problems. This group in- 
cludes some representative papers concerned with a 
group of machines subject to various failure laws and 
servicing arrangements. 

S. System reliability. 

D. Methods applicable to sensitivity data: including 
the fitting of dosage-mortality curves. 

In the bibliography the references are presented in an 
alphabetical order by authors with the appropriate group 
classification (indicated by the above capital letters) shown 
in parentheses at the end of each reference. Some of the 
groups are obviously subtopics of others but the system of 
classification has been kept relatively simple, e.g., references 
under censored sampling are not represented under order 
statistics. Some of the groups are represented only by some 
of the more important papers which may be consulted for 
additional references. In general, the main body of con- 
tinental literature is represented only by papers revealed by 
the references given at the end of other papers. 


(Wm. R. Buckland) 


11.5 (6.1) 


Calculation of multiple correlations by means of electronic computers—In German 


Metrika (1958) 1, 148-153 (1 ref.) 


The paper deals with a method of calculating regression 
coefficients and coefficients of determination by means of an 
electronic computer (IBM 604). The paper discusses the 
calculation of the covariance matrix and indicates the number 
of required products. For the solution of the system of 
linear equations necessary to obtain the regression co- 
efficients the method of Gauss is used. (The method of 
Gauss-Doolittle, which is not mentioned by the author, 
would allow for a remarkable reduction of the computa- 
tion-work by making use of the symmetry of the covariance 
matrix.) The paper further indicates a method of calculating 
the coefficients of determination by means of solving a 
system of linear equations. The author shows that the 
amount of computation is smaller than in the case of 
using direct calculation. For a more detailed discussion 
of the method of Gauss-Doolittle and the use of punchcard 
equipment the reader is referred to ‘Practice of multi- 
variate correlation calculations”, Reihe Mathematische 
Statistik, edited by Institut fiir Praktische Mathematik, 
Prof. Dr. A. Walther, Technische Hochschule, Darmstadt, 
Parts’, lland IN. 
(J. Roppert) 
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SIOTANI, M. & OZAWA, M. (Institute of Statistical Mathematics, Tokyo) 


11.1 (5.4) 


Tables for testing the homogeneity of k independent binomial 
experiments on certain event based on the range—In English 
Ann. Inst. Statist. Math., Tokyo (1958) 10, 47-63 (1 ref., 15 tables) 


This paper gives tables for testing the homogeneity of 
independent binomial experiments. The test criterion is 
based upon the range of the observed number of occurrences 
v;in the k experiments of N trials. In this paper is considered 
the case where N is not large enough to apply the Normal 
approximation to the single succession of trials but where, 
under the null hypothesis: H:p,=p. = ... py=p (say), Nk 
is large enough to regard the estimate $ = L»;/Nk as the 
true value of p. 

Let v; (¢ = 1, 2, ..., k) be the observed number of occur- 
rences, and p;(¢=1, 2,..., k) the probability of a certain 
event in the 7th binomial experiment of N trials. The test 
criterion is defined as: R,(N,p) = vy—v2 where », and v2 
are the largest and smallest values, respectively, in the set 
of k numbers 1, v2, ..., vz and p is the probability of the event 
under the null hypothesis. 

The tables in the paper were computed on an auto- 
matic relay computer for N = 10(1)20, 22, 25, 27 and 30; 
k = 2(1)15; p= 0°50 (0°10) o:10. ‘The entries in the 
tables, while they are related to the probabilities 


P{R,(N, p) 2rz} 


are, in fact, the greatest value of 7; such that the probability 
P{R,(N, p) 21x} does not exceed ten specimen values of the 
level of significance (a) for fixed N, p and k. The tabulated 
values of « are 0'001, 0:005, 0-or (0°01) 0:06, 0°08, o-10. 
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The reason for adopting this device of tabulation is solely 
on account of clarity and economy of presentation. Attention 
is drawn to certain problems of rounding common with 
tables of this kind and a numerical illustration is given. 


(C. Hayashi) 


11.8 (-.-) 


Methods Used in Compiling the United Nations Price Indexes 
for Basic Commodities in International Trade—In English 


Statist. Papers ; Series M, No. 29 (1958) (66 pp.) 


In the issue for the last month of each quarter of its Monthly 
Bulletin of Statistics, the Statistical Office regularly publishes 
indexes showing different aspects of the price movements 
of basic commodities in international trade. Half the weight 
in certain of these is assigned to export prices (f.o.b.) and 
half to import prices (c.i.f.) in order to take account of the 
fact that the price determining markets for a number of 
primary commodities are in the importing countries and 
that, for them, export rather than import prices are subject 
to variation caused by change in freight rates and the fact 
that, for a number of important currents of trade, import 
prices are available and export prices are not. The con- 
struction of these indexes depends on the methodological 
innovation of dividing weights evenly within each commodity 
between export and import prices. This is accomplished by 
appropriately dividing between the export and import prices 
the weight of each current of trade for which both export 
and import prices are available. The manner of making 
the division is described in general and also specifically 
for each of the 47 primary commodities entering the indexes 


in question. The paper furthermore lists the price series 
entering the indexes, gives the weight of each and sets out, 
for each commodity, the economic considerations that 
determined the weighting pattern. 

French and Spanish editions are in preparation. 


(W. R. Leonard) 


a” oe 

uni@oly) wi) aoa. fe > eroremenel aM) antes to aiday. 4 
“Wi ‘ead ‘ire n ) pork ner ST wallets an Stam 
. ee diire o oreeihi aj ovseuccnndsted cof ne’ boreal ig 
r| ‘ tai Exe otefy Wo wahd) iver) earn o L .sqnrg alate od, are YW a8 

werroPS male elgg ot bgpsereais piahi pede 

~ranlw joe deay We COLSON fr he ott, o 


ol) AE Mh ey a, glares 


; vet AVA\ yi, em Fe eetaligviien oti: pein ot a 


Hin te select bavpoedo ont 90 18 
mah ®. te Lait lihadieere, arit iAus iP rete 
toa) of lL’ lamer Vie bedvces leinvewtid, 
ee ‘bee im YtSsjee en ¢ - “VA et 
H ‘lovee pier seb flan by 
: ; titey6-oeft | wrilide Liven WS 8 my Ratt hc, ally 

2 mmantioned I 

f -Ciik He f rei shake yew i | sult ab. 
Of bre Se 28 te other <= Vitel mene 02. 

: . oct mi eoriisrs ar. 76. tore) ao =e ry 4 
i . ; worhidadioag walycn taemiiss. core, Cond) lle 


: | det VO 


ywildedorq ay tech ron & qt to oulagr rewherp att a! 
i 10 Bayley ris tee | h es Hoes 20 7 t =e" 
baniludat aft toe & VM = = (ata ave 


“er 


OF Boao oa » d20°8) 1B » hehe am @ 
. = : 


i 
ii 


% 


se fnotytideal?) SVIOTTAM GB iTitAU- SH A aTRO 


l sind anoint Vt besbialy uf} geiliqmo;? ni lesutd ah he 
Wwilgt\ wl—sbarl’ lanoiterreetal ob aghiit Lamune peek qe i 
(4 68) > (Re ei) Ge a4 MV ono 7 ewe tek nae 


_ 
a, ‘ ; ; (att) Oe > sr} a fi vidinoba, est Fo Witty if yr its es ey ef titel = 
BE Ais bree Howe Yo) ; ey ersine wit yale  seufeidioy ehinloyerstits teste ly aperen 
tras tf ; ibovreerins wot ireutovort GM ol) Bev eee inert eke 
— iad imysteh tpiow sf UH udder To rrodaucresti§ tie -TTieoernct 
7 " / lame bern o ~ int Cato’) 27 ogee a lachemngiaames a onart 79 5 vite 
¢ . 1) eT Salen oe valve nb CAMS) seating, ste i 
; io cadena aol ateleet yah ene 
Piha P : : ie ate oe | yolWwiey hat het epee s(t ith ‘wt ron f 
a pelching ow by 9 ton ik Seely ‘rite ots 
: vet “i? Dk ested Urat of aggeie gah 
ROmnre ahurks Me ener trues Fora) yt 
no alt cern Tle sag TG a ‘By 
’ . lesiwolo beets lt an cogs ae sie 
: etibornuns oahials iive a rit ei Won 
+e tHolall aie GTS ey AL Ley t 
eating ® ee [eek Leen Proeyey inal 
athe » dbrewt) sloetae ge wiais bee 
nite i ie i a ‘ 


colroegs dale on fe Frontal 


rob od? yoked 


VAN YZEREN, J. (Netherlands Central Bureau of Statistics) 


11.8 (11.3) 


A note on the useful properties of Stuvel’s index numbers—In English 


Econometrica (1958) 26, 429-439 (1 nomogram) 


This paper develops the algebra of the formulae introduced 
by G. Stuvel [Econometrica (1957) 25, 123-132] for price 
and quantity index numbers. The principal feature of these 
indices is the synthesis of the additive and multiplicative 
way of decomposing value changes. Whereas Laspeyres’ 
indices generally have an upward bias, in the Stuvel indices 
the biases cancel out. Stuvel’s indices, unlike Laspeyres’ 
and Paasche’s, satisfy the time reversal test. 

The properties of Fisher, Edgeworth and Stuvel indices 
are considered so far as the aggregation of the two sub- 
aggregates into one index is concerned. Whereas, with the 
first two types of index number, the aggregation cannot be 
carried out using value data only, it is shown that there is a 
simple and symmetric way of aggregating Stuvel indices 
with the help of value weights. A further point is that if 
two Stuvel indices are equal, then the overall index will have 
the same value, whereas both the Fisher and Edgeworth 
formulae violate this requirement. 

It is shown that Stuvel’s indices are unique in that similar 
indices derived from Paasche, and not Laspeyres’ formulae, 
do not have the same derivable aggregative properties. 

Algebraic formulae are then developed in order to com- 
pare the values taken by Stuvel indices with those taken by 
Laspeyres, Paasche, Fisher and Edgeworth and Stuvel 
indices. Under given conditions, Fisher, Edgeworth and 
Stuvel indices are intermediate between those of Laspeyres 
and Paasche. 
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VINSKI, IVO. (Ekonomski Institut, Zagreb) 
Problems in Calculating National Wealth—In German 
Metrika (1958) 1, 158-172 (13 refs., 3 tables) 


After a historical introduction (I) the article discusses con- 
ceptual problems (II), the relationship between national 
income and national wealth (III), the rate of growth of 
national wealth (IV), and international comparisons (V). 

Since estimates of national wealth were in vogue during 
the eighteenth and nineteenth centuries no further methodo- 
logical progress has been made since the earliest thirties. 
But with the rapid progress of modern national accounting 
after World War II, national wealth again got attention in 
the professional world. 

National wealth contains natural resources and produced 
goods. Of the former, only cultivated soil is included since 
unused natural forces etc. cannot be evaluated. Of the 
latter, museums, monuments, materials for the armed forces, 
are excluded for want of an economic evaluation. Differ- 
ences between the theorists are much greater than between 
statisticians engaged in those calculations, since the statistical 
material available in most countries reduces too ambitious 
concepts. Formerly, economists preferred an indirect 
method of evaluation, i.e. by capitalisation of income flows. 
The author regards this method as obsolete and inaccurate. 
Apart from that it is not applicable if the greater part of 
wealth is not in private hands. Modern procedures are 
based on direct evaluations by surveys either in form of 
complete enumerations (census of agriculture, industry, 
trade, public utilities, etc.) or in the form of sample surveys 
(dwellings, furniture, stocks in households, etc.). If such 
surveys are not simultaneous, the results can be synchronised 


The position of Edgeworth and Fisher price indices with 
respect to Stuvel’s depends on the absolute value of the 
Edgeworth index and on the relative values of the cor- 
responding Laspeyres’ price and quantity indices. If there 
is no price change according to Edgeworth then there is 
none according to Stuvel either. 

An empirical test with some numerical data is carried 
out, which suggests that Fisher and Stuvel indices meet 
the circular test equally well. 

Finally, the proportionality test is discussed, which the 
Stuvel indices do not satisfy. A way out of this difficulty 
is suggested which involves the modification of the formulae 
proposed by Stuvel by the introduction of an additional 
parameter, which would normally take the value 1, but 
could, if (say) a monetary reform were carried through, be 
adjusted to allow for it. The introduction of this extra 
parameter does not destroy the other desirable properties 
of the Stuvel indices. 


(D. J. Day) 


11.0 (8.0) 


by Goldsmith’s “‘ perpetual inventory method ’”’, i.e. by 
adding net investment to capital stocks of back years. The 
crucial problem of all national wealth estimates is the 
principle of evaluation itself. Of the three possibilities : 
(1) original purchase values, (2) market prices (i.e. sale 
values) and (3) reproduction values (gross and net of 
depreciation), the author argues for the latter and presents 
—with due regard for the unreliability of depreciation—all 
his figures for Yugoslavia in two series (gross and net). 

The relationship of national income to national wealth 
seems not to be a plain one. It is not simply an indicator of 
capital productivity, but depends on natural conditions, 
labour force, terms of trade, etc. Different stages in eco- 
nomic development, therefore, produce relationships of 
different size. In underdeveloped as well as in highly 
developed countries the figure is relatively small (national 
wealth equals about two to three times the national income), 
in countries with forced industrial development (as Yugo- 
slavia) the figure is about four to five. In a similar way the 
rate of growth of national wealth depends on the stage of 
growth already arrived at: it seems to flatten out with the 
increasing growth and age of the economy. 

With regard to international comparisons the author 
recommends not to rely on per capita figures expressed in 
arbitrary exchange units but to compare the structures of 
country totals of national wealth. 'The comparison in the 
author’s table III is convincing. 


(L. Bosse) 
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